Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,060 G→C 6.7% R236G (CGC→GGC)  astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0600GC 6.7% 316.0 / 8.9 100R236G (CGC→GGC) astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/48);  new base C (5/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGC  >  NC_000913/1829911‑1830200
                                                                                                                                                     |                                                                                                                                            
gCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACcgcgcg                                                                                                                                              >  2:1034248/1‑150 (MQ=255)
gCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACcgcgcg                                                                                                                                              >  2:1107996/1‑150 (MQ=255)
 cccTGCGCAGGCTGCCCTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGAGGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  2:190668/1‑150 (MQ=255)
 cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             <  2:1468680/150‑1 (MQ=255)
 cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  2:697623/1‑150 (MQ=255)
 cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  1:1312308/1‑150 (MQ=255)
     gcgcAGGCTGCCCTTCTGCCACGTCCACCCGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCCCGCGTCGGCCCACCGGCACAGATGGCGCTATAGTTACGGGAGCGAAAACCGTCTGTCTCCAGCACCGCGCGGGCAg                                                                                                                                         >  2:7286/1‑150 (MQ=255)
       gcAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCTCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgc                                                                                                                                       >  2:1600511/1‑150 (MQ=255)
        cAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTCCGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgcg                                                                                                                                      <  2:849637/150‑1 (MQ=255)
         aGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgcgc                                                                                                                                     <  1:1429812/150‑1 (MQ=255)
          ggCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGcgcg                                                                                                                                    <  2:1825355/150‑1 (MQ=255)
             tGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGtt                                                                                                                                 <  2:744686/150‑1 (MQ=255)
               cccGTCTGCCACTTCCACCAGCCGACTTTTACGCAGGGCGCGCACGCGGTCGATGTCCCACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCACCACCGCGCCTGCCCGCGCGGtgct                                                                                                                               <  1:462459/150‑4 (MQ=255)
                ccTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCGTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTGTGc                                                                                                                              >  1:599810/1‑150 (MQ=255)
                 cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                             <  2:1830803/150‑1 (MQ=255)
                  ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                            <  2:236184/150‑1 (MQ=255)
                    cagCCACTTCCCCCAGCCGACTTTTACGGAGGGCGCGCACGCGGGCGAAGGCACACTCCAGCGGCGGCCCACCGTCAAAGATGTCGATAAAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                          <  1:1502934/148‑1 (MQ=255)
                    cTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACCCTCACGCGTCGGCCCCCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                          >  2:1397789/1‑150 (MQ=255)
                    cTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGGCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                          <  1:641232/150‑1 (MQ=255)
                    cTGCCACTTCCACCAGCCGACTTTTACGGAAGGCGCGCACGCGGGCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                          <  1:680233/150‑1 (MQ=255)
                         aCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                     >  1:660562/1‑150 (MQ=255)
                         aCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCATCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGGTTGCGGATGTAcc                                                                                                                     >  2:1299867/1‑150 (MQ=255)
                          cTTCCACCCGCCGACTTTTTCGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGGCGGCCCGCCGTCATAGATGTCGATATAGGTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCt                                                                                                                    >  2:68271/1‑150 (MQ=255)
                                atggCCGATTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGa                                                                                                              >  1:680482/4‑150 (MQ=255)
                                  aGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                                                                                                            <  2:1585261/150‑1 (MQ=255)
                                  aGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                                                                                                            <  2:746610/150‑1 (MQ=255)
                                   gCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGAGATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGCCCGATgc                                                                                                           >  1:932365/1‑149 (MQ=255)
                                     cGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACg                                                                                                         <  1:1715590/150‑1 (MQ=255)
                                     cGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACg                                                                                                         >  2:650618/1‑150 (MQ=255)
                                            tACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTggg                                                                                                  >  2:658854/1‑150 (MQ=255)
                                                     cgcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGTCTTCCTggg                                                                                         >  1:1521988/1‑150 (MQ=255)
                                                     cgcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTggg                                                                                         <  1:885124/150‑1 (MQ=255)
                                                       cgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTAGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAt                                                                                       >  2:497038/1‑150 (MQ=255)
                                                        gcACGCGGGCGATGGCACACTCAAGCGTCGGCCCACCGTCAAAGATGCCGAAATAGATACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                                                                                      <  2:1785135/150‑1 (MQ=255)
                                                         cccgcgGTCGACGTCACACTCTAGCGTCGGCCCACCGGCAAAGATGTCGATAAAGTTACGGTTGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGGGCGGTTTGCGGATGTACCTGACCGCTGACGTCCTGGGCTTCCTGGGATaa                                                                                     <  2:700067/148‑1 (MQ=255)
                                                           cgcgGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaaaa                                                                                   <  2:1312308/150‑1 (MQ=255)
                                                                  gATGTCTCACTCACGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGTCTTCCTGGGATAAAAAGTGGGt                                                                            <  1:296041/150‑1 (MQ=255)
                                                                       cacaCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa                                                                       >  1:966271/1‑150 (MQ=255)
                                                                            gccaGCGTCGGCCCACCGTCAAAGATGCCGATATCGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGAACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGa                                                                  <  2:1375722/147‑1 (MQ=255)
                                                                                gCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGtt                                                              <  1:1053101/150‑1 (MQ=255)
                                                                                  gTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTc                                                            <  1:497038/150‑1 (MQ=255)
                                                                                  gTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTc                                                            <  1:377733/150‑1 (MQ=255)
                                                                                    cGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCgg                                                          <  2:1815959/150‑1 (MQ=255)
                                                                                      gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGTATAAAAAGTGGGTATAGATCGGATGTTTCGGCa                                                        <  1:157855/150‑1 (MQ=255)
                                                                                      gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATATAAAGTGGGTATAGATCGGAGGTTTCGGcc                                                        >  1:1797106/1‑149 (MQ=255)
                                                                                        ccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATc                                                      >  2:686450/1‑150 (MQ=255)
                                                                                        ccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATc                                                      >  2:175685/1‑150 (MQ=255)
                                                                                          aCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAg                                                    <  1:1281769/150‑1 (MQ=255)
                                                                                           ccGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGt                                                   >  2:1350927/1‑150 (MQ=255)
                                                                                             gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                 <  2:1557288/150‑1 (MQ=255)
                                                                                             gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                 >  2:883841/1‑150 (MQ=255)
                                                                                                aaaGATGTCGATATAGTTACGGTAGCGACAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              >  1:1663834/1‑150 (MQ=255)
                                                                                                aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCTCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  1:934225/150‑1 (MQ=255)
                                                                                                aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  2:1669260/150‑1 (MQ=255)
                                                                                                aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  2:669287/150‑1 (MQ=255)
                                                                                                aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  1:1302304/150‑1 (MQ=255)
                                                                                                 aaGATGGCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCa                                             <  2:1435791/150‑1 (MQ=255)
                                                                                                  aGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCTCGGTTTTCTGATGTACCTGACCGATGACGGCCTGTGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                            <  1:225850/150‑1 (MQ=255)
                                                                                                  aGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                            <  2:1273093/150‑1 (MQ=255)
                                                                                                     tGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaa                                         >  1:1663200/1‑150 (MQ=255)
                                                                                                     tGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaa                                         <  2:776931/150‑1 (MQ=255)
                                                                                                        cGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAACTg                                      >  2:891224/1‑150 (MQ=255)
                                                                                                         gATATAGTTACGGTAGCGAAAACCTACTTTCACCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGACCTGGGCTTCCTGGGATAAAAAGGGGGTATAGCTCGGATGTATCGGCAACATTTCTGCAATAAATGc                                     >  2:1015825/1‑150 (MQ=255)
                                                                                                          atatAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCCCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGcc                                    <  2:1049375/150‑1 (MQ=255)
                                                                                                             taGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGAAGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCttt                                 <  2:655730/150‑1 (MQ=255)
                                                                                                               gTTACGGTTGCGCAAACCTTCTTTCTCCAGCAACGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGTCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGCTGGTTCGGCATCAGTTCAGCAATAAATGCCTTTTg                               <  2:1207472/150‑1 (MQ=255)
                                                                                                                ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGc                              <  1:190668/150‑1 (MQ=255)
                                                                                                                    ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                          >  2:924033/1‑150 (MQ=255)
                                                                                                                    ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                          <  2:680482/150‑1 (MQ=255)
                                                                                                                    ggTAGCGAAAACCATCTTTCTCCAGCACCGCGCTGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTGCCTGGGATAAAAATTGGGTATAGATCTTATGTTTCGGCATCATTTCTGCAATAAATGCCTTTTGCCCgg                          <  2:226831/150‑1 (MQ=255)
                                                                                                                         cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGTTAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                     >  2:22479/1‑150 (MQ=255)
                                                                                                                         cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                     >  1:695195/1‑150 (MQ=255)
                                                                                                                          gcaaaaCTTCTTTCTCCAGCACCGCGCGGGCACGCGCGTTTTTCGGATGTACCTGACCGCTGACGGCCTGGGATTCCTGGGATAAAAAGTGGGTATAGATCGGATGTGTCGGCATCAGGTCTGCAAAAAATGCCTGTTGCCCAGTgaagc                    >  2:154009/3‑146 (MQ=255)
                                                                                                                           aaaaCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCa                   <  1:1034248/150‑1 (MQ=255)
                                                                                                                               ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga               >  1:1702067/1‑150 (MQ=255)
                                                                                                                               ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga               >  1:1426124/1‑150 (MQ=255)
                                                                                                                               ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga               <  2:523216/150‑1 (MQ=255)
                                                                                                                                 ttctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaa             <  1:699/150‑1 (MQ=255)
                                                                                                                                  tctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGTGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaaa            <  1:1512324/150‑1 (MQ=255)
                                                                                                                                   ctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAAt           >  1:105591/1‑150 (MQ=255)
                                                                                                                                   ctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAAt           >  2:48046/1‑150 (MQ=255)
                                                                                                                                   ctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTGTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAAt           >  2:53496/1‑150 (MQ=255)
                                                                                                                                    tttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATc          >  1:741948/1‑150 (MQ=255)
                                                                                                                                    tttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATc          <  1:650618/150‑1 (MQ=255)
                                                                                                                                    tttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATc          >  1:841497/1‑150 (MQ=255)
                                                                                                                                     ttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCg         >  1:24867/1‑150 (MQ=255)
                                                                                                                                       ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgc       >  2:1154074/1‑150 (MQ=255)
                                                                                                                                        tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGTTGTACCTGACCGATGCGGGCCTGGGCTTCCTGGGAGACAAAGTGGGTATAGAGCGGATGTTTCGGCAGCAGTTCTGCTATAAGTGCCTGTTGCCCGGTGCCGCAGAGAAATgcggcg      >  2:127211/1‑150 (MQ=255)
                                                                                                                                        tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTTTAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGGAGAGAAAATCGgcg      >  1:733260/1‑150 (MQ=255)
                                                                                                                                        ccgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATAGGATGTTTCGGCATCAGATCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg      >  1:914915/7‑150 (MQ=255)
                                                                                                                                         cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc     >  2:1830468/6‑150 (MQ=255)
                                                                                                                                         cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc     >  1:1645072/6‑150 (MQ=255)
                                                                                                                                         ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGGATAGATCGGATGTTTCGTCATCAGTTCTGCAATAAACGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc     >  2:360575/1‑150 (MQ=255)
                                                                                                                                          cAGGACCGCGAGGGAAGGCGCGGTTTGTGGATGTTCCTGACGGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg    <  1:99080/150‑1 (MQ=255)
                                                                                                                                          cAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCGTCCTGGGATAATAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg    >  1:1828276/1‑150 (MQ=255)
                                                                                                                                           caaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAGAAAAAGTGGGTATAGATCGGTTGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCgg   >  1:339576/4‑150 (MQ=255)
                                                                                                                                           aGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAGAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCgg   >  1:778945/1‑150 (MQ=255)
                                                                                                                                            gCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCggg  >  2:556959/1‑149 (MQ=255)
                                                                                                                                            gCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc  >  2:564284/1‑150 (MQ=255)
                                                                                                                                            aaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc  >  2:1400328/3‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGC  >  NC_000913/1829911‑1830200

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