Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,419 G→A 5.7% intergenic (+58/‑150) 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,4190GA 5.7% 715.9 / 21.0 229intergenic (+58/‑150)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (110/105);  new base A (4/9);  total (115/114)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 241 total aligned reads displayed.

GCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAA  >  NC_002947/4901271‑4901567
                                                                                                                                                     |                                                                                                                                                       
gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                        <  1:4429661/150‑1 (MQ=255)
    gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                    <  2:2543818/150‑1 (MQ=255)
       tGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTccgc                                                                                                                                                 <  1:4202746/150‑3 (MQ=255)
       cttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑cc                                                                                                                                              >  2:3267893‑M1/95‑150 (MQ=255)
        gtgaagggtgtgggggcagattgccagcgactacaccgggtgggtttatacgggttaaacaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                             >  2:149350‑M1/94‑150 (MQ=255)
        gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                             <  2:4715657/150‑1 (MQ=255)
         tgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                                            >  1:3793189‑M1/93‑150 (MQ=255)
          ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                                           <  2:265099/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                                          <  1:2153495/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                                          <  1:4098004/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                                          <  1:65972/150‑1 (MQ=255)
           aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttccgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                                          >  1:5471240‑M1/91‑150 (MQ=255)
                tgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGc                                                                                                                                     >  2:2430608‑M1/86‑150 (MQ=255)
                 gttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                                               <  2:3909073‑M1/55‑1 (MQ=255)
                 gttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                                               >  1:3909073‑M1/85‑139 (MQ=255)
                  ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   >  1:1751165‑M1/84‑150 (MQ=255)
                  ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   >  2:5430303‑M1/84‑150 (MQ=255)
                  ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcg                                                                                                                                                         <  1:26384/131‑1 (MQ=255)
                  ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   <  2:1036332/150‑1 (MQ=255)
                  ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   <  2:91710/150‑1 (MQ=255)
                  ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCATGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   <  2:2483579/150‑1 (MQ=255)
                   cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                  <  1:1334292/150‑1 (MQ=255)
                    ggggcagattgccagcgactataccgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                 >  2:5621521‑M1/82‑150 (MQ=255)
                    acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                 >  2:1932340‑M2/1‑149 (MQ=255)
                     caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                   >  1:2410518/1‑147 (MQ=255)
                     caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                                                >  2:3558310‑M2/1‑148 (MQ=255)
                      ggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                                               >  1:3044594‑M1/80‑150 (MQ=255)
                       agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta                                                                                                                              >  1:935274‑M2/1‑146 (MQ=255)
                        gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGGCGCGTCTCGGGTGAATACCCGCCGACAGGCGCGGCCCCGGCGCGGGCC‑‑‑CCCCGAGGGGCGGGgtttc                                                                                                                             >  1:3316856/1‑148 (MQ=255)
                        gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                             >  2:1241079‑M2/1‑145 (MQ=255)
                          agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttaccc                                                                                                                           >  2:1518654‑M2/1‑143 (MQ=255)
                          agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGAGCGGGtttaccc                                                                                                                           >  2:5266905‑M2/1‑143 (MQ=255)
                            aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                                         >  1:4750169‑M2/1‑141 (MQ=255)
                              ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAAGACGCTCCTACAGGCGCGGACCCGGCGCGGTCC‑‑‑CCCCGAGGAGCGGGtttaaccgcga                                                                                                                       <  1:2768695‑M2/150‑12 (MQ=255)
                               ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCTGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                      >  1:3728324‑M1/71‑150 (MQ=255)
                               cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCATCGCGGGTAAATCCGCTCCGACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                                                      <  1:3249290‑M2/150‑13 (MQ=255)
                                aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATACCCTCCTACACGCGCCGGCCCTGCGCGGTGC‑‑‑CCCGTGGGGGCGGGtttaaccccgaag                                                                                                                     <  2:1394378‑M2/150‑14 (MQ=255)
                                 cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaga                                                                                                                    >  2:2876889‑M2/1‑136 (MQ=255)
                                  gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCTTAGGAGCGGGtttacccgcgaagag                                                                                                                   >  1:4362503‑M2/1‑135 (MQ=255)
                                  gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaagag                                                                                                                   >  1:5179891‑M2/1‑135 (MQ=255)
                                      aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCACCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                                               <  2:1642503‑M2/150‑20 (MQ=255)
                                             ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                                                        <  2:3633964/150‑1 (MQ=255)
                                                  tctgttctCGGTGTGGGAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGACGAGGCCGGCACAGGCAGaa                                                                                                   <  2:3310152/150‑1 (MQ=255)
                                                    tgttctCGCTGTGGGAATGGGTTGTTTCCCACAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGc                                                                                                 <  2:388118/150‑1 (MQ=255)
                                                      ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAATCCGCGCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                                               >  1:3722574/1‑150 (MQ=255)
                                                      ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                                               <  2:4232285/150‑1 (MQ=255)
                                                      ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCGCCTACAGGCGCGGCCCGTCCGCGGTCC‑‑‑CCCGGAGGAGCGGGTTTTACCGCGACGCGGCCGGCACAGGCAGAAGCag                                                                                               >  2:3475896/1‑149 (MQ=255)
                                                             tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCCGGTAAATACGCTCCTACCGGCGCGGACCGTGCGCGGGCC‑‑‑CCCCGAGGAGCGGGTTTACCCCCGACGGGGCCGGCCCAGGCCCAAGCAAAACCaaa                                                                                        >  2:2369917/1‑150 (MQ=255)
                                                                   aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCCGAAGCATAACCAAAGCAttt                                                                                  >  2:4312003/1‑150 (MQ=255)
                                                                      gggTTGTGTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGGGGGTAAATCCGCTCCTACAGGCGCGGACGGTGCGCGGTCC‑‑‑CCCGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                                               <  2:2364087/150‑1 (MQ=39)
                                                                       ggTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                                              <  2:3401664/150‑1 (MQ=255)
                                                                       ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCAACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAGAAGCATAAACAAAGCATTTTTTg                                                                              >  1:5070987/1‑150 (MQ=255)
                                                                        gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCAATTCGCGGGTAAATACGCTACTAAAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAACAAAGCATTTTTTGt                                                                             >  2:776701/1‑150 (MQ=38)
                                                                                cccaaGCGTGTTCGCGCCGGCCCCGGCCTCTGCGGGGGTAAATCCGCTCCAACTGGCGCGGACGGTGGGCGGTCC‑‑‑CCGTTAGGAGCGGGTTTCCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                                                     <  2:530059/147‑1 (MQ=37)
                                                                                ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCGTTAAc                                                                     >  2:4089843/1‑150 (MQ=255)
                                                                                ccAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGGAAATACGCGCCAACAGGCGCGGGCCCGGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCCAAGCAATTTTTGTTCTTTAAc                                                                     >  2:570898/1‑150 (MQ=255)
                                                                                ccAAAGCGTGATACCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGACC‑‑‑CCCGTAGGAGCGTTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                                                     <  2:2079380/150‑1 (MQ=255)
                                                                                 cAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGGGCGGGTTTTCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTGTTTTAcc                                                                    >  2:2949662/1‑150 (MQ=39)
                                                                                         gATTGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCCGGCGCGGACCGTGCGCGGCCC‑‑‑CCCGTTGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                                            <  2:1446989/150‑1 (MQ=255)
                                                                                              cgcgTGCGCGGGCCTCTTGGGGGTAAATCCCGCCCCTACCGGCGGGGACGGTGCGGGGCCC‑‑‑CCCGTCGGAGGGTTTTTACCCGCTAAGCGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccac                                                       <  1:2414019/150‑1 (MQ=255)
                                                                                                  tGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                             >  1:2161661‑M2/1‑71 (MQ=255)
                                                                                                  tGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                             <  2:2161661‑M2/76‑6 (MQ=255)
                                                                                                           CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGCCCGTGCGGGGTCC‑‑‑CCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                          >  1:4938942‑M1/1‑150 (MQ=255)
                                                                                                           CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGACGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                          <  2:5145544‑M1/150‑1 (MQ=255)
                                                                                                            TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                          >  2:5236745‑M1/1‑149 (MQ=255)
                                                                                                            TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                          <  1:5237134‑M1/149‑1 (MQ=255)
                                                                                                              TTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg                                       <  1:102031‑M1/150‑1 (MQ=255)
                                                                                                                  CGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt                                   <  2:4358658‑M1/150‑1 (MQ=255)
                                                                                                                   GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc                                  <  1:3323375‑M1/150‑1 (MQ=255)
                                                                                                                   GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc                                  <  1:1559613‑M1/150‑1 (MQ=255)
                                                                                                                     GTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg                                <  2:2797881‑M1/150‑1 (MQ=255)
                                                                                                                      TAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCATTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc                               <  2:4593551‑M1/150‑1 (MQ=255)
                                                                                                                       AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                              <  1:3628743‑M1/150‑1 (MQ=255)
                                                                                                                        aaTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcaaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttg                             <  1:4658208‑M2/150‑102 (MQ=255)
                                                                                                                        AATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTCTAGTCGCTGGCACTCTGccc                             <  2:5017562‑M1/150‑1 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaggggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                           <  2:3613774‑M2/150‑104 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTCTAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  2:3955130/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAGCCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  1:3833705/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  1:3434289/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  2:3661300/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  1:4916142/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  2:2636957/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  1:86876/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTGATAGTCGCTGGCAATCTGCCCCa                           >  1:2463892/1‑150 (MQ=255)
                                                                                                                          tCCGCGCCTGTAGGAGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                           >  2:5077666/1‑150 (MQ=38)
                                                                                                                            ccGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa                          >  2:2759949/1‑150 (MQ=255)
                                                                                                                              gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                        >  2:1159291/1‑150 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                    <  1:1125280‑M2/150‑111 (MQ=255)
                                                                                                                                   tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCACTCTGCCCCAACACCCtt                   >  2:3992040/1‑150 (MQ=255)
                                                                                                                                   tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattc                                                                          <  1:3955449‑M2/95‑57 (MQ=255)
                                                                                                                                     cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg                 <  1:5180307‑M2/150‑114 (MQ=255)
                                                                                                                                     cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttg                                                                                         <  1:1487984‑M2/78‑42 (MQ=255)
                                                                                                                                      aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggc                <  2:4111200‑M2/150‑115 (MQ=255)
                                                                                                                                       acgggggACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGTGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAg               >  2:1985372/6‑150 (MQ=255)
                                                                                                                                          gcgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc            >  2:2321727/1‑150 (MQ=255)
                                                                                                                                           gggaCCGCGCACGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatca           <  1:181878‑M2/149‑120 (MQ=255)
                                                                                                                                            ggACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACaa          >  2:3067838/1‑150 (MQ=255)
                                                                                                                                                cGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctg      <  1:1850738‑M2/150‑125 (MQ=255)
                                                                                                                                                  tGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGc    >  1:2944815/1‑150 (MQ=255)
                                                                                                                                                    cgcgGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCaa  >  2:614173/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAA  >  NC_002947/4901271‑4901567

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: