Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,153 A→T 5.8% K497* (AAA→TAA) ‡ fadH → 2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,1530AT 5.8% 399.7 / 9.3 134K497* (AAA→TAA) ‡fadH2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked
Reads supporting (aligned to +/- strand):  ref base A (67/58);  new base T (4/4);  total (72/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 142 total aligned reads displayed.

CACCGTGACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCTTAAGCCCG  >  NC_000913/3233005‑3233302
                                                                                                                                                    |                                                                                                                                                     
cacCGTGACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGcaacaa                                                                                                                                                      >  1:1394359/1‑150 (MQ=255)
 cccGGGACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACaaa                                                                                                                                                     <  2:738204/149‑1 (MQ=255)
    gTGACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGtt                                                                                                                                                  >  1:1298669/1‑150 (MQ=255)
     cgACGGCGGATCAGCTACAGGCTCTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTg                                                                                                                                                 <  1:1843214/149‑1 (MQ=255)
      gACGGCGTATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGc                                                                                                                                                <  2:1199769/150‑1 (MQ=255)
      gACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGc                                                                                                                                                <  2:242308/150‑1 (MQ=255)
       aCGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGATGGCAACAAAGTTGcc                                                                                                                                               >  1:1002855/1‑150 (MQ=255)
        cggcggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCca                                                                                                                                              >  1:421019/1‑150 (MQ=255)
        cggcggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCca                                                                                                                                              <  2:852688/150‑1 (MQ=255)
        cggcggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCca                                                                                                                                              >  2:584985/1‑150 (MQ=255)
         ggcggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCcat                                                                                                                                             >  2:1299928/1‑150 (MQ=255)
         ggcggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCcat                                                                                                                                             <  2:673210/150‑1 (MQ=255)
           cggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACTTTGTTGCcaacc                                                                                                                                           <  2:117205/150‑4 (MQ=255)
           cggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACTTTGTTGCcaacc                                                                                                                                           <  2:1496266/150‑4 (MQ=255)
           cggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCcatca                                                                                                                                           >  1:601023/1‑150 (MQ=255)
            ggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCcatcat                                                                                                                                          >  1:905832/1‑150 (MQ=255)
            ggATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGATGGCAACTTTGTTGCcaaccg                                                                                                                                          <  1:224245/150‑5 (MQ=255)
             gATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCatcatc                                                                                                                                         <  2:1386978/150‑1 (MQ=255)
                  gTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTtg                                                                                                                                    >  1:1826587/1‑150 (MQ=255)
                    tACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGCCGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCTCCGGTTGGCAACAAAGTTGCCATCATCGGTtgtg                                                                                                                                  >  1:1281327/1‑150 (MQ=255)
                     aCAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCCTCGACGGGATCGATCATCCGAAGGTATTGAGTTAACTCGATGTACTGCTCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTgg                                                                                                                                 >  2:316806/1‑150 (MQ=255)
                     aCAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTgg                                                                                                                                 >  2:1102638/1‑150 (MQ=255)
                      cAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGc                                                                                                                                <  1:451143/150‑1 (MQ=255)
                      cAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGc                                                                                                                                <  1:1307968/150‑1 (MQ=255)
                        ggCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCgg                                                                                                                              >  2:39670/1‑150 (MQ=255)
                        ggCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCgg                                                                                                                              >  1:572923/1‑150 (MQ=255)
                        ggCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCgg                                                                                                                              >  2:46050/1‑150 (MQ=255)
                        ggCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCAATCGACGGGATCGATCATCCGAAGGTATTGAGGTATCTCGATGTACTGCGCGACGAAGCGCCGGTTGGCAACTAAGTTGCCATCATCGGTTGTGGCgg                                                                                                                              >  1:1091176/1‑150 (MQ=255)
                        ggCGTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCCCGCACGCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCgg                                                                                                                              <  1:168613/150‑1 (MQ=255)
                         gCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCggg                                                                                                                             <  2:38848/150‑1 (MQ=255)
                          cTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCgggt                                                                                                                            >  1:1941357/1‑149 (MQ=255)
                           tttCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGAt                                                                                                                           >  2:1040936/1‑150 (MQ=255)
                             tCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTg                                                                                                                         <  2:1085602/150‑1 (MQ=255)
                              cGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGCCGGGATCGATCATCCGAAGGTAGTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTgg                                                                                                                        >  1:433725/1‑150 (MQ=255)
                                aTGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGtt                                                                                                                      >  2:1187297/1‑150 (MQ=255)
                                 tGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGttt                                                                                                                     >  2:188221/1‑150 (MQ=255)
                                  gAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAATGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGtttt                                                                                                                    <  2:514138/150‑1 (MQ=255)
                                  gAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGTGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGgttt                                                                                                                    <  1:85470/150‑1 (MQ=255)
                                   aaaCGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTg                                                                                                                   <  1:483124/150‑1 (MQ=255)
                                   aaaCGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTg                                                                                                                   <  2:1616444/150‑1 (MQ=255)
                                   aaaCGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACAGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTg                                                                                                                   >  1:1178120/1‑150 (MQ=255)
                                    aaCGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGa                                                                                                                  >  1:1061153/1‑150 (MQ=255)
                                    aaCGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGACGGTTGGCAACAAAGTTGCCATCATCGGTTTTGGCGGGATTGGTTTTGa                                                                                                                  >  1:1995422/1‑150 (MQ=255)
                                       gATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATAc                                                                                                               >  2:823273/1‑150 (MQ=255)
                                          ccTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGc                                                                                                            >  2:694402/1‑150 (MQ=255)
                                          ccTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGc                                                                                                            >  2:989048/1‑150 (MQ=255)
                                           cTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGAGTGGTTTTGATACGGCg                                                                                                           >  2:580234/1‑150 (MQ=255)
                                              gCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATg                                                                                                        <  1:1325857/150‑1 (MQ=255)
                                               ccAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGt                                                                                                       <  1:766896/150‑1 (MQ=255)
                                                cAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTAGTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGGGGCGGGCTTGGTTTTGATACGGCGATGTa                                                                                                      >  1:247172/1‑150 (MQ=255)
                                                cAGGGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTa                                                                                                      <  1:35288/150‑1 (MQ=255)
                                                  gTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTAtt                                                                                                    <  1:1999158/150‑1 (MQ=255)
                                                  gTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTAtt                                                                                                    >  2:1122589/1‑150 (MQ=255)
                                                    gggATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTa                                                                                                  <  1:975081/150‑1 (MQ=255)
                                                     ggATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTaa                                                                                                 <  2:1947513/150‑1 (MQ=255)
                                                      gATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAg                                                                                                <  2:380247/150‑1 (MQ=255)
                                                        tCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGGTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTc                                                                                              <  1:1692134/150‑1 (MQ=255)
                                                                cgcACTCCGCCCATCGCCGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCAGCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGc                                                                                      >  1:163446/1‑150 (MQ=255)
                                                                    cTCCGCCCATCGACTGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAAt                                                                                  <  1:1863441/150‑1 (MQ=255)
                                                                     tCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATc                                                                                 <  2:2014899/150‑1 (MQ=255)
                                                                      ccGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATcc                                                                                >  1:1666159/1‑150 (MQ=255)
                                                                      ccGCCCATCCACGGGATAGAGCATACGAAGGTATTTAGTTATATCGATGTACTGCGCGCCAAAGCGCCGGGTTGCAACCAAGTTGCCATAATCGGTTTTTGCGGGATTGGTTTTGATACGGCGATTTATTTAAGTCAGCCGGGCGAAtac                                                                                >  2:1146025/1‑148 (MQ=255)
                                                                                 cGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGaa                                                                     <  1:1742840/150‑1 (MQ=255)
                                                                                   ggATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAata                                                                   <  1:1177256/150‑1 (MQ=255)
                                                                                        gatcATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGcc                                                              >  1:1623484/1‑150 (MQ=255)
                                                                                        gatcATCCGAAGGTATTCAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTCGTTTTGATACGGCGATGAATTTAAGTCAGCCGGTCGAAACCACCAGACAGAATATCGcc                                                              >  1:149626/1‑150 (MQ=255)
                                                                                         atcatcCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCg                                                             <  1:714720/150‑1 (MQ=255)
                                                                                           catcCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCggg                                                           >  2:1515389/1‑150 (MQ=255)
                                                                                           catcCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCggg                                                           >  2:1943109/1‑150 (MQ=255)
                                                                                           catcCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCggg                                                           >  2:1067359/1‑150 (MQ=255)
                                                                                              ccGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTc                                                        >  2:634175/1‑150 (MQ=255)
                                                                                               cGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATAATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCt                                                       >  1:1754555/1‑150 (MQ=255)
                                                                                                gAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTg                                                      <  1:700512/150‑1 (MQ=255)
                                                                                                         gAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTaatgaatg                                             >  1:1355820/1‑150 (MQ=255)
                                                                                                         gAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTaatgaatg                                             >  1:529025/1‑150 (MQ=255)
                                                                                                          aGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATgg                                            >  2:1196903/1‑150 (MQ=255)
                                                                                                            ttctctcGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAACGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATgggg                                          <  2:1855359/147‑1 (MQ=255)
                                                                                                             tATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATggggg                                         <  2:344269/150‑1 (MQ=255)
                                                                                                              aTCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCCCCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGa                                        >  2:80829/1‑150 (MQ=255)
                                                                                                                ctcGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATc                                      <  2:1259097/150‑1 (MQ=255)
                                                                                                                 tcGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCg                                     <  2:471355/150‑1 (MQ=255)
                                                                                                                      gTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGt                                >  2:13214/1‑150 (MQ=255)
                                                                                                                       tACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTa                               >  1:1536358/1‑150 (MQ=255)
                                                                                                                         gtGCGCGACATAGAGACGGGTGGCAACAAAG‑TGCCATCATCGGGTGTGGCGGGATGGGTGTTGATACGGCGTTGTATTTTAGTCAGTCGGGCGAATCCACCAGTCAGAATATCGCCGGGTTCTGTAATGTATGGGGGATCGACAGTAGcc                            <  1:269088/149‑1 (MQ=255)
                                                                                                                            cgcgACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTa                          <  1:846392/150‑1 (MQ=255)
                                                                                                                              cgACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATTAATGGGGGATCGACAGTAGCCTaca                        >  2:1619386/1‑150 (MQ=255)
                                                                                                                              cgACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGGTCTGTAATGAATGGGGGATCGACAGAAGCCTaca                        >  2:944420/1‑150 (MQ=255)
                                                                                                                                aCAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTacaac                      <  2:1485862/150‑1 (MQ=255)
                                                                                                                                  aaaGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAg                    >  1:148118/1‑150 (MQ=255)
                                                                                                                                    ggcgcCGGTTGGCAACAAAGTTGCCATCATCGGGTGGGGCGGGAGTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGACTCCACCAGCCAGAATATCGCCGGGGTCTGTAATGACTGGGGGATCGTCAGTAGCCGACAACAGGc                  >  1:1420931/2‑150 (MQ=255)
                                                                                                                                    aGCGCCGGTTGGCAACAAAGTTGCCATCATCGCTTGAGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCGCGAATATCGCCGGGTTCAGAAATGAATCGGGGATCGACAGTAGCCTACAACAGGc                  >  1:295280/1‑150 (MQ=255)
                                                                                                                                       gcCGGGTGGCAACAAAGTTGCCATCATCGGTGGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGTATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCtgg               <  2:1091176/150‑1 (MQ=255)
                                                                                                                                         cGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCtggtg             >  2:772049/1‑150 (MQ=255)
                                                                                                                                          ggTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCtggtgg            <  2:635566/150‑1 (MQ=255)
                                                                                                                                            ttGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCt          <  2:1370222/150‑1 (MQ=255)
                                                                                                                                            atGGCAACTTTGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCt          >  2:1781725/2‑150 (MQ=255)
                                                                                                                                              ggCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCTTa        >  2:1161864/1‑150 (MQ=255)
                                                                                                                                               gCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCTTaa       <  2:1014794/150‑1 (MQ=255)
                                                                                                                                                caacaaAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCTTAAg      >  2:1201006/1‑150 (MQ=255)
                                                                                                                                                caacaaAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGAAGCCTACAACAGGCTGGTGGCTTAAg      >  1:195393/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCGTGACGGCGGATCAGTTACAGGCTTTCGATGAAACGATCCTCGCCAGTGGGATCGTGCCGCGCACTCCGCCCATCGACGGGATCGATCATCCGAAGGTATTGAGTTATCTCGATGTACTGCGCGACAAAGCGCCGGTTGGCAACAAAGTTGCCATCATCGGTTGTGGCGGGATTGGTTTTGATACGGCGATGTATTTAAGTCAGCCGGGCGAATCCACCAGCCAGAATATCGCCGGGTTCTGTAATGAATGGGGGATCGACAGTAGCCTACAACAGGCTGGTGGCTTAAGCCCG  >  NC_000913/3233005‑3233302

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: ‑