Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 631,584 T→C 5.0% M564T (ATG→ACG)  cstA → carbon starvation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,5840TC 5.0% 242.8 / 2.6 79M564T (ATG→ACG) cstAcarbon starvation protein A
Reads supporting (aligned to +/- strand):  ref base T (41/34);  new base C (2/2);  total (43/36)
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.67e-01

ACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCC  >  NC_000913/631435‑631729
                                                                                                                                                     |                                                                                                                                                 
tctgcctgcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGGGGGGCGACTTCCTCCATCAGGGCGGGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACGAGATGCTGGCAGGGATGGCGCTGAt                                                                                                                                                   <  1:1085919/149‑1 (MQ=255)
      tgcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCtg                                                                                                                                             >  1:299693/1‑150 (MQ=255)
       gcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCtgt                                                                                                                                            >  1:593277/1‑150 (MQ=255)
         tAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGc                                                                                                                                          <  1:1594396/150‑1 (MQ=255)
              tgctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg                                                                                                                                     <  1:157400/150‑1 (MQ=255)
              tgctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg                                                                                                                                     <  2:1609830/150‑1 (MQ=255)
              tgctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg                                                                                                                                     <  1:1445645/150‑1 (MQ=255)
              tgctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg                                                                                                                                     <  1:474885/150‑1 (MQ=255)
              tgctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg                                                                                                                                     <  1:1174929/150‑1 (MQ=255)
                ctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtggtg                                                                                                                                   >  1:304377/1‑150 (MQ=255)
                    ccaccgcgcTGTCCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                                                                                                               <  1:707309/145‑1 (MQ=255)
                    cAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                                                                                                               >  2:1293799/1‑150 (MQ=255)
                    cAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                                                                                                               >  1:935934/1‑150 (MQ=255)
                    cAACAGCGCTGTGCGTACTGGCGTGGCGCTACTTCCTCCATCAGCGCGTGGTCGATCCGCTCGGCGCCAATAACACTCTGTGGCCGCTGATTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                                                                                                               >  1:891832/1‑150 (MQ=255)
                     aaCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGgtgt                                                                                                                                  >  2:311589/1‑146 (MQ=255)
                     aaCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGgtgt                                                                                                                                  <  1:311589/146‑1 (MQ=255)
                     aaCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTgttgtt                                                                                                                              <  1:355284/150‑1 (MQ=255)
                               tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                                                                                                                    >  1:1081758/1‑150 (MQ=255)
                                gCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAAc                                                                                                                   >  2:34839/1‑150 (MQ=255)
                                     cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCGGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCaa                                                                                                              >  1:82557/1‑150 (MQ=255)
                                       ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACg                                                                                                            >  2:1222598/1‑150 (MQ=255)
                                        gCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGt                                                                                                           <  2:1368319/150‑1 (MQ=255)
                                               gCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCt                                                                                                    >  1:1369590/1‑150 (MQ=255)
                                                cTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTg                                                                                                   <  2:454800/150‑1 (MQ=255)
                                                  aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTggg                                                                                                 >  2:679443/1‑150 (MQ=255)
                                                  aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTggg                                                                                                 <  1:1516875/150‑1 (MQ=255)
                                                    ttCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTg                                                                                               >  1:1279132/1‑150 (MQ=255)
                                                      cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTggg                                                                                                 <  1:1398740/146‑1 (MQ=255)
                                                      cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTggg                                                                                                 >  2:1398740/1‑146 (MQ=255)
                                                      cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGgc                                                                                             >  2:1292865/1‑150 (MQ=255)
                                                        tccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGgcgc                                                                                           <  1:234363/150‑1 (MQ=255)
                                                          cATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTg                                                                                         >  2:68912/1‑150 (MQ=255)
                                                            tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGt                                                                                       <  2:253350/150‑1 (MQ=255)
                                                            tCAGGGCGGGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCGGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGt                                                                                       <  1:1238341/150‑1 (MQ=255)
                                                             cAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTa                                                                                      >  1:99001/1‑150 (MQ=255)
                                                             cAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTa                                                                                      >  2:223316/1‑150 (MQ=255)
                                                               gggCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTAcc                                                                                    <  1:245505/150‑1 (MQ=255)
                                                                       tCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCt                                                                            >  1:1294697/1‑150 (MQ=255)
                                                                         gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTgg                                                                          >  1:627081/1‑150 (MQ=255)
                                                                          aTCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGgc                                                                         >  1:434442/1‑150 (MQ=255)
                                                                          aTCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGACGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGgc                                                                         >  1:953375/1‑150 (MQ=255)
                                                                             cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGgctgc                                                                      >  2:385190/1‑150 (MQ=255)
                                                                                 gggCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGAt                                                                  >  2:1196485/1‑150 (MQ=255)
                                                                                 gggCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATCGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGCTGGCGCTGGTACCAACGGCCTGGCTGCTGAt                                                                  >  1:1104835/1‑150 (MQ=255)
                                                                                  ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGAtt                                                                 <  1:223316/150‑1 (MQ=255)
                                                                                    cggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTg                                                               <  1:34839/150‑1 (MQ=255)
                                                                                      gcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTa                                                             >  2:1496937/1‑150 (MQ=255)
                                                                                       cATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTAc                                                            <  2:833734/150‑1 (MQ=255)
                                                                                         ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTAccc                                                          <  2:475004/150‑1 (MQ=255)
                                                                                             cacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGAc                                                      >  2:782097/1‑150 (MQ=255)
                                                                                                       ccGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTgg                                            >  2:930860/1‑150 (MQ=255)
                                                                                                        cGCTGTTTGGCATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGc                                           <  1:362264/150‑1 (MQ=255)
                                                                                                         gCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCa                                          <  1:1196485/150‑1 (MQ=255)
                                                                                                            ggttGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGaa                                       >  2:708640/3‑150 (MQ=255)
                                                                                                                ggTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCg                                   >  2:1335050/1‑150 (MQ=255)
                                                                                                                 gTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGt                                  >  2:722695/1‑150 (MQ=255)
                                                                                                                  tATTGCCAACGAGATGCTGGCAGGGATGGCGCTGATGAGCTGTGCGGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGGCTGGGTGGCGCTGGTAGCAAGGGCCTGGCTGCTGATTTGTACCGTGACCGCAGGCTGGGAGAAAGCGtt                                 <  2:507588/150‑1 (MQ=255)
                                                                                                                  tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCa                                                  >  2:928709/1‑133 (MQ=255)
                                                                                                                  tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCa                                                  <  1:928709/133‑1 (MQ=255)
                                                                                                                   aTTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGttt                                <  1:1222598/150‑1 (MQ=255)
                                                                                                                        cAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGccc                           >  1:1326647/1‑150 (MQ=255)
                                                                                                                        cAACCAGATGCTGGCAGGGATGGCGCTGACGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGccc                           >  2:147793/1‑150 (MQ=255)
                                                                                                                          aCCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCgg                         >  1:957951/1‑150 (MQ=255)
                                                                                                                           ccAGATGCTGGCAGGGTTGGCGCTGATGCTCTGTGGCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGa                        <  2:891832/150‑1 (MQ=255)
                                                                                                                           ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGa                        >  2:1516688/1‑150 (MQ=255)
                                                                                                                           ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCAGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGa                        >  2:928087/1‑150 (MQ=255)
                                                                                                                             aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTATTGCTTTGTGGATATGGTACCAACGGCCTTGCTGCTGATTTTTACCCTGACCGCATGCTGTCAGAAAGAGTTTAGgccttatg                      <  2:362838/150‑6 (MQ=255)
                                                                                                                             aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATg                      >  2:869350/1‑150 (MQ=255)
                                                                                                                                tGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGa                   >  1:1089125/1‑150 (MQ=255)
                                                                                                                                 gCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGaa                  >  2:1031852/1‑150 (MQ=255)
                                                                                                                                   tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAg                <  1:722695/150‑1 (MQ=255)
                                                                                                                                    ggCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGt               <  1:679443/150‑1 (MQ=255)
                                                                                                                                     gCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTc              <  2:593277/150‑1 (MQ=255)
                                                                                                                                        gggATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGc           <  1:1292865/150‑1 (MQ=255)
                                                                                                                                         ggATGGCGCTGACGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCt          <  2:953375/150‑1 (MQ=255)
                                                                                                                                           aTGGCGCTGCTGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCAGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTc        >  2:1526866/1‑150 (MQ=255)
                                                                                                                                            tGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTcc       >  1:526872/1‑150 (MQ=255)
                                                                                                                                                gcTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGc   <  2:304377/150‑1 (MQ=255)
                                                                                                                                                gcTGACGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGc   <  1:147793/150‑1 (MQ=255)
                                                                                                                                                 cTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGcc  <  2:935934/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCC  >  NC_000913/631435‑631729

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: