Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,126,490 T→C 15.2% E234G (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,4900TC15.2% 164.6 / 3.3 75E234G (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 (36/25);  new base C (3/10);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 3.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACTGGTAAGTGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTC  >  NC_000913/2126257‑2126727
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
aCCACTGGTAAGTGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGGGCCAGCCCGCGGGTAGGGCCGTCAAAGCCCGTGCCGCGGTTAATGTGCGCAAAAACCGGCGGGGGGGTCCccag                                                                                                                                                                                                                                             >  2:168252/1‑238 (MQ=255)
aCCACTGGTAAGTGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTccagc                                                                                                                                                                                                                                            <  1:858244/239‑1 (MQ=255)
          aGTGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCCCGGTAAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCACCAACACCGGCTGGGGGTGTCCCAGCCAAACCTCag                                                                                                                                                                                                                                  >  2:434625/1‑238 (MQ=255)
            tGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCTTCAAAAACCACCCGGCCTTTGTAACCCACCCCTTTGCCGATGGTTTGCGCCACCTCCCGGATAGTGCAGGCAAGCCCCGTCCCGCCCTTAATGTCGCACAACCTCGGCAGGGGTGTGCCCAGCCAAACTCAATgc                                                                                                                                                                                                                                >  2:43310/1‑239 (MQ=255)
            tGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATg                                                                                                                                                                                                                                 >  1:949776/1‑238 (MQ=255)
            tGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGCCACCGACCGTGCCGACGTTAATGTGCGACAACATCGGCTCGGGGTTCTCCAGCCAGACTTAATg                                                                                                                                                                                                                                 >  1:1124971/1‑238 (MQ=255)
            tGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGGGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGCTAGTGCAGTCAACGCCCGGGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGCATTCAtc                                                                                                                                                                                                                                 >  1:138519/1‑237 (MQ=255)
               tGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGTGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGAAAACATCGGCTGGGTGGTCTCCAGCCAGACTTCATgcgc                                                                                                                                                                                                                              >  1:291664/1‑238 (MQ=255)
               cggaaaGCCCCGCTTCCCGGGGGATTTCGGGCATCCCGCCAAGCTGGATCCGGCGGCGCCCATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCCTCAAAAACCCCCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGcc                                                                                                                                                                                                                             <  2:264530/235‑1 (MQ=255)
                  cAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAg                                                                                                                                                                                                                           <  1:1216194/238‑1 (MQ=255)
                    aGCCCCGGTTCCAGTGAGATTGCGTGATACCAGCCAAGCTTATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCAGTTTGCTGGCATCAAAAACCCCCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTc                                                                                                                                                                                                                        <  1:655756/239‑1 (MQ=255)
                    aGCCCCGCCTCCAGTGAGATTTCGGGATACCCGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTc                                                                                                                                                                                                                        <  1:811054/239‑1 (MQ=255)
                      ccccGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGGGCGACAAAACGGGCTGGGGGTCCTCCACCCCAACTTCATggccccccccca                                                                                                                                                                                                                      >  2:362875/1‑228 (MQ=255)
                           gtttccggtatgtttctttaacccAGCCAAGCTGGTGCAGGCTCTTCCCATCCAGCGATTTGCGCGGCGTGCCAGCCGGTTTGGTGGCCTCAAAAACCCCCCGGCCTTTGTAACCCACCACTTTGGCGAGGGTTTGGGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCcatgaca                                                                                                                                                                                                                 <  2:228117/217‑1 (MQ=255)
                            ttccgtggggttttcgggtccccAGCCAAGCTGATGCAGGCCCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCCCCCGGCCTTTGGAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCcatgacat                                                                                                                                                                                                                <  2:745545/218‑1 (MQ=255)
                                       tttctgtatcccACCCAAGCTGTGCCGGGCCCTTCACATCCAGAATTTTGCGGGGCGTGCCTCCGTGTTTGCGGGCATCAAAACCCCCCGGGCTTTTGAACCCCACCCCTTTGGCGAGGTTTTGCGCCAGCTCGCGGATGGTGCAGTCAACCCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTcgcc                                                                                                                                                                                                     <  2:1101227/229‑1 (MQ=255)
                                                             aTGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCCCCATCTCCATGCCATGAATGCTCGCCCCCGCCCTATCATCGCCggaca                                                                                                                                                                               >  1:369957/1‑235 (MQ=255)
                                                                       gaaccaaccgcaggttctggcggggCGCCACCGGGTTGGCGGGCACAAAAAACCCCCCGGCTTTTTAACCCCACCCCGTTGGCGAGGTTGTGCGCAAGCCCGCGGATAGTGCAGTCACCGCCCGGGCCGACGTAAATGTGCGACAACTGCGGTTGGGTGTTCTCCAGCCAGACTTAATGCGCCAGCTCCATGACATGAAGGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTcgcg                                                                                                                                                                     <  2:1098358/220‑1 (MQ=255)
                                                                       gTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGCCAACATCGGCTGGGAGTTCTCCAGCCAGCCTTCATGCGCCAGCTCCATGACATGAATGCTCGCCCCCGCCATATCATCGCCCTGCAGAAATTcgc                                                                                                                                                                      >  1:1190036/1‑238 (MQ=255)
                                                                        tacccccccggaggttGGGCCGCGGGCCCCCCCGTTTGCTGGCAACCAAAACCCCCCCGGCCTTTGTAACCCCCCACTTTGGCGAGGGTTGGCGCCAGCTCGCGGATAGGGCAGTCAACGCCCGGGCCGACGTTAATGGGCGACAACCACGGCTGGGTTTTCTCCAGCCAGACTTCATGCGCCAGCTCCATAACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTcgcg                                                                                                                                                                     <  2:1190036/225‑1 (MQ=255)
                                                                        tCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGGGGAGAAATTCgcgc                                                                                                                                                                    >  1:1087528/1‑239 (MQ=255)
                                                                        tCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGAGTTTTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCCGAAATTCgcgc                                                                                                                                                                    >  1:100491/1‑239 (MQ=255)
                                                                          ccctcccggcgttttgggggcgTCCCTTTCGTTGTGCGGGCAATAAAACCCCCCCGGCTTTTGAAACCCCCCGCTTTGGCGATGGTTGGCGCCAGCCCGCGGATGGGTCAGTCAACGCCCGTGCCGACGTTAATGTGCGAAAACATCGGCTGGGTGTTCCCCAGCCAGACTTCATGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCAt                                                                                                                                                                  <  2:421167/221‑1 (MQ=255)
                                                                                  cagTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGAGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCCTATCATCCCCGTGCAGAAATTCGCCCATCGGGGTa                                                                                                                                                           >  1:257677/1‑238 (MQ=255)
                                                                                      ttGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATACCCCAGCTCCATGACCTGAATGTTCGCCCCCGCCAAACCATCGACGTGCCGAAATTCCCGACTCggggaccccc                                                                                                                                                       >  2:79375/1‑230 (MQ=255)
                                                                                       tGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCcat                                                                                                                                                                                                                     >  2:64588/1‑175 (MQ=255)
                                                                                       tGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCcat                                                                                                                                                                                                                     <  1:64588/175‑1 (MQ=255)
                                                                                         cgcgGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGGGTACCGCTGc                                                                                                                                                    >  1:1159668/1‑238 (MQ=255)
                                                                                                    ccccccgtgtgtcgggaacaaaaaaCCCCCCGGCCTTTTACACCCACCCTTTTGGGGTTGGTTGGCGCCAGCCGGGGGTGAGTGCAGTCACCCCCCGTGCCGACGTTAATGTGCGAAAACATCGGCGGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccac                                                                                                                                         <  2:949776/219‑1 (MQ=255)
                                                                                                        cgttttgggggCATAAAAAACCCCCGGGCTTGTGAACCCCCCCGCTTTGGCGAGTGTTTGGGCCAGCCGGGGGATGGTGCAGTCAACCCCCGTGCCGACGTTAATGCGCGAAAACACCGGCTGGGTGTCCCCCAGCCAGACTTCATCCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGt                                                                                                                                    <  2:1105076/230‑1 (MQ=255)
                                                                                                             tGCTGGCATCAAAAACCACCCGGCATTTGTAACCCACCACTTTAGTGATGGTTTGCGCAAGCCCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGCGTTCCCCATCCCAACTTCATGCGCCAGCTCCATGAAATGAAAGCTCGCTGCCGCCATGTCAGCGACGTGCAGAAATTCACGCACCAGGGTACCGTCGCCCCACACCACCACGTCCgcg                                                                                                                               >  2:310975/1‑237 (MQ=255)
                                                                                                                      tCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCCTATCACCCCCCGTCAGAAATTCCCTCACCCGGGTACCGCGTCCCCAACCACCCACCTCCGGACCATTCtt                                                                                                                       >  2:363791/1‑238 (MQ=255)
                                                                                                                           acccaccCGGCCTTTGTAACCCCCCACTTTGGCGATGGTTTGCGCCAGCTCGGGGATAGTGCAGTCACCGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCCCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCg                                                                                                                  <  1:362875/237‑1 (MQ=255)
                                                                                                                               ccccGGCTTTTGCAACCCACCAGTTGGGGGGTGGTTTTCGCCTGCCCCGGTGAGGTACAGCAAACCCCGGTCGCGAGTTTAATTGTCGACACAAACGGCTGGGTGTTCTCCAGCCCGACTTCAGGCGCCAGTTCCATGACGTGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGc                                                                                                              <  2:940010/238‑1 (MQ=255)
                                                                                                                               accCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATCCTGTGCCGTcc                                                                                                               >  1:37542/1‑237 (MQ=255)
                                                                                                                                        ttGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCCGAAATTCGCCCATCGGTGTACCGCTGCCCCATACCCCCACGCCCGGCCCCTTCTGtcccctcccccctgggaa                                                                                                     >  2:14616/1‑231 (MQ=255)
                                                                                                                                               ccaccaCTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTCCCGTCGCCTCGTGGAAGCGACgc                                                                                             >  1:483394/1‑239 (MQ=255)
                                                                                                                                                 ccccgtggggggtgtgTTTGGCCCGGCTGGCGGGTGGGGCAGTCACCGCCGGCCCCGGCTTTATGGCGCAACAACCCCGGCTGGGTGTTCCCCAGCCAGACTTCATGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgcag                                                                                           <  2:37542/224‑1 (MQ=255)
                                                                                                                                                 cccccTTTGGCTATGGTTTGCGCCAGCTCGCGGATAGGGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgcag                                                                                           <  1:500184/235‑1 (MQ=255)
                                                                                                                                                 accacTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACCTGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCCCGCATCGGGGCACCCCTTCCCCATACCACCCCGTCCGGGGCATTCTCTGCCGACGCCCCGTGGAGGAGAccacg                                                                                           >  2:127094/1‑235 (MQ=255)
                                                                                                                                                 accacTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCCCATTCTGTCCCGTCGCCTCGTGGAAGCGACgcag                                                                                           >  1:519196/1‑239 (MQ=255)
                                                                                                                                                 accacTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCCCGCATCGGGGTACCGCTGCCCCATACCACCACGCCCGGCGCATTCGTTGCCGCCCCCCCGGGGGAGCgccccca                                                                                           >  2:906481/1‑233 (MQ=255)
                                                                                                                                                  cctcTTTGGCGATGGTTTGCGCCAGCGCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGTTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgcagc                                                                                          <  1:906481/236‑1 (MQ=255)
                                                                                                                                                   cacTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgcagc                                                                                          >  1:1233884/1‑238 (MQ=255)
                                                                                                                                                   cacTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGCGCGGAAATTCCCGCATCGGTGTACCGCTGCCCCCTACCACCCCGTCCGGCCCCTTCTTGCCCGTCGCCCCGGGGAAACgccccacca                                                                                         >  2:412985/1‑231 (MQ=255)
                                                                                                                                                      ttgggCGATGGTTTGCGCCAGCTCGCGGATAGGGCAGTCAACGCCCGTGCCGACGTTAATGGGCGACAACATCGGCTGGGTGTTCCCCAGCCAGACTTCATGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATg                                                                                      <  1:1075319/236‑1 (MQ=255)
                                                                                                                                                             aTGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCCgggg                                                                                >  1:485608/1‑237 (MQ=255)
                                                                                                                                                              ttgTTTGCGCCAGCTCGGGGATAGTCCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGAt                                                                               <  1:434625/236‑1 (MQ=255)
                                                                                                                                                                           cccgggggtAGTGCACTCACCCCCCGTGGCGACGTTAATGGGCGACAACAGGGGGTGGGGTCTCTCCAGACAGACTTCATGCGCCCGCTCCATGACATGAATGCGCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTAc                                                                 <  2:892110/231‑1 (MQ=255)
                                                                                                                                                                           cTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCACCGGGGTCCCCCTGCCCCCACCCACCACGTCCGCCGCATTCTTGGCCTCGCCCCACGGGAAGCGCCGCAGACATGCCGGGAACACAGGCGGATTa                                                                  >  2:822751/1‑238 (MQ=255)
                                                                                                                                                                            tCGGGGATAGTGCCGTCACCGCCGGTGCCGACGTTAATGTGCGCAAACATCGGCTGGGTGTTCTCAAGCCAGACTTCATGCGCCAGCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACt                                                                <  1:498657/239‑1 (MQ=255)
                                                                                                                                                                            ccgggggttgggcgatcagccccccGTGCCGCGTTGTATGGGCGAAAACTGCGGGGGGGTTTCCTCCCGCCAGACTTCCGGCGCCCTCTCCATGCCATGAATCCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTAc                                                                 <  2:387829/217‑1 (MQ=255)
                                                                                                                                                                                       tGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGc                                                                                                                                                       >  2:265165/1‑144 (MQ=255)
                                                                                                                                                                                       tGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGc                                                                                                                                                       <  1:265165/144‑1 (MQ=255)
                                                                                                                                                                                         cggTCACCGCCCGTGCGCAGGTTAAGGGGCAACAACCTCGGTTGGGTTTTCCCCAGCCCGGCTTCCGGCCCCAGCTCCATGAAATGAATGCCCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGtt                                                    <  1:704959/237‑1 (MQ=255)
                                                                                                                                                                                          ggccAACCCCCTGGCCGCCGTTAATGGGGGACAACATGGGTTGGGTTTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGtt                                                    <  1:363791/235‑1 (MQ=255)
                                                                                                                                                                                             cAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAGTTCGCGCATCGGTGTACCGCTGCCCCATACCCCCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGGAACGCCCCAGCCATGCCGGGAACCCATGCCAATTACTCGGGGGGGAGCTCTCg                                                >  2:1325784/1‑239 (MQ=255)
                                                                                                                                                                                                 gCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCCACGCAGCAATGCTGGGCTCACATGCGAATTACTCGGGTGGAAGGTGGCGtgtg                                            >  2:1032458/1‑239 (MQ=255)
                                                                                                                                                                                                 gCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGACCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACgcagca                                                                                         >  1:53350/1‑194 (MQ=255)
                                                                                                                                                                                                 gCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGGGTTCTCCAGCCAGACTTCCTGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCAAATCATCGACGTGCAGAAATTCCCGCATCGGTGTCCCGCTGCCCCCGACCCCCAAGTCCGCCCCTTCCTGTCCCTCTGCCTCGTGCAGGCGACGCACCAGACCCAGGACCACATGCGATTTACCCGGGGGGGAGTGGTCggtgg                                            >  2:881625/1‑235 (MQ=255)
                                                                                                                                                                                                   ccGTGCCGACTGTAATGGCGCACAACATGGGCTGGGTGTCCCCCAGCCAGACTTCATGCGCCAGCCCCATGAAATGAATGCTCCCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGc                                          <  1:412985/239‑1 (MQ=255)
                                                                                                                                                                                                    gggcccGACTTTAGTGGGCAACCACATCGGCTGGGTCTCTCCCACCCAGACTCCATGCACCATCTCCATGACATGAATGCCCGCCTCGCCCATATCATCGACGTGCAGAAATTGGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGcc                                         <  2:915721/235‑1 (MQ=255)
                                                                                                                                                                                                      tGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGCGGAAGCGACCCAGCACGCCTGGGACCACCTGCGGATTACCCGGGGGGAAGTTGTCGTGGGGCCCg                                       >  2:704959/1‑239 (MQ=255)
                                                                                                                                                                                                        ccGACGTTAAGGTGCGACAACATCGGGTGGGTGTTCTCCAGCCAGACTTCATGCGCGAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGt                                      <  1:79375/238‑1 (MQ=255)
                                                                                                                                                                                                            cGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCCCCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCCCATGCGAATTCCCCGGGGGGAGGTTGTCGTGTGGCCCGTTCaag                                 >  2:939542/1‑237 (MQ=255)
                                                                                                                                                                                                             gTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCCTACCACCACGTCCGGCGCATTCTGTGCCGCCGCCTCGTGGAAGCCACGCACCACTGCTGGGATCACATCGGAATTACTCGGGGGGGAGTTTTCGTGGGGCCGGCACGGGt                                >  2:1377817/1‑239 (MQ=255)
                                                                                                                                                                                                               tAATGTGCTACAACCTCGGCTGGGTGTCCTCCAGCCAGACTTCGTGCGCCAGCTCCATGGCATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTg                              <  1:573522/239‑1 (MQ=255)
                                                                                                                                                                                                                aaagggccGAAAACATGGGCTTGGTTGTCTCCAGCCAGCCTTCATGCGCAAGCTCCATGACCTGAAGGCTCGCGGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTg                              <  1:1046325/232‑1 (MQ=255)
                                                                                                                                                                                                                       cGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccac                                                                                                                                         >  2:776370/1‑126 (MQ=255)
                                                                                                                                                                                                                       cGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATaccac                                                                                                                                         <  1:776370/126‑1 (MQ=255)
                                                                                                                                                                                                                       cGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTCGGATCACATGCGAATTACTCGGGTGGAAGTTGCCGTGGGGCCCGTACCAGTTGGTCGGCa                       >  1:675815/1‑239 (MQ=255)
                                                                                                                                                                                                                        gACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTAc                                                                 >  1:773293/1‑196 (MQ=255)
                                                                                                                                                                                                                              caTCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCCCATGCGACTTACTCGGGTGGACGTTGTCGGGGGGCCCGTAAAGGTTgggggggcggcatg                 >  2:185496/1‑227 (MQ=255)
                                                                                                                                                                                                                                         gtgtTCCCCCGCCCGCCTTCATGGCCCATCCCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAAt        <  1:822751/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gtgtTCCCCAGCCAGCATTCATGCGCCGTCCCCAGGACATGAACCCCCCCCCCCGGCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATGCCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAAt        <  1:658048/238‑1 (MQ=255)
                                                                                                                                                                                                                                              ctcCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCACTGCCGGGATCACATGCGAATTACTCGGGGGGAAGTTGTCGTGGGGCCCGGCCAGGTTGGTGGGCCTggatgcgggggaaaggactc  >  2:656511/1‑220 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
ACCACTGGTAAGTGCTGGCAAGCCCCGCTTCCAGTGAGATTTCGTGATACCAGCCAAGCTGATGCAGGCGCGTCACATCCAGCAGTTTGCGCGGCGTGCCATCCGGTTTGCTGGCATCAAAAACCACCCGGCCTTTGTAACCCACCACTTTGGCGATGGTTTGCGCCAGCTCGCGGATAGTGCAGTCAACGCCCGTGCCGACGTTAATGTGCGACAACATCGGCTGGGTGTTCTCCAGCCAGACTTCATGCGCCAGCTCCATGACATGAATGCTCGCCGCCGCCATATCATCGACGTGCAGAAATTCGCGCATCGGTGTACCGCTGCCCCATACCACCACGTCCGGCGCATTCTGTGCCGTCGCCTCGTGGAAGCGACGCAGCAATGCTGGGATCACATGCGAATTACTCGGGTGGAAGTTGTCGTGTGGCCCGTACAGGTTGGTCGGCATGACTGAGCGGTAATCGCGTC  >  NC_000913/2126257‑2126727

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: