Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,392 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,898,3920GA 5.7% 499.4 / 13.8 158intergenic (+58/‑150)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (75/74);  new base A (3/6);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 164 total aligned reads displayed.

GCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGC  >  NC_002947/4898244‑4898538
                                                                                                                                                    |                                                                                                                                                         
gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                          <  1:581225/150‑1 (MQ=255)
gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                          <  2:3874884/150‑1 (MQ=255)
gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                          <  1:3123156/150‑1 (MQ=255)
gcaTGAGTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                          <  2:4069513/150‑1 (MQ=255)
 caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGcgcg                                                                                                                                                         >  1:1848507/1‑150 (MQ=255)
    gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                      <  1:3740188/150‑1 (MQ=255)
    gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACTTCGCGGGTAAATACGCGCCAAAAGGCGCGGACCCGGCGCGGTc                                                                                                                                                      >  1:329407/1‑150 (MQ=255)
     aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTcc                                                                                                                                                     <  1:172815/150‑1 (MQ=255)
     aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTcc                                                                                                                                                     <  1:326604/150‑1 (MQ=255)
     aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTcc                                                                                                                                                     <  1:2584774/150‑1 (MQ=255)
     aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCGCCTACCGGCGCGGACCGGGCGCGGTcc                                                                                                                                                     >  2:3954684/1‑150 (MQ=255)
      ttGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                                                 >  2:2593062/1‑150 (MQ=255)
        gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c‑‑cc                                                                                                                                             <  2:397961/150‑1 (MQ=255)
         cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCg                                                                                                                                           >  1:4445094/1‑150 (MQ=255)
          ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGt                                                                                                                                          >  2:3992869/1‑150 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCTGGGGGCGGCGTTCTGTTCTCGCTGTGGGCATGGGTTGTTTCCCAAAGCGTGATAGCGCGTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTccgcgcct                                                                                                                                               <  1:3312484/150‑7 (MQ=38)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGTCTCTTCGCGGGTAAATCCGCTCCCACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  1:4509600/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTGGTCGCGGGGGAATCCGGGCCTCCACGCGCGGCCCGGGGGCGGGGC‑‑‑C‑‑CCGtg                                                                                                                                         >  1:2334587/1‑149 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  2:4161339/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  1:3414141/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  1:3978804/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  1:3978833/150‑1 (MQ=255)
           cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTa                                                                                                                                         <  2:1325985/150‑1 (MQ=255)
            gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGGTCTCGCTGTGGGAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAg                                                                                                                                        <  2:1865534/150‑1 (MQ=255)
            agggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAg                                                                                                                                        >  1:587378‑M1/90‑150 (MQ=255)
              ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGa                                                                                                                                      <  2:2488502/150‑1 (MQ=255)
                  ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCgg                                                                                                                                 >  2:4382245‑M1/84‑150 (MQ=255)
                  ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCgg                                                                                                                                 >  1:2229850‑M1/84‑150 (MQ=255)
                  ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGctg                                                                                                                                  <  2:2306591/150‑3 (MQ=255)
                    ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGt                                                                                                                               >  1:805386‑M1/82‑150 (MQ=255)
                    acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAAAACCCAACAAAAGGCGCCGACCGTGCGCCGGCC‑‑‑C‑‑CCGGAGGGGCgggg                                                                                                                               >  1:2344988/1‑149 (MQ=255)
                     gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtt                                                                                                                              >  1:549156‑M1/81‑150 (MQ=255)
                     gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtt                                                                                                                              >  1:3769476‑M1/81‑150 (MQ=255)
                     gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtt                                                                                                                              >  1:2987203‑M1/81‑150 (MQ=255)
                      aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGttt                                                                                                                             >  1:1809551‑M2/1‑147 (MQ=255)
                       agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGttta                                                                                                                            >  2:3498316‑M2/1‑146 (MQ=255)
                       agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGttta                                                                                                                            >  2:2373891‑M2/1‑146 (MQ=255)
                        cagattgccagcgactataccgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTAc                                                                                                                           >  1:2574656‑M1/78‑150 (MQ=255)
                         agattgccagcgactagaacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcgtctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTAcc                                                                                                                          >  2:4456995‑M1/77‑150 (MQ=255)
                          agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCCTAGGAGCGGGtttaccc                                                                                                                         <  1:3335109‑M2/150‑8 (MQ=255)
                           gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGC‑‑CTGTAGGAGCGGAtttacccg                                                                                                                        >  1:1549042‑M2/1‑142 (MQ=255)
                           gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGGCGCGGGGGAAAGCCCTCCTAAAGGCGCGGGCCCGGCGCGGGGC‑‑‑C‑‑CCCAAGGAGCGGGGTTTAccc                                                                                                                         >  1:498138/1‑150 (MQ=255)
                              ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcga                                                                                                                     >  2:750041‑M2/1‑139 (MQ=255)
                               cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaa                                                                                                                    >  1:2145809‑M2/1‑138 (MQ=255)
                                    gccgggggCGGCGTTCTGTTCTCGCTGTGGGAAGGGGTTGTTTCCCAACGCGTGATCGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCGCCTGTAGGAGCGGATTTACCCGCGAaga                                                                                                                  <  2:2697040/147‑1 (MQ=25)
                                    gCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggc                                                                                                               >  2:2737997‑M2/1‑133 (MQ=255)
                                      aGGGGGCGGCGTTCTGTTCGCGCTGTGGGATTGGGTTGTGTCCCAAAGCGTGATAGCGCCGGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCGACCGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                                             >  2:3422088‑M2/1‑131 (MQ=255)
                                      aGGGGGCGGCGTTCTGTACTCGCTGTGGGAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                                             >  2:2121444‑M2/1‑131 (MQ=255)
                                        ggggCGGCGTTCTGTACTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTGTCCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGTtttacccgcgaagaggccggc                                                                                                           >  1:679768‑M2/1‑129 (MQ=255)
                                             ggcgTTCTGTTCTCGCTGTGGGATGGGGTTGTTGCCCAACGCTTTATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGGCC‑‑‑C‑‑CCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAgg                                                                                                      <  2:3899705/150‑1 (MQ=38)
                                                       tctcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATa                                                                                            <  2:1075527/150‑1 (MQ=255)
                                                         tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAc                                                                                          >  2:4057247/1‑150 (MQ=255)
                                                            cTGTGGTATTGGGTTTTTTCCCAATGCGTGATAGCGCGTGCGCCGGCCTCTTGGGGGGTAAATCCGCTCCTACAGGCGGGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaa                                                                                       <  1:3727839/150‑1 (MQ=38)
                                                             tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCCTGCGCGGTCC‑‑‑C‑‑CCGGAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                                                      >  2:2034637/1‑150 (MQ=255)
                                                                    aTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGTAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttt                                                                               <  2:3079168/150‑1 (MQ=255)
                                                                            tttCCCAAAGCGAGATAGCGCCTGCGCCGGCCTCTTAGCGGGTAAATACGCTCCAACAGGGGCGGACCGGGCGCGGTCC‑‑‑C‑‑CCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGGACAGGGAGAAGCATAACCAAAGCATTTTTTGTTCtt                                                                       >  1:4030613/1‑150 (MQ=38)
                                                                                ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                                                   <  2:1578168/150‑1 (MQ=255)
                                                                                ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCGGGGAAATCCGCTCCTACAGGCGCGGACCCTGCGCGGTCC‑‑‑C‑‑CCCGAGGAGCGGGGTTACCCCCGAAGAGGCCGGCCCCGGCAGAAGCAGAACCAAAGCAATTTTTGTTCTTTAAc                                                                   >  2:607255/1‑150 (MQ=255)
                                                                                        tGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaa                                                           >  2:4497605/1‑150 (MQ=255)
                                                                                                ccTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccaccc                                                   <  2:2571217/150‑1 (MQ=255)
                                                                                                   gcgcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGtt                                                <  2:2961645/150‑1 (MQ=255)
                                                                                                      CCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGTTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                                             >  2:893870‑M1/1‑150 (MQ=255)
                                                                                                            TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                       <  2:3901804‑M1/150‑1 (MQ=255)
                                                                                                            TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                       <  2:2548209‑M1/150‑1 (MQ=255)
                                                                                                              TTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg                                     <  1:344766‑M1/150‑1 (MQ=255)
                                                                                                                CGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCa                                   <  1:4377219‑M1/150‑1 (MQ=255)
                                                                                                                        AATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc                           >  1:3992869‑M1/1‑150 (MQ=255)
                                                                                                                         aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt                          <  2:2386097‑M2/150‑103 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttc                              >  1:1758912‑M2/1‑47 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                         <  2:2552223‑M2/150‑104 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:4517156/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  1:974157/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  1:3552916/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:224993/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:3852909/1‑150 (MQ=255)
                                                                                                                          tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:3852955/1‑150 (MQ=255)
                                                                                                                            cGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggc                       <  1:4012737‑M2/150‑106 (MQ=255)
                                                                                                                            cGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggc                       <  1:1029145‑M2/150‑106 (MQ=255)
                                                                                                                             gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                      <  1:4106462‑M2/150‑107 (MQ=255)
                                                                                                                             gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                      <  2:1395664‑M2/150‑107 (MQ=255)
                                                                                                                             gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcactatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                      <  2:1657610‑M2/150‑107 (MQ=255)
                                                                                                                             gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                      >  2:3068818/1‑150 (MQ=255)
                                                                                                                             gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgcttt                                                                                        <  2:336026‑M2/84‑41 (MQ=255)
                                                                                                                               tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  1:2792248‑M2/150‑109 (MQ=255)
                                                                                                                               tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  1:2452971‑M2/150‑109 (MQ=255)
                                                                                                                                 cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                  <  2:41815‑M2/150‑111 (MQ=255)
                                                                                                                                  tACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                 >  1:1317583/1‑150 (MQ=255)
                                                                                                                                  tACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                 >  2:228877/1‑150 (MQ=255)
                                                                                                                                    cAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg               <  1:2577287‑M2/150‑114 (MQ=255)
                                                                                                                                     tacgggggACCGCGCACGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa              >  2:3734405/7‑150 (MQ=255)
                                                                                                                                     aGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa              >  1:731513/1‑150 (MQ=255)
                                                                                                                                     aGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa              <  2:329407/150‑1 (MQ=255)
                                                                                                                                     aGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa              >  1:2216188/1‑150 (MQ=255)
                                                                                                                                       gcgcGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa            >  1:1707699/1‑150 (MQ=255)
                                                                                                                                       gcgcGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAAACCACCCGTTATAGTCGCTTGCAATCTGCCCCAACACCCATAAAGa            >  1:141262/1‑150 (MQ=255)
                                                                                                                                       cgggggACCGCGCACGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTGCAAGa            >  2:4282076/5‑150 (MQ=255)
                                                                                                                                        cgcgGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaat           <  1:1058902‑M2/150‑118 (MQ=255)
                                                                                                                                         ggggACCGTGCGGGGTCC‑‑‑C‑‑CCGTTGGAGCGTGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc          <  1:2756056/148‑1 (MQ=255)
                                                                                                                                         ggggACCGTCCGGGGCCC‑‑‑C‑‑CCGTTGGTGCGGTTTGACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc          <  1:4197290/148‑1 (MQ=255)
                                                                                                                                                 tGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGc  >  2:1640380/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                         
GCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑C‑‑CCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGC  >  NC_002947/4898244‑4898538

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: