Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,603 A→G 6.2% A17A (GCA→GCG)  tisB → toxic membrane persister formation peptide, LexA‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,6030AG 6.2% 599.0 / 20.7 186A17A (GCA→GCG) tisBtoxic membrane persister formation peptide, LexA‑regulated
Reads supporting (aligned to +/- strand):  ref base A (82/91);  new base G (5/7);  total (88/98)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 194 total aligned reads displayed.

TCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGG  >  NC_000913/3853332‑3853873
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
tCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGCCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTTGGTATCGCCATTCTTATCCCCAAACTCATTGTTCCAGCACt                                                                                                                                                                                                                                                                              >  4:469298/1‑274 (MQ=255)
      tACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCCCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCAAAGGAGAAGCGTATGAACCT‑GGTGGATAACGCCATTCTTATCCTCAAACTCATTGTTGCAGACCTTCAACTg                                                                                                                                                                                                                                                                       >  4:434242/1‑275 (MQ=255)
             tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACCGGCAGCGGCAAGGTACGTCCGCTGGCCATTCCCCTGAACCACAGGCACCGCATATGACCCTGGGTGATTATC‑CCATTCTTACCCTCAAACTCATTGTTGAAACACCGCCACTTCTTggt                                                                                                                                                                                                                                                                 >  4:6865/1‑272 (MQ=255)
             tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  4:334687/275‑1 (MQ=255)
             gtggtgaTGCGGGGGCCGCATCGGTACCCGGCGCACTAAGTCCTGTGTGAACCGGGTGGGGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  4:485872/269‑1 (MQ=255)
                         gggCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAAt                                                                                                                                                                                                                                                    >  2:449739/1‑275 (MQ=255)
                                            cgcACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAATTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGAtt                                                                                                                                                                                                                                 >  1:123172/1‑275 (MQ=255)
                                                           cTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACCGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCAATTCTTATCCCAAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTTAAATACCTGAAGTAATTCACATTCAAGTCGTCCCaaa                                                                                                                                                                                                                   >  4:530958/1‑274 (MQ=255)
                                                           cTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTc                                                                                                                                                                                                                                >  4:173913/1‑261 (MQ=255)
                                                                          ccGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCg                                                                                                                                                                                                   >  2:341843/1‑275 (MQ=255)
                                                                             tCAGCGCCTTAACCCCGCTTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                                                                                <  2:202899/275‑1 (MQ=255)
                                                                               aGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTTCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTccc                                                                                                                                                                                              <  1:72617/275‑1 (MQ=255)
                                                                                              cgTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCt                                                                                                                                                                                                                                                >  2:215110/1‑210 (MQ=255)
                                                                                              cgTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGACATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGACTCAAGTCGCACCAAAGGGGAGCGGGAAACCGATCCCCTTTTATATTTAGCg                                                                                                                                                                               >  2:514641/1‑275 (MQ=255)
                                                                                                     acacTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGt                                                                                                                                                                        <  4:341842/275‑1 (MQ=255)
                                                                                                                     acaaGCACAGCGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACCCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcgat                                                                                                                                                        >  4:120326/1‑275 (MQ=255)
                                                                                                                          cccaACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGAAACAGGCAGCGGAAAGGTACGTCACCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGa                                                                                                                                                   <  2:485977/273‑1 (MQ=255)
                                                                                                                                     tcCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       <  6:183385/148‑1 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCtt                                                                                                                                                                                                                                                                    >  8:1048993/1‑148 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCtt                                                                                                                                                                                                                                                                    >  7:752478/1‑148 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCCCAAACTCATTGTTGCAGCACTGCa                                                                                                                                                                                                                                                                           <  6:1019729/141‑1 (MQ=255)
                                                                                                                                         tGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   <  6:77958/148‑1 (MQ=255)
                                                                                                                                              tACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGa                                                                                                                                                                                                                                                                  <  6:137894/144‑1 (MQ=255)
                                                                                                                                                 cgcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGtt                                                                                                                                                                                                                                                           >  8:241473/1‑148 (MQ=255)
                                                                                                                                                 cgcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtc                                                                                                                            >  2:481462/1‑275 (MQ=255)
                                                                                                                                                  gcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          <  6:948558/148‑1 (MQ=255)
                                                                                                                                                  gcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          <  7:247774/148‑1 (MQ=255)
                                                                                                                                                      gCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  5:951016/148‑1 (MQ=255)
                                                                                                                                                       cacaACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTCGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          >  8:591770/1‑143 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGt                                                                                                                                                                                                                                                            >  8:633982/1‑139 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  6:247217/145‑1 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      >  5:247217/1‑145 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAATGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   >  5:1013585/1‑131 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGACCCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGCTGCTGACGCATCTTATCCGCCCTTACATCTCTccctg                                                                                                                    >  1:113985/1‑272 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTC‑TATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  7:87175/148‑1 (MQ=255)
                                                                                                                                                            caGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATa                                                                                                                                                                                                                                                   <  7:1071375/145‑1 (MQ=255)
                                                                                                                                                                  gCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTa                                                                                                           >  2:62666/1‑275 (MQ=255)
                                                                                                                                                                     aCAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGAt                                                                                                                                                                                                                                                                 <  8:845500/122‑1 (MQ=255)
                                                                                                                                                                     aCAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGAt                                                                                                                                                                                                                                                                 >  7:845500/1‑122 (MQ=255)
                                                                                                                                                                      cAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTc                                                                                                                                                                                                                                      <  8:1008035/148‑1 (MQ=255)
                                                                                                                                                                          cgggtgcaaGGTACGTCAGATGGCATTGGTCCTGCGCCCCCGGCGACGCGTATGAACCG‑GGGGGTTTTCGCCTTTTTNTTCCTCCAGTTCTTTGTTGCAGCGTTGCAACTGCTTGAAGCTGTTCTCAAAGACTTGAAATCATTCAGATTCAAGTCGCTCCAAATGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCTGGTTGTTGGCGGATGGGATGCTGACGCATCTTATCCGCCGTACCATCTCTCCCGGCAACAGTTATTGCCGCt                                                                                                   <  3:113985/268‑1 (MQ=255)
                                                                                                                                                                          cAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAcaa                                                                                                                                                                                                                                                                       >  5:936336/1‑109 (MQ=255)
                                                                                                                                                                          cAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAcaa                                                                                                                                                                                                                                                                       >  5:791543/1‑109 (MQ=255)
                                                                                                                                                                             cGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCa                                                                                                                                                                                                                               <  8:176351/148‑1 (MQ=255)
                                                                                                                                                                             cGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAAcc                                                                                                                                                                                                    >  4:350806/1‑175 (MQ=255)
                                                                                                                                                                               gAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       >  7:606488/1‑106 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAca                                                                                                                                                                                                                                                                        <  5:570895/104‑2 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGa                                                                                                                                                                                                                                              <  7:939242/130‑1 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTg                                                                                                                                                                             <  2:158239/195‑1 (MQ=255)
                                                                                                                                                                                 aaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                                                                                <  2:350807/175‑1 (MQ=255)
                                                                                                                                                                                     tACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTaa                                                                                                                                                                                                                                         >  5:586536/1‑130 (MQ=255)
                                                                                                                                                                                         tCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  6:520372/113‑1 (MQ=255)
                                                                                                                                                                                          cAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTCTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  <  5:216523/148‑1 (MQ=255)
                                                                                                                                                                                          cAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:216523/1‑148 (MQ=255)
                                                                                                                                                                                              tGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCa                                                                                                                                                                                                                                                                           >  6:932560/1‑87 (MQ=255)
                                                                                                                                                                                              tGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCa                                                                                                                                                                                                                                                                           <  5:932560/87‑1 (MQ=255)
                                                                                                                                                                                                   gTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTc                                                                                                                                                                                                                                      >  6:835019/1‑119 (MQ=255)
                                                                                                                                                                                                      cTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAc                                                                                                                                                                                                                                                  <  7:690359/104‑1 (MQ=255)
                                                                                                                                                                                                       tCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAATTACTTCAggta                                                                                                                                                                                                                                 <  7:588341/120‑4 (MQ=255)
                                                                                                                                                                                                        ccTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAAt                                                                                                                                                                                                                                        <  5:1035303/112‑1 (MQ=255)
                                                                                                                                                                                                        ccTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      <  7:1057299/146‑1 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  7:53660/87‑1 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                >  8:53660/1‑87 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaaa                                                                                                                                                                                                                   <  8:1057775/132‑1 (MQ=255)
                                                                                                                                                                                                            aaCCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTccc                                                                  <  2:335783/274‑1 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaaa                                                                                                                                                                                                                                                     >  4:104607/1‑91 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 >  5:921547/1‑95 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  6:921547/95‑1 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTa                                                                                                                                                                                                                                          >  8:718090/1‑102 (MQ=255)
                                                                                                                                                                                                                    gagaCGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  2:104608/91‑1 (MQ=255)
                                                                                                                                                                                                                        cgcgTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTaa                                                                                                                                                                                                                                         >  6:498663/1‑95 (MQ=255)
                                                                                                                                                                                                                            gtcgtaaccgtggggATATCTCCATTTTTATCCTCAAACTACTGGTTGCAGCACTGGAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGAAGCAAGTCGCACCAAATGGGAGCGGGAAACCGCTCCCTTTGTATATTTATCGTGCGGGTTGGTTTCGGATGCGATGCTGATGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTatat                                                  <  4:514640/262‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       >  8:215043/1‑60 (MQ=255)
                                                                                                                                                                                                                             aTGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  5:511958/77‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGa                                                                                                                                                                                                                                              >  5:484220/1‑85 (MQ=255)
                                                                                                                                                                                                                             aTGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGc                                                                                                                                                                                                                         <  5:1108677/106‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    <  7:228972/111‑1 (MQ=255)
                                                                                                                                                                                                                                aaCCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      <  5:225271/122‑1 (MQ=255)
                                                                                                                                                                                                                                aaCCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      >  6:225271/1‑122 (MQ=255)
                                                                                                                                                                                                                                  ccT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  5:140715/109‑1 (MQ=255)
                                                                                                                                                                                                                                    t‑ggtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAg                                                                                                                                                                                                                                    >  6:1094208/1‑88 (MQ=255)
                                                                                                                                                                                                                                      ggtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   <  7:174064/56‑1 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTgcagca                                                                                                                                                                                                                                                                                <  4:362530/39‑1 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCGCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 >  6:363591/1‑102 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCGCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  5:363591/102‑1 (MQ=255)
                                                                                                                                                                                                                                           atatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:63982/1‑100 (MQ=255)
                                                                                                                                                                                                                                            tatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCa                                                                                                                              <  3:69893/183‑1 (MQ=255)
                                                                                                                                                                                                                                              tCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAAATGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGAGGGAAAc                                                                                                                                                                                                     >  6:576780/1‑110 (MQ=255)
                                                                                                                                                                                                                                               cGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  8:1096250/97‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:767855/1‑86 (MQ=255)
                                                                                                                                                                                                                                                         aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGAtt                     >  4:498928/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    >  6:450697/1‑79 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    <  5:450697/79‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGc                                                                                                                                                                                                  <  8:88942/97‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTCACATTCATCAATGTGATAGATTCCtctc               >  2:426240/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 attCACTGTTGAAGCACGTCAACTGCTTGAGGCTGTACTGAAATACCTGAAGTAATGTAGCTTCAAGTCGCACCAGAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTGGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGATAAACAATTCCCCATGCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTccc             <  4:13099/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cagcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATggg  <  2:120327/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 >  5:76313/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                               caCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAg                                                                                                                                                                                                             <  7:204648/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCg                                                                                                                                                                                                           >  7:1041858/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
TCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCT‑GGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGG  >  NC_000913/3853332‑3853873

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑