Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,555,334 G→A A132V (GCC→GTC)  eutK ← putative ethanol utilization carboxysome structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,555,3340GA98.9% 347.3 / NA 90A132V (GCC→GTC) eutKputative ethanol utilization carboxysome structural protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (48/41);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGGCAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCACTG  >  NC_000913/2555186‑2555479
                                                                                                                                                    |                                                                                                                                                 
gTTAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACa                                                                                                                                                 <  1:1324289/151‑1 (MQ=255)
  tAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGc                                                                                                                                               >  2:475881/1‑151 (MQ=255)
  tAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGc                                                                                                                                               >  1:1024300/1‑151 (MQ=255)
  tAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGc                                                                                                                                               >  1:189538/1‑151 (MQ=255)
   aaTGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGcc                                                                                                                                              >  1:1456691/1‑151 (MQ=255)
            tgttgtAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCt                                                                                                                                     >  1:1307975/1‑151 (MQ=255)
             gttgtAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTg                                                                                                                                    >  2:1272735/1‑151 (MQ=255)
               tgtAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGcc                                                                                                                                  >  2:730962/1‑151 (MQ=255)
               tgtAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGc                                                                                                                                   <  1:936864/150‑1 (MQ=255)
                     acacCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCGCTTCTCCTGCCGTCAt                                                                                                                             <  2:1335700/150‑1 (MQ=255)
                      cacCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCAtt                                                                                                                            <  1:1056646/150‑1 (MQ=255)
                      cacCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTc                                                                                                                           >  2:670907/1‑151 (MQ=255)
                      cacCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTc                                                                                                                           <  1:982586/151‑1 (MQ=255)
                       acCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCAtt                                                                                                                            >  1:1244970/1‑149 (MQ=255)
                       acCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTcc                                                                                                                          >  2:1306591/1‑151 (MQ=255)
                         ccGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCt                                                                                                                        >  1:165113/1‑151 (MQ=255)
                         ccGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCt                                                                                                                        >  1:607695/1‑151 (MQ=255)
                            gCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGAc                                                                                                                     >  2:1200660/1‑151 (MQ=255)
                            gCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGAc                                                                                                                     <  2:189513/151‑1 (MQ=255)
                             cACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGGCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACg                                                                                                                    >  1:1001167/1‑151 (MQ=255)
                              aCCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGt                                                                                                                   >  1:85722/1‑151 (MQ=255)
                                  ccGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTg                                                                                                               >  2:170916/1‑151 (MQ=255)
                                  ccGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTg                                                                                                               >  2:959094/1‑151 (MQ=255)
                                   cGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTg                                                                                                               >  1:1127090/1‑150 (MQ=255)
                                     aCAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATg                                                                                                            >  2:333530/1‑151 (MQ=255)
                                      cAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGt                                                                                                           >  2:654040/1‑151 (MQ=255)
                                        ggTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTa                                                                                                         <  2:44971/151‑1 (MQ=255)
                                         gTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTaa                                                                                                        <  2:305765/151‑1 (MQ=255)
                                          ttAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGATCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGGGGGCAGCCACTTCTCCCGCCGTCATTCCCTGACGTACTGATGTTAAc                                                                                                       >  2:643099/1‑151 (MQ=255)
                                            aaTGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACa                                                                                                      <  1:381976/150‑1 (MQ=255)
                                            aaTGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAg                                                                                                     <  1:1009713/151‑1 (MQ=255)
                                                    ttGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCcaac                                                                                              >  1:166908/1‑150 (MQ=255)
                                                     tGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCcaaca                                                                                             <  1:1425426/150‑1 (MQ=255)
                                                     tGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCcaaca                                                                                             <  2:473128/150‑1 (MQ=255)
                                                      gagaCGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCcaacaa                                                                                            <  1:232981/150‑1 (MQ=255)
                                                       agaCGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACt                                                                                          >  1:1281075/1‑151 (MQ=255)
                                                        gaCGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTc                                                                                         <  1:1365368/151‑1 (MQ=255)
                                                         aCGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCg                                                                                        >  2:1153102/1‑151 (MQ=255)
                                                           gATAACGACTACTGCGTTTACGTAACGTCCCGTCAGAAAAGAGCTGTTCGTGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTCCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTc                                                                                      <  2:1471986/151‑1 (MQ=255)
                                                                cGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTcagcag                                                                                 <  1:269895/151‑1 (MQ=255)
                                                                  actactGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGAt                                                                               <  2:176883/151‑1 (MQ=255)
                                                                        gCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCcgc                                                                          <  2:293199/150‑1 (MQ=255)
                                                                         cGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGcga                                                                        <  1:825395/151‑1 (MQ=255)
                                                                          gTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGcgac                                                                       >  2:1047283/1‑151 (MQ=255)
                                                                              aCGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGAt                                                                    >  1:568306/1‑150 (MQ=255)
                                                                               cGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATa                                                                   >  2:301806/1‑150 (MQ=255)
                                                                                 tAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACt                                                                >  2:1032895/1‑151 (MQ=255)
                                                                                  aaCGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACt                                                                <  2:1242861/150‑1 (MQ=255)
                                                                                  aaCGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTg                                                               >  2:1189754/1‑151 (MQ=255)
                                                                                  aaCGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTg                                                               >  1:227187/1‑151 (MQ=255)
                                                                                    cGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGgc                                                             >  2:316416/1‑151 (MQ=255)
                                                                                       cccGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGgcgc                                                           <  2:622598/150‑1 (MQ=255)
                                                                                        ccGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGTCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGt                                                         <  2:145007/151‑1 (MQ=255)
                                                                                        ccGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGt                                                         <  1:1014295/151‑1 (MQ=255)
                                                                                        ccGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGt                                                         >  1:1173680/1‑151 (MQ=255)
                                                                                             aGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGt                                                    >  1:1087644/1‑151 (MQ=255)
                                                                                              gAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGtt                                                   <  1:1272839/151‑1 (MQ=255)
                                                                                              gAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGtt                                                   <  1:67013/151‑1 (MQ=255)
                                                                                                aaaGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTcc                                                 >  1:1394557/1‑151 (MQ=255)
                                                                                                aaaGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTcc                                                 >  1:120422/1‑151 (MQ=255)
                                                                                                   gagCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCtt                                              <  2:1337748/151‑1 (MQ=255)
                                                                                                     gCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTa                                             <  2:963922/150‑1 (MQ=255)
                                                                                                        gTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACggg                                         <  2:728050/151‑1 (MQ=255)
                                                                                                        gTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACggg                                         >  1:807991/1‑151 (MQ=255)
                                                                                                         ttCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGt                                        >  1:1395132/1‑151 (MQ=255)
                                                                                                           cGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTg                                      <  1:670980/151‑1 (MQ=255)
                                                                                                             aGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCt                                    <  1:1200757/151‑1 (MQ=255)
                                                                                                              gcgcATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCtt                                   <  1:706635/151‑1 (MQ=255)
                                                                                                               cgcATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTc                                  >  2:364894/1‑151 (MQ=255)
                                                                                                                 cATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCgg                                >  2:223697/1‑151 (MQ=255)
                                                                                                                    ttCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCt                              <  1:1352045/150‑1 (MQ=255)
                                                                                                                        ggCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTgg                            <  1:1153198/148‑1 (MQ=255)
                                                                                                                        ggCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCgg                         >  2:1129650/1‑151 (MQ=255)
                                                                                                                        ggCTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCg                          <  1:607374/150‑1 (MQ=255)
                                                                                                                          cTTTTTCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGtt                       >  1:714398/1‑151 (MQ=255)
                                                                                                                            ttttCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGtt                       <  1:857292/149‑1 (MQ=255)
                                                                                                                             tttCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTaa                     <  1:1306697/150‑1 (MQ=255)
                                                                                                                             tttCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTaa                     >  1:1251247/1‑150 (MQ=255)
                                                                                                                             tttCCAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTaa                     >  2:428510/1‑150 (MQ=255)
                                                                                                                                 cAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGcc                 >  2:273762/1‑150 (MQ=255)
                                                                                                                                 cAGCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCa                <  1:1354878/151‑1 (MQ=255)
                                                                                                                                   gCGGCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAg               <  2:436287/150‑1 (MQ=255)
                                                                                                                                     ggCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTa             >  1:1419365/1‑150 (MQ=255)
                                                                                                                                     ggCCAGCCAAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTa             <  1:1047379/150‑1 (MQ=255)
                                                                                                                                         agccaAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAg        <  2:85711/151‑1 (MQ=255)
                                                                                                                                            caAAGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGcc      >  1:76407/1‑150 (MQ=255)
                                                                                                                                             aaaGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCAc    >  2:788715/1‑151 (MQ=255)
                                                                                                                                              aaGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCAc    >  1:1076119/1‑150 (MQ=255)
                                                                                                                                              aaGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCACt   >  1:803893/1‑151 (MQ=255)
                                                                                                                                               aGTGGACAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCACTg  <  2:1307869/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTTAATGCCGGGTGTTGTAGGACACCCGGCACCTCCGACAGGTTAATGGGGCTTGAGACGATAACGACTACTGCGTTTACGTAACGTCCCGGCAGAAAAGAGCTGTTCGAGCGCATTTCTGGCTTTTTCCAGCGGCCAGCCAAAGTGGGCAGCCACTTCTCCTGCCGTCATTCCCTGACGTACTGATGTTAACAGCGCCAACAACTCGTCAGCAGATTCCGCGACGATAACTGGCGCGTCGGGTTCCCTTACGGGTTGCTTCGGCTGGCGGTTAAAGCCAGTAACCAGCCACTG  >  NC_000913/2555186‑2555479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: