Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,166 G→A 7.7% intergenic (+77/‑58) 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,1660GA 7.7% 314.3 / 10.1 104intergenic (+77/‑58)lrp/ftsKleucine‑responsive global transcriptional regulator/DNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base G (53/43);  new base A (2/6);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 107 total aligned reads displayed.

CAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATT  >  NC_000913/933017‑933313
                                                                                                                                                     |                                                                                                                                                   
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:10819/150‑1 (MQ=255)
           ggACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTc                                                                                                                                          >  1:541624/1‑150 (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:517281/150‑1 (MQ=255)
                           ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                          >  1:824165/1‑150 (MQ=255)
                           ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                          >  2:508144/1‑150 (MQ=255)
                            gaagaaGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGCCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAACTGGTACTGTTGTCCGTTTTAGCATCGGGCa                                                                                                                         >  1:647539/1‑150 (MQ=255)
                              agaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGAAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAgg                                                                                                                       <  1:1555514/150‑1 (MQ=255)
                                    cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                                 >  2:1706819/1‑150 (MQ=255)
                                    cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                                 >  2:1195816/1‑150 (MQ=255)
                                     aaGCAGAGTAATCGTCTGGTTATTAAGCCGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTACCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                                >  2:281585/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:1498038/150‑1 (MQ=255)
                                       gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                              <  1:62792/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)
                                                                         cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAACGCCTGTAACCTGGAGTGCGTTTCTTTCGCCAGGAATACa                                                                            >  2:125729/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                                                                        <  1:1399129/150‑1 (MQ=255)
                                                                             gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                        <  2:671285/150‑1 (MQ=255)
                                                                             gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                        >  2:1139047/1‑150 (MQ=255)
                                                                              aaCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGaa                                                                       >  2:1457410/1‑150 (MQ=255)
                                                                               aCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAg                                                                      <  1:685073/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)
                                                                                                gTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTg                                                     <  1:1432411/150‑1 (MQ=255)
                                                                                                  aTTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGACTTTCTTTAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGAc                                                   >  1:283997/1‑150 (MQ=255)
                                                                                                           acacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                          <  1:1325389/150‑1 (MQ=255)
                                                                                                              cTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                       >  1:1008790/1‑150 (MQ=255)
                                                                                                              cTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                       >  1:485271/1‑150 (MQ=255)
                                                                                                                  tGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                   <  1:810742/150‑1 (MQ=255)
                                                                                                                       aTCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccg                              >  1:591176/1‑150 (MQ=255)
                                                                                                                          cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                           >  1:1641971/1‑150 (MQ=255)
                                                                                                                          cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                           <  2:541624/150‑1 (MQ=255)
                                                                                                                            tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                         <  2:719473/150‑1 (MQ=255)
                                                                                                                            tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                         <  1:508144/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)
                                                                                                                                   aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAGAGAAGTCACATTGACAAAGTTAAGTAGCGGACGCCGCCTTCTGGAGGc                  >  1:1521577/1‑150 (MQ=255)
                                                                                                                                   aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGc                  >  1:104160/1‑150 (MQ=255)
                                                                                                                                   aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCg                 <  1:503053/150‑1 (MQ=255)
                                                                                                                                    cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCg                 >  2:655946/1‑150 (MQ=255)
                                                                                                                                      gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGATAAGTAGCGGCCGCCGCCTTCTGGAGGCGtt               >  1:681074/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)
                                                                                                                                                  ccTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAGGCGTTGCTGATCCTTAt   <  2:1184088/150‑1 (MQ=255)
                                                                                                                                                  ccTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt   >  1:1715532/1‑150 (MQ=255)
                                                                                                                                                  ccTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt   >  2:1414681/1‑150 (MQ=255)
                                                                                                                                                   cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGATGCGTTGCTGATCCTTAtt  >  2:1227988/1‑150 (MQ=255)
                                                                                                                                                   cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAtt  >  2:160360/1‑150 (MQ=255)
                                                                                                                                                   cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAtt  >  1:1125049/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATT  >  NC_000913/933017‑933313

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: