Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,159 G→T 7.9% intergenic (+70/‑65) lrp → / → ftsK leucine‑responsive global transcriptional regulator/DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,1590GT 7.9% 289.6 / 11.1 101intergenic (+70/‑65)lrp/ftsKleucine‑responsive global transcriptional regulator/DNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base G (49/44);  new base T (2/6);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 106 total aligned reads displayed.

CTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCC  >  NC_000913/933010‑933308
                                                                                                                                                     |                                                                                                                                                     
cTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg                                                                                                                                                       >  1:584609/1‑150 (MQ=255)
cTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGAAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg                                                                                                                                                       >  1:1122891/1‑150 (MQ=255)
  ggCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTa                                                                                                                                                     <  1:1673103/150‑1 (MQ=255)
     gTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAcc                                                                                                                                                  >  2:1332134/1‑150 (MQ=255)
      tCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCt                                                                                                                                                 <  1:1674393/150‑1 (MQ=255)
       cAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTg                                                                                                                                                <  1:475060/150‑1 (MQ=255)
        aaTGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTgg                                                                                                                                               <  1:898390/150‑1 (MQ=255)
        aaTGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTgg                                                                                                                                               <  2:507260/150‑1 (MQ=255)
          tGACACACGGACATACTTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTa                                                                                                                                             <  2:418677/150‑1 (MQ=255)
              acacGGTCTTACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAAGACGGATCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTTATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgt                                                                                                                                         <  2:1152388/150‑1 (MQ=255)
               cacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGATCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgtt                                                                                                                                        <  1:1178159/150‑1 (MQ=255)
               cacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgtt                                                                                                                                        <  2:728143/150‑1 (MQ=255)
               cacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgtt                                                                                                                                        <  1:6680/150‑1 (MQ=255)
               cacGGACATACGTTGCTATGGAAGAAGTCAAGCAGAGGAATCGTCTGGTTATTAAGACGCGCTACCACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACtgtt                                                                                                                                        <  2:1321897/150‑1 (MQ=255)
                acGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgttg                                                                                                                                       <  1:521763/150‑1 (MQ=255)
                  ggACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTc                                                                                                                                     >  1:541624/1‑150 (MQ=255)
                  ggACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTc                                                                                                                                     <  1:10819/150‑1 (MQ=255)
                       tACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGttt                                                                                                                                >  2:297990/1‑150 (MQ=255)
                       tACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgttgctgtatg                                                                                                                                <  2:1263801/150‑8 (MQ=38)
                          gttgttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGTTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAg                                                                                                                             <  2:137203/150‑1 (MQ=255)
                           ttgttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGc                                                                                                                            >  2:190626/1‑150 (MQ=255)
                             gttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCAt                                                                                                                          >  1:417967/1‑150 (MQ=255)
                              ttATGGACGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATc                                                                                                                         <  1:1401770/150‑1 (MQ=255)
                              ttATGGAAGAAGACAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATc                                                                                                                         >  1:1532954/1‑150 (MQ=255)
                               tATGGAAGACGGCAAGCAGAGTACTCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCg                                                                                                                        <  1:203341/150‑1 (MQ=255)
                               tATGGAAGAAGTCAAGCAGTGTAATCGTCTGGTTATTAAGACGCGCTAACACGTAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACTGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCg                                                                                                                        <  1:309837/150‑1 (MQ=255)
                                aTGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCgg                                                                                                                       >  1:1518989/1‑150 (MQ=255)
                                 tGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGACCAGGTGCAAACTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGATGTCCGTTTTAGCAGCggg                                                                                                                      >  2:454402/1‑150 (MQ=255)
                                 tGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCggg                                                                                                                      >  2:1163507/1‑150 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGTGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                     >  2:499015/1‑150 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                     >  1:824165/1‑150 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                     >  2:508144/1‑150 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                     <  1:517281/150‑1 (MQ=255)
                                   gaagaaGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGCCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAACTGGTACTGTTGTCCGTTTTAGCATCGGGCa                                                                                                                    >  1:647539/1‑150 (MQ=255)
                                     agaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGAAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAgg                                                                                                                  <  1:1555514/150‑1 (MQ=255)
                                           cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                            >  2:1195816/1‑150 (MQ=255)
                                           cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                            >  2:1706819/1‑150 (MQ=255)
                                           cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                            >  2:1401626/1‑150 (MQ=255)
                                            aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                           <  1:1633390/150‑1 (MQ=255)
                                            aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                           <  1:911446/150‑1 (MQ=255)
                                              tcAGAGTAATCGTCTGGTTATTAAGCCGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGTTTACTCTCCTGTTTATCCATACAGCAACAGTACTTGGGTAAACTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                         <  1:174145/149‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                         <  1:62792/150‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                         <  1:1498038/150‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                         <  2:1079739/150‑1 (MQ=255)
                                                agagTAATCGTCTGGTTATTAAGACGCTCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAAACCATAAATCAACAGTACTGGGGTAACCTGGTAGTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGACTGTa                                                                                                       >  1:5268/1‑150 (MQ=255)
                                                     aaTCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTg                                                                                                  >  2:1139946/1‑150 (MQ=255)
                                                         gTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTAAGCTTCGGGCAGAAAAAGCATGTACCCTGgaca                                                                                              >  1:257947/1‑148 (MQ=255)
                                                           cTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGc                                                                                            <  1:1121497/150‑1 (MQ=255)
                                                              gTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCtt                                                                                         <  1:73944/150‑1 (MQ=255)
                                                                tattaAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTc                                                                                       >  1:487879/1‑150 (MQ=255)
                                                                 attaAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCt                                                                                      <  2:1225218/150‑1 (MQ=255)
                                                                   taAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTg                                                                                    >  2:1661332/1‑150 (MQ=255)
                                                                          cgcTAACACGGAACATGTGCAAATTGGGCGTATTTTGATTACACTCCTGTTAATCGTTACAGCAACTGTTCTGGGGTAAGCTGGTACTGTTGTTCGTTTTAGCATCGGGCAGGAAAAGACTGTTACCTGGAGAGCCTTTCTTGAGCCAgg                                                                             <  1:1630023/150‑1 (MQ=255)
                                                                             tAACACGGAACAGGGGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGCGTAACCGGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCATGAAt                                                                          <  2:819022/150‑1 (MQ=255)
                                                                              aaCACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATa                                                                         >  2:1621628/1‑150 (MQ=255)
                                                                                cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                       >  2:376046/1‑150 (MQ=255)
                                                                                cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                       >  2:16755/1‑150 (MQ=255)
                                                                                cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                       >  1:489217/1‑150 (MQ=255)
                                                                                 acGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACAt                                                                      >  1:732249/1‑150 (MQ=255)
                                                                                  cGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACAtt                                                                     >  2:1324520/1‑150 (MQ=255)
                                                                                    gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                   <  2:671285/150‑1 (MQ=255)
                                                                                    gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                   >  2:1139047/1‑150 (MQ=255)
                                                                                    gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                   <  1:1399129/150‑1 (MQ=255)
                                                                                     aaCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGaa                                                                  >  2:1457410/1‑150 (MQ=255)
                                                                                      aCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAg                                                                 <  1:685073/150‑1 (MQ=255)
                                                                                      aCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAg                                                                 <  1:560145/150‑1 (MQ=255)
                                                                                         ggTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACa                                                              >  1:719473/1‑150 (MQ=255)
                                                                                          gTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAACAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACaa                                                             >  2:1271496/1‑150 (MQ=255)
                                                                                            gCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaag                                                           <  2:1469871/150‑1 (MQ=255)
                                                                                              aaaaTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaagaa                                                         <  1:1332134/150‑1 (MQ=255)
                                                                                               aaaTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaagaag                                                        <  1:1423224/150‑1 (MQ=255)
                                                                                                  tCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTca                                                     <  2:737948/150‑1 (MQ=255)
                                                                                                     gCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACAt                                                  <  2:545464/150‑1 (MQ=255)
                                                                                                       ggaTTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTg                                                <  2:1097686/148‑1 (MQ=255)
                                                                                                       gTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTg                                                <  1:1432411/150‑1 (MQ=255)
                                                                                                         aTTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGACTTTCTTTAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGAc                                              >  1:283997/1‑150 (MQ=255)
                                                                                                           tttGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTAATGTTGTCCGTTTTAGCATAGGCCAAGAAAAGCCTGTAACCTGGAGAGCCTGTCTAGAGCCAGTAATACATTGAAGACAAAGAAGTCACATTGACaa                                            >  1:227386/1‑150 (MQ=255)
                                                                                                                  acacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                     <  1:1325389/150‑1 (MQ=255)
                                                                                                                     cTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                  >  1:485271/1‑150 (MQ=255)
                                                                                                                     cTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                  >  1:1008790/1‑150 (MQ=255)
                                                                                                                         tGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                              <  1:810742/150‑1 (MQ=255)
                                                                                                                              aTCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccg                         >  1:591176/1‑150 (MQ=255)
                                                                                                                                 cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                      <  2:541624/150‑1 (MQ=255)
                                                                                                                                 cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                      >  1:1641971/1‑150 (MQ=255)
                                                                                                                                   tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                    <  1:508144/150‑1 (MQ=255)
                                                                                                                                   tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                    <  2:719473/150‑1 (MQ=255)
                                                                                                                                     ggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                  >  2:1647486/4‑150 (MQ=255)
                                                                                                                                     cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                  >  2:1458173/1‑150 (MQ=255)
                                                                                                                                      gacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg                 >  1:852167/3‑150 (MQ=255)
                                                                                                                                        cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg               >  2:629598/1‑150 (MQ=255)
                                                                                                                                        cAACAGTACTGGGGTAAACTGGTACTGTTGTCCGTTATAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTATCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGAAGCGGCCGCCGCCTTCTGGAg               >  1:1537710/1‑150 (MQ=255)
                                                                                                                                         aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAgg              <  2:1532954/150‑1 (MQ=255)
                                                                                                                                          aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGc             >  1:104160/1‑150 (MQ=255)
                                                                                                                                          aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCg            <  1:503053/150‑1 (MQ=255)
                                                                                                                                           cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCg            >  2:655946/1‑150 (MQ=255)
                                                                                                                                                cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCt       >  1:195385/1‑150 (MQ=255)
                                                                                                                                                 tGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg      <  1:190626/150‑1 (MQ=255)
                                                                                                                                                    ggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATc   <  1:1163507/150‑1 (MQ=255)
                                                                                                                                                    ggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAAAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATc   >  1:1337834/1‑150 (MQ=255)
                                                                                                                                                     gTAACGTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGTAGAGCCTTTCCTGAGCCAGGATTACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATcc  <  1:1118643/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCC  >  NC_000913/933010‑933308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: