Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,126,514 T→C 18.0% E226G (GAG→GGG)  wcaG ← bifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase/ GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,126,5140TC18.0% 111.8 / 7.4 61E226G (GAG→GGG) wcaGbifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase/ GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase
Reads supporting (aligned to +/- strand):  ref base T (32/16);  new base C (3/9);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTCCGTACTGGCGGTTGT  >  NC_000913/2126279‑2126742
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
ccccGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGCCGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCA‑TGCGCCCGCCCca                                                                                                                                                                                                                                    >  1:1067623/1‑239 (MQ=255)
ccccGCTTCCAGGGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCTGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCca                                                                                                                                                                                                                                    <  1:1254269/239‑1 (MQ=255)
           gTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGCCAACATCGGCTGGGGGGTCTCCAGCCCAACCTCC‑TGCGCCAGCCCcacgcaacgaaa                                                                                                                                                                                                                          >  2:703137/1‑228 (MQ=255)
           gTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGGGCCGACGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCA‑TTGGCCAGCTCcagtcaatgaagg                                                                                                                                                                                                                         >  1:678422/1‑228 (MQ=255)
                aTTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCCCGCGGATAGTGCCGTCAACCCCCCTGCCGCGCTTAATGTGCGACAACACCGTCTGGGGGTTCCCCAGCCA‑ACATCACTTGCCCAGCCCCCAGACAAGAATGCTcgc                                                                                                                                                                                                                    >  2:1293222/1‑239 (MQ=255)
                    cGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTcgccgc                                                                                                                                                                                                                 >  1:1290154/1‑238 (MQ=255)
                    cGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGGGCAGTCAACGCCCGTGCCGCCGTTAATGTGCGACACCATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTcgcccc                                                                                                                                                                                                                 >  1:405383/1‑236 (MQ=255)
                    cGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGATTGCG‑CCAGCTCGCGGCTAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGCCAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTcgccgc                                                                                                                                                                                                                 >  1:456432/1‑238 (MQ=255)
                        aTACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCAGTGCCGCCGTACATGGGAGACAACATCGCCTGGGTGTACTCCAGCAAGCCTTCA‑CGGGCCCGCTCCAGTACCTGAATGCtgccccacccca                                                                                                                                                                                                            >  2:701013/1‑228 (MQ=255)
                                       aTGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCa                                                                                                                                                                                             >  1:1234198/1‑239 (MQ=255)
                                        tGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATAGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGCCATGAATGCTCGCCGCTGCCATATCATCGACGTGCAg                                                                                                                                                                                            >  1:696405/1‑239 (MQ=255)
                                             gcgcGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGCCGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCCAATCATCGACGTGCAGAAAt                                                                                                                                                                                        >  1:1076109/1‑238 (MQ=255)
                                                              gtgtGCGCGGGCTGCCTCCCTGTTTGTTGCCCTCAAAACCCGCCCGGCC‑TTGTAACCCCACCCCTTGGGCGGTGGTTGGCG‑CCAGCTCGCGGTTAGTGCAGTCAACGCCCGTGCCGACTTTAATGTGCGAAAACATCGGTTGGGTGTCTCCCAGCCAGACTTCA‑TGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCg                                                                                                                                                                      <  2:167643/236‑1 (MQ=255)
                                                               tttGCGGGGCGTGCCATCCGGTTTGGTGGCAACAAAAACCCCCCGGCCTTTGTAC‑CCCACCCCTTTGGCGATGGTTTGCG‑CCAGCCCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCg                                                                                                                                                                      <  1:876109/238‑1 (MQ=255)
                                                                ttGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAACCCCCCCGGCCTTTGTAC‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCg                                                                                                                                                                      <  1:812183/237‑1 (MQ=255)
                                                                 tGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCGGACTTCA‑TGCGCCTGCTCCCTGACATGAATGATCGCCGTCGCCATATCATCGACGTGCAGCAATTCGCGCATCGGTGTACCGCt                                                                                                                                                                    >  1:601137/1‑238 (MQ=255)
                                                                            ctacccGTTTTGCGGGTCTTAAACACCCCCCGTCTTTTGTCA‑CCCCCCACTTTGGGGATGTTTTGGG‑CCTGTCCGGGGAGTGGGCGTGCAACGCCTGTCCCGACGTTAATGTGCGACAACCTCGGGGGGGTGTTCTCCCCGCAGACTTCA‑GGCCCAAGCTCCATGACGGGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccac                                                                                                                                                        <  2:423645/235‑1 (MQ=255)
                                                                             cATCCGGTTTGCTGGCATCAAAACCCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccac                                                                                                                                                        <  1:1293222/238‑1 (MQ=255)
                                                                                 cgttttGCTGGCATCAAAAACCCCCCGGCCTTTGTAA‑CCCTCCCCTTTGGCGATGGTTTGCG‑CCAGCTCGGGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAAAATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGt                                                                                                                                                   <  1:964704/236‑1 (MQ=255)
                                                                                 ccgttgggggtcATCAAAACCCCCCCGGCCTTTGTAC‑CCCACCCCTTGGGCGATGG‑TGGCGCCCAGCCCGCGGATCGGGGAGTGCGCGCCCGTGCGGACGTTAATGGGCGAAAACATCGGCTGGGTGTTCCCCAGCCGGACTTCA‑TGCGCGCGCTCAATGACATGAATGCTCGCCGCCGCCATCTCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATCCCACCACGt                                                                                                                                                   <  2:1227948/228‑1 (MQ=255)
                                                                                 cGGTTTGCTGTGATAAAACCCCCCCCGCCCTTTCTCC‑CCCACCACTTGGGCGGTGTTTTGCG‑CCACCGCGCGGAGTGGGCCGCAAACGCCCGGGCGGACGTTAGTGTGCGACAACTTGGGCTGGGTGTTCCCCCGCCAGTCTTCA‑TGCGCACGCTCCATGACATGAATGCTCGCCCCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGt                                                                                                                                                   <  2:1234198/239‑1 (MQ=255)
                                                                                    tttGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTcc                                                                                                                                                 >  1:158103/1‑238 (MQ=255)
                                                                                      tGCTGGCACAAAAAACCCCCCGGCCTTTTTAC‑CCCCCCCTTTGGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGgc                                                                                                                                              <  1:888837/239‑1 (MQ=255)
                                                                                      tGCCGGCACCAAAACCCCCCCGGCCTTTGTAA‑CCCACCCCTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGGGCAGTCACCGCCCGTGCCGACGTTAATGTGCGCCAACATCGGCTGGGGGTTCTCCAGCCAGACTTCA‑TGCGCCAGTTCCATGACATGAATGCTCGCGGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGgc                                                                                                                                              <  1:243365/239‑1 (MQ=255)
                                                                                         tGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TTCGCCAGCTCCATGACAGGAATGCTCGCCGCCCCCATATCCTCCACGTGCAGAAATTCGCCCATCGGGGGACCCCCGCCCCAAACCCCCACCCCCGgcacc                                                                                                                                           >  2:387381/1‑236 (MQ=255)
                                                                                                     ccaccCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGCGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACACCGTCCGGCCCATCCTGTGcctt                                                                                                                                >  1:1200074/1‑236 (MQ=255)
                                                                                                       accCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCCGCCAGACTTCA‑TGCGCCAGCTCCAGGACATGAATGCTCGCCTCCGCACAATCATCGACGTGCATAAATTCGCCCATCGGGGTCCCCCTGCCCCATACCAACACGCCCCGCGCATTCTGTGCCGTCg                                                                                                                              >  2:629549/1‑238 (MQ=255)
                                                                                                            gCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGt                                                                                                                        >  1:420329/1‑239 (MQ=255)
                                                                                                                     accccACCCTTTGGGCTATGGTTTGCG‑CCAGCTGGGGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGAAAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGa                                                                                                                <  1:996261/237‑1 (MQ=255)
                                                                                                                         caccacTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCCCGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgc                                                                                                             >  1:258958/1‑238 (MQ=255)
                                                                                                                            cacTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGGGCGACGACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGCGTCCCGCTGCCCCCTACCCCCCCGGCCCGCGCATTCTTTGCCCGCGCCCCCAGGAAGCGACggcaga                                                                                                         >  2:510635/1‑234 (MQ=255)
                                                                                                                              cTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATAC‑‑CCACGTCCGGCGCATTCTGTGCCGTCGCCTCGGGGAAGCGACGCAGAAATGc                                                                                                     >  1:507127/1‑239 (MQ=255)
                                                                                                                                           tggcg‑cCAGCTGGGGGAGCGTGCAGTCAACCCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGTGGTTCTCCAGCCAGACTTCA‑TGCGCCATCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATcaca                                                                                          <  1:1214905/236‑1 (MQ=255)
                                                                                                                                               g‑cCAGCCCGCGGATAGGGCAGTCAACGCCCGTGCCGACGTTAATGGCCAACAACATCGGCTGGGTGTTCCCCAGCCAGCCTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCg                                                                                      <  1:703137/238‑1 (MQ=255)
                                                                                                                                                          ggATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCt                                                                                                                           <  1:587527/192‑1 (MQ=255)
                                                                                                                                                              aGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCCGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGATCGCCCCCCCCCTATCCTCACCGTGCAGAAATTCGCGCTACGGTTTACCGCTGCCCCATACCACCACGTCCGCCGCATCTCGTCCCGTCGCCTCGTGGAGCCCACCCACCAACGC‑TGGGCTCACATGGGCAttcaccggggagc                                                                       >  2:1066818/1‑228 (MQ=255)
                                                                                                                                                                gcaaGTCAACCCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTTTTCCCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAg                                                                     <  1:387381/236‑1 (MQ=255)
                                                                                                                                                                    gTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGGGTACCGCTGCCCCATACCACCACGCCCGGCGCATTCTGTGCCGTCGCCCCGTGGAAGCg                                                                                                                 >  2:587527/1‑192 (MQ=255)
                                                                                                                                                                       aaCGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCTCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGGTTCCCGCCGCCCCATACCACCACGTCCGGCGCCTTCTTTGCCGTCGCCTCATGGCAGCGCCGCAGCAATGC‑TGGGATACCCTGCGCATTACTGgggggggaagtgtggt                                                              >  2:943282/1‑229 (MQ=255)
                                                                                                                                                                          gCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGCCGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCCCCTGGACGCGACGCACCAATGC‑TGGGATCACCTGCGAATTCCCCGGGGGGAAGtatggggggg                                                           >  2:1004662/1‑230 (MQ=255)
                                                                                                                                                                               tGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGCCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGGATTACTCGGGGGGAAGTTGTCGTGGGGCCCg                                                      >  1:1055158/1‑239 (MQ=255)
                                                                                                                                                                                   gccGTTAATGTGCGACAACATCGGCTGGGCGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCCTGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCTCGCATCGGTGTACCGCTGCCCCATACCACCACCTCCGGCGCATTCTGTCCCTGTGCCTCCGGGAAGCGA‑GCAGCAACACGTGGGACCACTTGCGAATGACTCGGGGGGCAATGGACGGGTGGCCCGTGCa                                                  >  2:1006313/3‑239 (MQ=255)
                                                                                                                                                                                        tAATGGGCGACAACCTCGGCTGGTTTTTCGCCAGCCAGACTTCA‑TGCCCCAGCCCCATTACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGGCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTg                                             <  1:128639/239‑1 (MQ=255)
                                                                                                                                                                                        tAATGGGCGACAACAGCGGCTGGGTGTTCCCCAGCCAGACTTCA‑GGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTg                                             <  1:568200/239‑1 (MQ=255)
                                                                                                                                                                                               cGACAACATCGGCTGGGTGTTTTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTTGCCGCCCCCATATCATCGACGCGCAGAAATTCGCGCATCGGTGTCCCGCTGCCCCATACCACCACGTCCGGCGCGTTCTGTGCCGTCTCCTCGCGGAACCCACTCAGACATGC‑TGGGCTCTCACCCACATTACTCGGATGGACGATTTGCTGTGGCCCGCTCAggcggggacgcg                                      >  2:232919/1‑229 (MQ=255)
                                                                                                                                                                                               cGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATa                                                                                                                                                            >  1:748200/1‑122 (MQ=255)
                                                                                                                                                                                                gACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGCATTACTCGGGTGGAAGGTGTCGTGTGGCCCGTACAGGTTGGTGGGCAt                                     >  1:634750/1‑239 (MQ=255)
                                                                                                                                                                                                gACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCAt                                     >  1:845926/1‑239 (MQ=255)
                                                                                                                                                                                                caaaaactgggCTGGGGTTTCCCCAGCAAGCCTTAA‑GGCGCCAGCCCCATGCACGGAATGCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCAt                                     <  1:1172947/230‑1 (MQ=255)
                                                                                                                                                                                                   actcgggctggggTCTTCCCCGCCCAAACTTCA‑TGCCCCAGTACCATGCCAGGACTGCTCGCCGCCGCCATATCCTCGACCTGGAGAAATTCGCGAGTTGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGAc                                  <  2:271361/229‑1 (MQ=255)
                                                                                                                                                                                                               gttttCCCCAGCCAGCCTTCA‑GGCCCCAGCTCCATGCCATGAATGCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAAt                       <  1:356047/235‑1 (MQ=255)
                                                                                                                                                                                                                   tccccagccagCCTTCA‑TGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGt                  <  1:1004662/236‑1 (MQ=255)
                                                                                                                                                                                                                   tccccagccagCCTTCA‑TGCCCCAGCCCCATGAAATGAATGCTCCCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGt                  <  1:629549/236‑1 (MQ=255)
                                                                                                                                                                                                                   tccccagccagACTTCA‑GGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGt                  <  1:510635/236‑1 (MQ=255)
                                                                                                                                                                                                                     tcCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCCTCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAAAGC‑TGGGCTCACCAGCGCATTACTCGGGGGGAAGTAGTCGGGAGGCCCGGACAGTGTGGGCCGCGTCACCGGCCGGAAATCGGGTcc                >  2:1162841/1‑239 (MQ=255)
                                                                                                                                                                                                                     tcCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTAACGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGACGCGACGCAGCACAGC‑Tggg                                                                                                >  1:807059/1‑159 (MQ=255)
                                                                                                                                                                                                                        agccagACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccacc                                                                                                                                                       >  1:399009/1‑102 (MQ=255)
                                                                                                                                                                                                                          ccagACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTCGTCGTGTGGCCCGTACAGGTTGGTCGGCATGCCTGACCGGTAATCGCGTCCgg              >  1:660013/1‑235 (MQ=255)
                                                                                                                                                                                                                          ccagACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGAAGCAATCC‑TGGGATCACATGCGAATTACTCGGGGGGCAGTTGTCGGTGGGCCCGTAACGGGTGGGCGGCATGACTGGGCGGTAAAcgccgccggaac           >  2:948204/1‑230 (MQ=255)
                                                                                                                                                                                                                          ccagACTTCA‑GGCGCCAGGCCCATAAACGGAACGCCCCCCCCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTCCGTACt           <  1:455451/239‑1 (MQ=255)
                                                                                                                                                                                                                                   a‑TGCGCCAGCCCCGTGACATTATTGCTCGCAGCCGCCATATCATCGACGTGCAGAAATTCGCGCAGCGGTGTACCGCTGCCCCATACCACCACGTCCGGCCCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTCCGTACTGGCGGTTGt  <  1:1162841/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
CCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAA‑CCCACCACTTTGGCGATGGTTTGCG‑CCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCA‑TGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGC‑TGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTCCGTACTGGCGGTTGT  >  NC_000913/2126279‑2126742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑