Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 6.7% intergenic (‑38/+35) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9510AC 6.7% 470.9 / 17.3 159intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (72/76);  new base C (3/8);  total (75/84)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 160 total aligned reads displayed.

CGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGG  >  NC_000913/4470803‑4471100
                                                                                                                                                    |                                                                                                                                                     
cGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGCGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAt                                                                                                                                                      <  2:654265/150‑1 (MQ=255)
cGGATTTACCGGGATCTGACCGGAGCTGATGATCATATTGCCCACATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGATAATGATATCTCCTTTATTACAGCGGTAAAAAGTCTGGATCAt                                                                                                                                                      >  2:1361177/1‑150 (MQ=255)
  gATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATta                                                                                                                                                    >  2:47595/1‑150 (MQ=255)
  gATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTATCCTTTATTACAGCGGTAAAAAGTCTGGCTCATta                                                                                                                                                    >  2:1256622/1‑150 (MQ=255)
      tACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAgg                                                                                                                                                >  1:16636/1‑150 (MQ=255)
      tACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGGAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAGTATAgg                                                                                                                                                >  1:278860/1‑150 (MQ=255)
       aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                                                                               >  1:912580/1‑150 (MQ=255)
       aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                                                                               >  1:699632/1‑150 (MQ=255)
        ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                              >  1:660015/1‑150 (MQ=255)
            gATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGAACGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGACa                                                                                                                                          >  1:112965/1‑150 (MQ=255)
             aTCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAg                                                                                                                                         <  2:1728015/150‑1 (MQ=255)
               cTGCCCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                                                                       <  1:484118/150‑1 (MQ=255)
               cTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                                                                       <  1:1440956/150‑1 (MQ=255)
               cTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                                                                       >  1:818045/1‑150 (MQ=255)
                 gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGTATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                                                                     >  2:1703837/1‑150 (MQ=255)
                 gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCAGTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTGTATTACAGCGGTAAAAAGTCTGGCGCATTATAGGGAGCCAGCCtt                                                                                                                                     >  1:1374890/1‑150 (MQ=255)
                    cGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                  <  2:2059351/150‑1 (MQ=255)
                     ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                                                                 <  2:1044108/150‑1 (MQ=255)
                      gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatt                                                                                                                                >  1:1854412/1‑150 (MQ=255)
                      gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatt                                                                                                                                <  2:870038/150‑1 (MQ=255)
                           gatgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                                                                           >  2:256522/1‑150 (MQ=255)
                             tgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                                                                         >  1:761244/1‑150 (MQ=255)
                               atCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                                                                       >  2:1135093/1‑150 (MQ=255)
                                gcATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                                                                      <  1:1309223/149‑1 (MQ=255)
                                   tatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                                                                   >  1:1699163/1‑150 (MQ=255)
                                    atTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                                                                  <  2:1742285/150‑1 (MQ=255)
                                    atTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                                                                  <  2:1384855/150‑1 (MQ=255)
                                      tGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                >  2:46663/1‑150 (MQ=255)
                                      tGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                >  2:1170187/1‑150 (MQ=255)
                                       gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                               >  2:262477/1‑150 (MQ=255)
                                       gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                               >  2:1886996/1‑150 (MQ=255)
                                       gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                               >  2:1085373/1‑150 (MQ=255)
                                       gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                               >  2:1539939/1‑150 (MQ=255)
                                         ccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                             <  1:1782539/150‑1 (MQ=255)
                                           aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGAAATTAATTGGcc                                                                                                           >  2:1239634/1‑150 (MQ=255)
                                            gATCAACGCCCTTTACGTAAGGACCGATAGCTGCCGGTGCATTTTACGTCGCTATAGTTTTGCTCATGATTTCTCCTTTCTTCCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                                                                          <  2:109316/150‑1 (MQ=255)
                                            gATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                                                                          <  2:384172/150‑1 (MQ=255)
                                             aTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                         <  2:1237395/150‑1 (MQ=255)
                                             aTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                         <  2:548774/150‑1 (MQ=255)
                                               cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                                                                       <  1:1135093/150‑1 (MQ=255)
                                                aaCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                                                                      <  1:1971508/150‑1 (MQ=255)
                                                 aCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacc                                                                                                     <  1:760494/150‑1 (MQ=255)
                                                  cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcgcca                                                                                                    >  2:1985561/1‑150 (MQ=255)
                                                         tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                             <  1:1861602/150‑1 (MQ=255)
                                                           cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                                                                           >  1:1821938/1‑150 (MQ=255)
                                                            gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                                                                          <  1:1342279/150‑1 (MQ=255)
                                                              aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                        >  2:2110857/1‑150 (MQ=255)
                                                               aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                                                                       <  1:39697/150‑1 (MQ=255)
                                                                  aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                                    >  1:1071836/1‑150 (MQ=255)
                                                                   ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                                                                   >  1:607694/1‑150 (MQ=255)
                                                                   ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                                                                   <  2:2170204/150‑1 (MQ=255)
                                                                   ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                                                                   >  1:396554/1‑150 (MQ=255)
                                                                     gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                                                                 >  1:566133/1‑150 (MQ=255)
                                                                      aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                <  2:670315/150‑1 (MQ=255)
                                                                         gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                             >  2:97465/1‑150 (MQ=255)
                                                                          cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGCGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                                                            <  1:20267/150‑1 (MQ=255)
                                                                           tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                                                           <  1:713310/150‑1 (MQ=255)
                                                                              cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGACAGCACCACATTATGGGAAAACTCTTTTTAACAGTa                                                                        >  2:1937931/1‑150 (MQ=255)
                                                                                gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                                                                      >  2:2094113/1‑150 (MQ=255)
                                                                                 tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                                                                     >  1:1741855/1‑150 (MQ=255)
                                                                                  gCATTTTCCGTCTCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTAATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATAATGGGAAAACTCTTTTTCACAGTATTTg                                                                    <  2:1399349/150‑1 (MQ=255)
                                                                                  gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                                                                    <  2:809889/150‑1 (MQ=255)
                                                                                   cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                                                                   >  1:749239/1‑150 (MQ=255)
                                                                                         ccGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                                                             >  1:1244116/1‑150 (MQ=255)
                                                                                         ccGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTATGGGAAAACTCTTTTTCAAAGTATTTGCATTTGa                                                             >  2:16830/1‑150 (MQ=255)
                                                                                          cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                            <  2:449622/150‑1 (MQ=255)
                                                                                           gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc                                                           <  1:94276/150‑1 (MQ=255)
                                                                                              gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                        <  1:512929/150‑1 (MQ=255)
                                                                                               cgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                       >  1:645543/1‑150 (MQ=255)
                                                                                                    gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCa                                                  >  2:2136956/1‑150 (MQ=255)
                                                                                                    gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCa                                                  <  2:621999/150‑1 (MQ=255)
                                                                                                     ttttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGTCCAGCTCCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                 >  1:110717/1‑150 (MQ=255)
                                                                                                     ttttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTAACAGTATTTGCATTTGAGCGCGATATCAt                                                 >  1:116182/1‑150 (MQ=255)
                                                                                                      tttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                <  2:655119/150‑1 (MQ=255)
                                                                                                          cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                            >  2:797087/1‑150 (MQ=255)
                                                                                                            cATGATTTCTACTTTATTACAGCGGTAAAAAGTCTGGATCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGTATTTGAGCGCGATATCATTGGcgcg                                          <  2:352890/150‑1 (MQ=255)
                                                                                                                     tcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCAc                                 <  1:47595/150‑1 (MQ=255)
                                                                                                                     tcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTAACAGTATTTGCAGTTGAGCGCGATATCATTGGCGCGTGTTCGCAc                                 >  1:1009384/1‑150 (MQ=255)
                                                                                                                       cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAACGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTAGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                               >  2:1245015/1‑150 (MQ=255)
                                                                                                                        tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                              <  2:1356193/150‑1 (MQ=255)
                                                                                                                         ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCa                             <  2:1760039/150‑1 (MQ=255)
                                                                                                                             taCAGCGGGAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCTCGATATCATTGGCGCGTTTTCGCACGGCAAAGc                         <  1:1335247/150‑1 (MQ=255)
                                                                                                                              aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCt                        >  2:558323/1‑150 (MQ=255)
                                                                                                                                aGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                      <  2:1071836/150‑1 (MQ=255)
                                                                                                                                  cGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGAt                    <  2:699632/150‑1 (MQ=255)
                                                                                                                                     ataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCCACCGCAATTAATTGGCCAGCACCACATTATGGGACACCTCTTTTTCACAGTATTTGCATTTGATCGCGATATCCTTCGCGCGTTTTCGCACGGCAAAGCTGGTTGaa                 >  1:100150/3‑150 (MQ=37)
                                                                                                                                       aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTAGCACGGCAAAGATGGATGAAAc               >  1:1452452/1‑150 (MQ=255)
                                                                                                                                          aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg            <  1:1170187/150‑1 (MQ=255)
                                                                                                                                           gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt           >  2:1619318/1‑150 (MQ=255)
                                                                                                                                            tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGtt          <  1:2140375/150‑1 (MQ=255)
                                                                                                                                             cTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc         >  1:1051692/1‑150 (MQ=255)
                                                                                                                                              ttgCTCATTATAGGGAGCAAGTCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg        <  1:1361177/148‑1 (MQ=255)
                                                                                                                                              tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg        <  1:314254/150‑1 (MQ=255)
                                                                                                                                              tGGCTCATAATACGGAGCCAGACTTTTATTACCAACCTCAATTAAATGGCCAGCACCACATTATTGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg        <  1:1368233/150‑1 (MQ=255)
                                                                                                                                               ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAATATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg       <  1:1482572/150‑1 (MQ=255)
                                                                                                                                                 cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa     >  1:623148/1‑150 (MQ=255)
                                                                                                                                                  tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt    <  2:749239/150‑1 (MQ=255)
                                                                                                                                                   cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg   >  1:1302799/1‑150 (MQ=255)
                                                                                                                                                   cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg   >  2:1603267/1‑150 (MQ=255)
                                                                                                                                                    aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg  <  2:890696/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGG  >  NC_000913/4470803‑4471100

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