Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,042 Δ1 bp 5.9% coding (287/1323 nt) atoE → short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,325,0420A. 5.9% 255.9 / 9.3 85coding (287/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base A (47/33);  new base . (3/2);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCC  >  NC_000913/2324894‑2325184
                                                                                                                                                     |                                                                                                                                              
gTGAAAATGTGGGGTGACGGTTGCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAg                                                                                                                                               <  1:296144/150‑1 (MQ=255)
       tgtgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCa                                                                                                                                        <  1:467334/150‑1 (MQ=255)
         tgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTcttt                                                                                                                                      <  2:668596/150‑4 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                    <  2:1093637/150‑6 (MQ=255)
             gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                  >  2:1302464/1‑150 (MQ=255)
               gACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTc                                                                                                                                >  1:573668/1‑150 (MQ=255)
                  ggTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGTCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACt                                                                                                                             >  1:1505992/1‑150 (MQ=255)
                  ggTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACt                                                                                                                             >  1:517176/1‑150 (MQ=255)
                        tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                       <  2:1695509/150‑1 (MQ=255)
                        tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                       <  1:974398/150‑1 (MQ=255)
                         ggAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCAGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCg                                                                                                                      >  1:336056/1‑150 (MQ=255)
                             cTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTc                                                                                                                  <  1:996634/150‑1 (MQ=255)
                             cTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTc                                                                                                                  <  1:410159/150‑1 (MQ=255)
                                gCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGt                                                                                                               <  1:67319/150‑1 (MQ=255)
                                   ggCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGc                                                                                                            >  1:1402866/1‑150 (MQ=255)
                                    gCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCTCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCt                                                                                                           >  1:955143/1‑150 (MQ=255)
                                          ggTATGCTGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtc                                                                                                     <  1:1666567/150‑1 (MQ=255)
                                          ggTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCCGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtc                                                                                                     >  1:1045392/1‑150 (MQ=255)
                                            tATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcat                                                                                                   <  1:1028560/150‑1 (MQ=255)
                                             aTGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCACAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatc                                                                                                  >  2:981584/1‑150 (MQ=255)
                                             aTGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCAGCTCTGCTCCGGTGAAAAGTTTTCTGCGTACTTCCGCCTACGCCGCAAAGACGCCCGTACAGGGCCTAATGCTGGTCACTTTATTAGGTTAAGTCGCTTGTGtcatca                                                                                                 >  2:692492/1‑150 (MQ=255)
                                              tGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatca                                                                                                 <  2:817884/150‑1 (MQ=255)
                                                     ggCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTgggg                                                                                          <  2:1520616/150‑1 (MQ=255)
                                                      gcgcTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGa                                                                                         >  1:104302/1‑150 (MQ=255)
                                                          ttATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttg                                                                                     <  2:647389/150‑1 (MQ=255)
                                                           tATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCACCGCCGCAAAGACGCCCGGACAGGGCGTCCTGCTGGTCTCTTTCTTCGGTTCCGTCGCTTGTGTCATCAACTGGGGAtttgg                                                                                    >  1:246042/1‑150 (MQ=255)
                                                             tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCAGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGTtttggtt                                                                                  >  1:24544/1‑150 (MQ=255)
                                                              catcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTTCAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                 >  2:644540/1‑150 (MQ=255)
                                                              catcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                 >  2:561067/1‑150 (MQ=255)
                                                              aatcGTTACCGGTCAGGCCCTTGCCAGCTCTGCTCGGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                 <  1:1293011/149‑1 (MQ=255)
                                                                tcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtttgg                                                                               <  1:598298/150‑1 (MQ=255)
                                                                 cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATttggtttggt                                                                              <  2:1605719/150‑1 (MQ=255)
                                                                    tACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGt                                                                           >  2:1171912/1‑150 (MQ=255)
                                                                     aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                          >  2:1012991/1‑150 (MQ=255)
                                                                       cGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCgg                                                                        >  1:252896/1‑150 (MQ=255)
                                                                        ggTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgc                                                                       >  1:1503144/1‑150 (MQ=255)
                                                                        ggTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCGCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgc                                                                       <  1:1692136/150‑1 (MQ=255)
                                                                            aTGCCCTTGCCAGCTCTGCTCCGGTGAAAAGGTGGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACGAACTGCGCTTCAGTCGCTTGTGGCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCaa                                                                   <  1:1039345/150‑1 (MQ=255)
                                                                              gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATg                                                                 >  2:1579126/1‑150 (MQ=255)
                                                                                ccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTg                                                                            >  1:981393/1‑137 (MQ=255)
                                                                                ccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTg                                                                            <  2:981393/137‑1 (MQ=255)
                                                                                ccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGtt                                                               <  2:252896/150‑1 (MQ=255)
                                                                                 cTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGttt                                                              >  2:598577/1‑150 (MQ=255)
                                                                                  ttGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTg                                                             <  1:1302464/150‑1 (MQ=255)
                                                                                            tGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGt                                                   <  2:605480/150‑1 (MQ=255)
                                                                                              cTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCg                                                 >  2:768678/1‑150 (MQ=255)
                                                                                               tCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGc                                                >  1:800662/1‑150 (MQ=255)
                                                                                                ccGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGct                                               >  1:1048539/1‑149 (MQ=255)
                                                                                                ccGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGcc                                               >  2:349144/1‑150 (MQ=255)
                                                                                                ccGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGcc                                               <  1:1008501/150‑1 (MQ=255)
                                                                                                   gTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCgg                                            >  2:542517/1‑150 (MQ=255)
                                                                                                    tGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGc                                           >  1:1170893/1‑150 (MQ=255)
                                                                                                     gAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCg                                          >  2:780686/1‑150 (MQ=255)
                                                                                                     gAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCg                                          >  1:553465/1‑150 (MQ=255)
                                                                                                      aaaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGa                                         >  2:496261/1‑150 (MQ=255)
                                                                                                       aaaGTTTTCTGCTTACTGCCGCCTCCGCCGCAAAGAAGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                        <  2:1095662/150‑1 (MQ=255)
                                                                                                       aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                        >  1:376637/1‑150 (MQ=255)
                                                                                                       aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                        <  1:980607/150‑1 (MQ=255)
                                                                                                       aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                        <  1:349144/150‑1 (MQ=255)
                                                                                                         aGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTGCGGTTCAGTCGCTTGGGTCAGCAGCTGGGGATTTGGTTTGGTGGTCGGCGCAATGTATGCCCGTGAAGTCGCCCGGCgcggc                                      >  1:189710/1‑146 (MQ=255)
                                                                                                           tttGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTccc                                    <  1:1565246/150‑1 (MQ=255)
                                                                                                               ctgcGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGt                                >  1:658130/1‑150 (MQ=255)
                                                                                                               ctgcGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGt                                >  2:1729391/1‑150 (MQ=255)
                                                                                                                 gcGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTc                              <  2:1474335/150‑1 (MQ=255)
                                                                                                                    tACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGa                           <  2:1045392/150‑1 (MQ=255)
                                                                                                                     aCTGCCGCCTCCGCCGCAAAGACGCCCGTACAGCGCGTCATGCGGGGCACTTTCTTCGGTTCAGTCGCTTGTGTCTTCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGATTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGAt                          <  1:536789/150‑1 (MQ=255)
                                                                                                                      cTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGGTCTGAtt                         >  1:1681896/1‑150 (MQ=255)
                                                                                                                         ccgccTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGGCCCTTTCTGCGGTTCCGTCGCTGGTGGCATCACCTGGGGAGTGGGGGTGGTTGTCGGCGCAAGGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTGCTGAGTATc                      >  1:58900/1‑150 (MQ=255)
                                                                                                                          ggccTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCAGTTTCTTCGGTTCAGTCGCTTGTGTCATCATCTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATcc                     <  1:1460717/149‑1 (MQ=255)
                                                                                                                            ccTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGt                   >  2:794866/1‑150 (MQ=255)
                                                                                                                            ccTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGt                   >  1:11875/1‑150 (MQ=255)
                                                                                                                             cTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGTTTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGAGTATCCGtt                  >  1:611761/1‑150 (MQ=255)
                                                                                                                                 gccgcAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCGGATTATCCGTTGCTc              >  1:1274604/1‑150 (MQ=255)
                                                                                                                                    gcAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCAtt           >  1:1691821/1‑150 (MQ=255)
                                                                                                                                    gcAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCAtt           >  2:602007/1‑150 (MQ=255)
                                                                                                                                      catgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc         >  2:644263/4‑150 (MQ=255)
                                                                                                                                      catgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc         >  2:399784/4‑150 (MQ=255)
                                                                                                                                      aaaGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc         <  1:692492/150‑1 (MQ=255)
                                                                                                                                       atgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc        >  1:1471935/3‑150 (MQ=255)
                                                                                                                                       aaGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc        >  2:566877/1‑150 (MQ=255)
                                                                                                                                       aaGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTGCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc        <  1:378604/150‑1 (MQ=255)
                                                                                                                                         gACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTg      <  2:605914/150‑1 (MQ=255)
                                                                                                                                           cGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg    >  2:1638254/1‑150 (MQ=255)
                                                                                                                                             cccGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc  >  2:667365/1‑150 (MQ=255)
                                                                                                                                             cccGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc  >  1:1709573/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCC  >  NC_000913/2324894‑2325184

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: