Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,159 G→T 5.8% intergenic (+70/‑65) lrp → / → ftsK DNA‑binding transcriptional dual regulator Lrp/cell division DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,1590GT 5.8% 336.4 / 5.3 99intergenic (+70/‑65)lrp/ftsKDNA‑binding transcriptional dual regulator Lrp/cell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base G (47/46);  new base T (3/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 104 total aligned reads displayed.

CTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATC  >  NC_000913/933010‑933307
                                                                                                                                                     |                                                                                                                                                    
cTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg                                                                                                                                                      >  2:1206336/1‑150 (MQ=255)
cTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg                                                                                                                                                      >  1:744851/1‑150 (MQ=255)
cTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg                                                                                                                                                      <  2:407724/150‑1 (MQ=255)
 tGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGt                                                                                                                                                     >  2:473694/1‑150 (MQ=255)
  ggCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTa                                                                                                                                                    >  1:459493/1‑150 (MQ=255)
  ggCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTa                                                                                                                                                    >  2:178261/1‑150 (MQ=255)
     gTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAcc                                                                                                                                                 >  1:1409345/1‑150 (MQ=255)
      tCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCt                                                                                                                                                >  1:1027711/1‑150 (MQ=255)
       cAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTg                                                                                                                                               <  2:1504016/150‑1 (MQ=255)
             cacacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACtg                                                                                                                                         <  1:608647/150‑1 (MQ=255)
              acacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACtgt                                                                                                                                        <  1:1598464/150‑1 (MQ=255)
              acacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATCTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACtgt                                                                                                                                        <  2:1627355/150‑1 (MQ=255)
                acGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGGGCTAACACGGATCAGGTTCAAAATTGGCGTATTTTGATTAGACTCCTGTTAATCCATACATCAACAGTACCAGGTTACCCCAGTACtgttg                                                                                                                                      <  2:157593/150‑1 (MQ=255)
                  ggACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgttgct                                                                                                                                    <  1:1526500/150‑3 (MQ=255)
                       tACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGttt                                                                                                                               >  2:1399771/1‑150 (MQ=255)
                        aCGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTTGTTATTAAGACGCGCTAACACGGAACAGTTTCAAAATAGGCGTATTTTGATTAGACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTTGTACTGTTGTCCGtttt                                                                                                                              <  2:829859/150‑1 (MQ=255)
                        aCGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGtttt                                                                                                                              >  2:1348475/1‑150 (MQ=255)
                        aCGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGtttt                                                                                                                              <  1:1529042/150‑1 (MQ=255)
                        aCGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGtttt                                                                                                                              >  1:1176406/1‑150 (MQ=255)
                          gttgttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACAGGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAg                                                                                                                            <  1:1237218/150‑1 (MQ=255)
                               tATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCg                                                                                                                       >  1:223824/1‑150 (MQ=255)
                                 tGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCggg                                                                                                                     >  2:1360209/1‑150 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                    <  2:974047/150‑1 (MQ=255)
                                  ggAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGc                                                                                                                    >  1:1517739/1‑150 (MQ=255)
                                  ggAAGAAGGCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTGCGTTTTAGCATCGGGc                                                                                                                    <  2:221388/150‑1 (MQ=255)
                                   gaagaaGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCa                                                                                                                   >  2:403407/1‑150 (MQ=255)
                                   gaagaaGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGCAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCa                                                                                                                   >  1:1611187/1‑150 (MQ=255)
                                    aagaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTTCAATATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAg                                                                                                                  <  2:305678/150‑1 (MQ=255)
                                     agaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTCTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAgg                                                                                                                 >  1:588387/1‑150 (MQ=255)
                                        agTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGaaa                                                                                                              >  2:341055/1‑150 (MQ=255)
                                            aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                          <  2:1068362/150‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACCGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                        <  2:1531176/150‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                        <  2:1536989/150‑1 (MQ=255)
                                              gCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTg                                                                                                        >  2:1467027/1‑150 (MQ=255)
                                               cAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGt                                                                                                       >  2:293740/1‑150 (MQ=255)
                                               cAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGt                                                                                                       <  2:313387/150‑1 (MQ=255)
                                                         gTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGgaga                                                                                             >  1:1315742/1‑150 (MQ=255)
                                                             ggTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCt                                                                                         >  2:644572/1‑150 (MQ=255)
                                                             ggTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCt                                                                                         <  2:1593808/150‑1 (MQ=255)
                                                             ggTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCt                                                                                         >  2:461959/1‑150 (MQ=255)
                                                              gTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCtt                                                                                        <  2:670783/150‑1 (MQ=255)
                                                              gTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCtt                                                                                        >  2:1337502/1‑150 (MQ=255)
                                                              gTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTAGAGAGCCtt                                                                                        <  1:614280/150‑1 (MQ=255)
                                                                tattaAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTc                                                                                      <  2:10321/150‑1 (MQ=255)
                                                                tattaAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTc                                                                                      <  2:170766/150‑1 (MQ=255)
                                                                   taAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTg                                                                                   >  1:107835/1‑150 (MQ=255)
                                                                   taAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTg                                                                                   <  1:473694/150‑1 (MQ=255)
                                                                   taAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTg                                                                                   <  1:1348475/150‑1 (MQ=255)
                                                                    aaGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGa                                                                                  >  1:103237/1‑150 (MQ=255)
                                                                      taCGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGc                                                                                <  1:425303/149‑1 (MQ=255)
                                                                      gACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGc                                                                                <  1:404581/150‑1 (MQ=255)
                                                                      gACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGc                                                                                >  2:27020/1‑150 (MQ=255)
                                                                        cgcgCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCa                                                                              >  1:1316783/1‑150 (MQ=255)
                                                                        cgcgCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCa                                                                              <  1:644572/150‑1 (MQ=255)
                                                                          cgcTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAgg                                                                            >  2:166637/1‑150 (MQ=255)
                                                                             tAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAAt                                                                         <  2:107835/150‑1 (MQ=255)
                                                                                cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                      >  1:610806/1‑150 (MQ=255)
                                                                                cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                      >  1:1371678/1‑150 (MQ=255)
                                                                                 acGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACAt                                                                     <  1:293740/150‑1 (MQ=255)
                                                                                  cGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACAtt                                                                    <  2:344269/150‑1 (MQ=255)
                                                                                   ggAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTg                                                                   >  2:1326771/1‑150 (MQ=255)
                                                                                    gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGa                                                                  >  1:1544404/1‑150 (MQ=255)
                                                                                         ggTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACa                                                             >  1:839972/1‑150 (MQ=255)
                                                                                         ggTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACa                                                             >  1:1327691/1‑150 (MQ=255)
                                                                                          gTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACaa                                                            <  2:628479/150‑1 (MQ=255)
                                                                                            gCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaag                                                          >  1:1226515/1‑150 (MQ=255)
                                                                                               aaaTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaagaag                                                       <  2:1160670/150‑1 (MQ=255)
                                                                                                    ggCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTcaca                                                  <  2:1567956/150‑1 (MQ=255)
                                                                                                    ggCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTcaca                                                  <  2:1409345/150‑1 (MQ=255)
                                                                                                     gCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACAt                                                 >  2:101984/1‑150 (MQ=255)
                                                                                                            ttGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACaaa                                          >  1:1180025/1‑150 (MQ=255)
                                                                                                             tGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAg                                         >  2:576581/1‑150 (MQ=255)
                                                                                                                ttACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                      <  2:1544404/150‑1 (MQ=255)
                                                                                                                  acacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                    <  2:103237/150‑1 (MQ=255)
                                                                                                                   cacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTca                                                    >  1:1380827/1‑133 (MQ=255)
                                                                                                                   cacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTca                                                    <  2:1380827/133‑1 (MQ=255)
                                                                                                                     cTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                 >  2:576659/1‑150 (MQ=255)
                                                                                                                      tCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGc                                <  2:223824/150‑1 (MQ=255)
                                                                                                                           ttAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccc                           >  1:192010/1‑149 (MQ=255)
                                                                                                                                ccATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                      >  1:1339037/1‑150 (MQ=255)
                                                                                                                                ccATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                      <  1:994828/150‑1 (MQ=255)
                                                                                                                                 cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                     >  2:711135/1‑150 (MQ=255)
                                                                                                                                    cggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                  >  2:43196/5‑150 (MQ=255)
                                                                                                                                    aCAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                  >  1:142827/1‑150 (MQ=255)
                                                                                                                                     ggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                   >  1:518185/4‑148 (MQ=255)
                                                                                                                                     ggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                   <  2:518185/145‑1 (MQ=255)
                                                                                                                                     cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                 <  1:403407/150‑1 (MQ=255)
                                                                                                                                     cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                 >  2:419708/1‑150 (MQ=255)
                                                                                                                                      aGCAACAGTACTGGCGTAACCTGGTACTGTTGTCCGTTTTAGCCTCGGGCAGGAAAAGCCTCTAACCTGGAGAGCCTTTCTTGAGCCAGGAATCCATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg                <  2:1076175/150‑1 (MQ=255)
                                                                                                                                        cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg              >  2:702215/1‑150 (MQ=255)
                                                                                                                                        cAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg              >  2:332841/1‑150 (MQ=255)
                                                                                                                                         aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAgg             <  2:1428448/150‑1 (MQ=255)
                                                                                                                                          aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGc            <  1:61360/150‑1 (MQ=255)
                                                                                                                                          aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGc            <  1:461959/150‑1 (MQ=255)
                                                                                                                                             gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCATGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGtt         <  2:1611187/150‑1 (MQ=255)
                                                                                                                                               aCTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc       <  2:1517739/150‑1 (MQ=255)
                                                                                                                                               aCTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc       <  2:204301/150‑1 (MQ=255)
                                                                                                                                                cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCt      <  2:610806/150‑1 (MQ=255)
                                                                                                                                                   gggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGAt   <  2:459493/150‑1 (MQ=255)
                                                                                                                                                    ggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATc  >  2:1035226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATC  >  NC_000913/933010‑933307

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