Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,393 T→C 8.8% intergenic (+32/‑176) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,901,3930TC 8.8% 556.2 / 30.5 196intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (72/103);  new base C (10/7);  total (85/111)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

Only 100 of 218 total aligned reads displayed.

CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4901244‑4901542
                                                                                                                                                     |                                                                                                                                                         
cAATGCCATGCGTGGGGAGCCGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                           >  2:3754685/1‑150 (MQ=255)
    gCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                       <  1:1216588/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                     <  2:697696/150‑1 (MQ=255)
                 agcagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGcgcg                                                                                                                                         <  2:4551075/150‑1 (MQ=255)
                   cagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                       <  1:1548375/150‑1 (MQ=255)
                    agcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGAc                                                                                                                                      <  2:3344915/150‑1 (MQ=255)
                    agcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCGCCTGTAGGAGCGGat                                                                                                                                      <  2:1576912/150‑2 (MQ=255)
                        gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTg                                                                                                                                  <  2:700267/150‑1 (MQ=255)
                         cagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgc                                                                                                                                 <  2:545281/150‑1 (MQ=255)
                           gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                               <  1:4429661/150‑1 (MQ=255)
                            caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGGCGCATGGCGGGTGAATCCGCGCCGACAGGCGCGGACCGGGcgcg                                                                                                                              >  1:1446989/1‑150 (MQ=255)
                            caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                              <  2:4126091/150‑1 (MQ=255)
                               gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCt                                                                                                                                                      >  2:2644313/1‑124 (MQ=255)
                               gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                           <  2:3569456/150‑1 (MQ=255)
                               gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                           <  2:2544032/150‑1 (MQ=255)
                                  cttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑cc                                                                                                                     >  2:3267893‑M1/95‑150 (MQ=255)
                                   gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCTTGCGCGGTcc‑‑‑ccc                                                                                                                    <  2:2820034/150‑1 (MQ=255)
                                   gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                    <  2:4715657/150‑1 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                   >  2:5289230/1‑150 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                   <  1:2410627/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCTCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:4098230/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:3316948/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:2153502/150‑1 (MQ=255)
                                      aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 >  1:942572‑M1/91‑150 (MQ=255)
                                      aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttccgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 >  1:5471240‑M1/91‑150 (MQ=255)
                                         ggtgttggggcagattgccagctactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                              >  2:3417721‑M1/88‑150 (MQ=255)
                                         ggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttccgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                              >  1:3263481‑M1/88‑150 (MQ=255)
                                         ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                              <  1:2125727/150‑1 (MQ=255)
                                         ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                              <  1:5435490/150‑1 (MQ=255)
                                            gttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                      >  1:3909073‑M1/85‑139 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          >  1:4272084‑M1/84‑150 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          >  1:1751165‑M1/84‑150 (MQ=255)
                                             gtggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          >  1:2428963‑M1/84‑150 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcg                                                                                                                                <  1:26384/131‑1 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  2:1036332/150‑1 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  1:3675170/150‑1 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCATGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  2:2483579/150‑1 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  2:5621611‑M1/82‑150 (MQ=255)
                                               acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTCGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATACGCTCCTAAAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        <  2:3934838‑M2/150‑2 (MQ=255)
                                                gggcagattgccagcgactataacgggtggggttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                       >  2:1448284‑M1/81‑150 (MQ=255)
                                                caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  2:2410518/147‑1 (MQ=255)
                                                   gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                    >  2:1241079‑M2/1‑145 (MQ=255)
                                                   gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                    >  2:5356861‑M2/1‑145 (MQ=255)
                                                   gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttac                                                                                                    >  1:5118224‑M2/1‑145 (MQ=255)
                                                   gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttac                                                                                                    >  2:1900120‑M2/1‑145 (MQ=255)
                                                   gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCCGGCGCGGGCC‑‑‑CCCGGAGGAGCGGGtttac                                                                                                    <  1:393964‑M2/150‑6 (MQ=255)
                                                     agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGAGCGGGtttaccc                                                                                                  >  2:5266905‑M2/1‑143 (MQ=255)
                                                       aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                >  1:4750169‑M2/1‑141 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATACGCGACTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgc                                                                                              <  2:61374‑M2/150‑12 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAAGACGCTCCTACAGGCGCGGACCCGGCGCGGTCC‑‑‑CCCCGAGGAGCGGGtttaaccgcga                                                                                              <  1:2768695‑M2/150‑12 (MQ=255)
                                                          cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  2:3545629‑M2/1‑138 (MQ=255)
                                                           aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATACCCTCCTACACGCGCCGGCCCTGCGCGGTGC‑‑‑CCCGTGGGGGCGGGtttaaccccgaag                                                                                            <  2:1394378‑M2/150‑14 (MQ=255)
                                                            cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaga                                                                                           >  2:3908094‑M2/1‑136 (MQ=255)
                                                                 aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCACCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                      <  2:1642503‑M2/150‑20 (MQ=255)
                                                                 aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTAGCGGGTAAAAACCCACCAACAAGCCCCGGCCCGGGCCGGGGCC‑‑‑CCCGAGGAGCCGGGTTAACCCCCAAGAGGccc                                                                                      >  1:3504116/1‑149 (MQ=255)
                                                                        ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                               <  2:3633964/150‑1 (MQ=255)
                                                                         ggtttctgttcgCGCTGTGGGATGGGGTTGTTTCCCATAGGGTGTTAGCGCCTGCGCCGGCCTCTTCGGGGGTAAATCCGCTCCTACCGGGGGGGACCGTGCGCGGCCC‑‑‑CCCGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGc                                                                              <  2:1497362/147‑1 (MQ=25)
                                                                                gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCATCCCCGGGAAATCCCCACCAACACGCGCCGACCGGGCGCCGTCC‑‑‑CCCCGAGGAGCGGGTTTACCCCCGAAGAGGCCCGCACCGGCACAAGCa                                                                       >  1:2293446/1‑150 (MQ=255)
                                                                                 ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                      <  2:4232285/150‑1 (MQ=255)
                                                                                       cTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaa                                                                <  2:2797689/150‑1 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCgcgca                                                                                                                                <  2:2348637/83‑2 (MQ=255)
                                                                                              aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCCGAAGCATAACCAAAGCAttt                                                         >  2:4312003/1‑150 (MQ=255)
                                                                                                 gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCGGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGTACGGTGCGCGGTCC‑‑‑CCCTTTGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                      <  2:2364286/150‑1 (MQ=38)
                                                                                                  ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                     <  2:4844493/150‑1 (MQ=255)
                                                                                                  ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCCCTCCTACAGGCGCGGCCCGTGCGCGGGCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCCAAGAGGCCGGCACAGGCAGAAGCATAAACAAAGCATTTTTTg                                                     >  1:5070762/1‑150 (MQ=38)
                                                                                                   gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCAATTCGCGGGTAAATACGCTACTAAAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAACAAAGCATTTTTTGt                                                    >  2:776701/1‑150 (MQ=38)
                                                                                                         tCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAGAAGCATAACCCAAGCATTTTTGGTTCTTTa                                              >  2:3872170/1‑150 (MQ=38)
                                                                                                           ccAAAGCGTGATACCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGACC‑‑‑CCCGTAGGAGCGTTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                            <  2:2079380/150‑1 (MQ=255)
                                                                                                            cAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGGGCGGGTTTTCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTGTTTTAcc                                           >  2:2949662/1‑150 (MQ=39)
                                                                                                             aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                                          <  2:905698/150‑1 (MQ=255)
                                                                                                                gCGTGATAGCGCCTGCGCGGGCCTGTTCGGGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGta                                       <  2:275223/150‑1 (MQ=255)
                                                                                                                    gATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                   <  1:3754685/150‑1 (MQ=255)
                                                                                                                             tGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                    <  2:2161661‑M2/76‑6 (MQ=255)
                                                                                                                               cgcCGGCCTCTTCGGGGTTAAATCCGCCCCTACCGGCGGGGACGGTGGGGGGCCC‑‑‑CCCGTCGGCGCGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                        <  1:2915096/150‑1 (MQ=255)
                                                                                                                                      CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGACGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 <  2:5145544‑M1/150‑1 (MQ=255)
                                                                                                                                       TCTTCGCGTTTAAATCCGCTCCTACCGGCGCGGCCGGGGCGCGGTCC‑‑‑CCGTTAGGCGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                >  1:696308‑M1/1‑150 (MQ=255)
                                                                                                                                       TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 <  1:5236745‑M1/149‑1 (MQ=255)
                                                                                                                                       TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 >  2:5236745‑M1/1‑149 (MQ=255)
                                                                                                                                       TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                <  1:913926‑M1/150‑1 (MQ=255)
                                                                                                                                         TTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg              <  1:102031‑M1/150‑1 (MQ=255)
                                                                                                                                         TTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCC‑‑‑GCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTATTTGTTCTTTAACCCGTATAAACCCACCCGGTATAGTCGCAggcta           <  2:407765‑M1/150‑3 (MQ=255)
                                                                                                                                             CGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt          <  2:4358658‑M1/150‑1 (MQ=255)
                                                                                                                                              GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc         <  1:1559613‑M1/150‑1 (MQ=255)
                                                                                                                                               GGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt        <  1:3211262‑M1/150‑1 (MQ=255)
                                                                                                                                                  AAATCCGCTCCTACCGGCGCGGACGGTGCGCGGCCC‑‑‑CCCGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc     >  1:3356963‑M1/1‑150 (MQ=255)
                                                                                                                                                   aaTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcaaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttg    <  1:4658208‑M2/150‑102 (MQ=255)
                                                                                                                                                   AATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTCTAGTCGCTGGCACTCTGccc    <  2:5017562‑M1/150‑1 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt   <  1:5521807‑M2/150‑103 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaggggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg  <  2:3613774‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattc                                                 <  2:2359450‑M2/103‑57 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattc                                                 >  1:2359450‑M2/1‑47 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAGCCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:3833705/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:4540227/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:4916142/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:372732/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:2540861/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:3661300/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:2636957/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTATACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAATTTTTGATCTTTAACCCGTATAAACCCACCCGTAATAGTCGCTGGCAATCTGCCCCa  >  2:3223496/1‑150 (MQ=255)
                                                                                                                                                     tCCGCGCCTGTAGGAGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:5077666/1‑150 (MQ=38)
                                                                                                                                                     |                                                                                                                                                         
CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4901244‑4901542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: