Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,061 G→T 6.5% A235A (GCC→GCA astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0610GT 6.5% 284.0 / 9.2 100A235A (GCC→GCAastAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/49);  new base T (5/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 104 total aligned reads displayed.

CCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGC  >  NC_000913/1829912‑1830200
                                                                                                                                                     |                                                                                                                                           
cccTGCGCAGGCTGCCCTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGAGGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  2:190668/1‑150 (MQ=255)
cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  2:697623/1‑150 (MQ=255)
cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             <  2:1468680/150‑1 (MQ=255)
cccTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                                                                                             >  1:1312308/1‑150 (MQ=255)
  cTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGc                                                                                                                                           <  1:697623/150‑1 (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)
              cccTTCTGCCACTTCCACCAGCCGACTATTACGGATGGCGCGCACGCGGTCGATGTCCCACTCACGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgct                                                                                                                               <  1:1715774/150‑4 (MQ=255)
              cccGTCTGCCACTTCCACCAGCCGACTTTTACGCAGGGCGCGCACGCGGTCGATGTCCCACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCACCACCGCGCCTGCCCGCGCGGtgct                                                                                                                               <  1:462459/150‑4 (MQ=255)
               ccTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCGTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTGTGc                                                                                                                              >  1:599810/1‑150 (MQ=255)
                cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATTTAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                             <  2:654015/150‑1 (MQ=255)
                cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                             <  2:1830803/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                                                                                                         >  2:650618/1‑150 (MQ=255)
                                    cGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACg                                                                                                         <  1:1715590/150‑1 (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)
                                                                                   cGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCgg                                                          <  2:1815959/150‑1 (MQ=255)
                                                                                     gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGTATAAAAAGTGGGTATAGATCGGATGTTTCGGCa                                                        <  1:157855/150‑1 (MQ=255)
                                                                                     gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATATAAAGTGGGTATAGATCGGAGGTTTCGGcc                                                        >  1:1797106/1‑149 (MQ=255)
                                                                                       ccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATc                                                      >  2:175685/1‑150 (MQ=255)
                                                                                       ccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATc                                                      >  2:686450/1‑150 (MQ=255)
                                                                                         aCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAg                                                    <  1:1281769/150‑1 (MQ=255)
                                                                                          ccGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGt                                                   >  2:1350927/1‑150 (MQ=255)
                                                                                            gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                 >  2:883841/1‑150 (MQ=255)
                                                                                            gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                 <  2:1557288/150‑1 (MQ=255)
                                                                                               aaaGATGTCGATATAGTTACGGTAGCGACAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              >  1:1663834/1‑150 (MQ=255)
                                                                                               aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCTCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  1:934225/150‑1 (MQ=255)
                                                                                               aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  2:1669260/150‑1 (MQ=255)
                                                                                               aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  1:1302304/150‑1 (MQ=255)
                                                                                               aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                              <  2:669287/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)
                                                                                                     gTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaaa                                        <  2:1172923/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)
                                                                                                               ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGc                              <  1:190668/150‑1 (MQ=255)
                                                                                                                   ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                          <  2:680482/150‑1 (MQ=255)
                                                                                                                   ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                          >  2:924033/1‑150 (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:1426124/1‑150 (MQ=255)
                                                                                                                              ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga               <  2:523216/150‑1 (MQ=255)
                                                                                                                              ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga               >  1:1702067/1‑150 (MQ=255)
                                                                                                                                ttctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaa             <  2:1776314/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)
                                                                                                                                                     |                                                                                                                                           
CCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGC  >  NC_000913/1829912‑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: