Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,393 T→C 9.6% 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 9.6% 582.5 / 35.9 201intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (77/102);  new base C (11/8);  total (90/111)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 215 total aligned reads displayed.

CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4901244‑4901542
                                                                                                                                                     |                                                                                                                                                         
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCTCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                           <  2:1158337/150‑1 (MQ=255)
 aaTGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATc                                                                                                                                                          <  1:2862110/150‑1 (MQ=255)
 aaTGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATc                                                                                                                                                          <  2:2084714/150‑1 (MQ=255)
     ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCt                                                                                                                                                     >  1:7301400/1‑150 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                    <  1:320697/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                    <  1:6081440/150‑1 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                                 <  2:3360753/150‑1 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                              >  2:1623208/1‑150 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                          <  2:5287169/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                          <  2:6158966/150‑1 (MQ=255)
                   cagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                       <  2:295696/150‑1 (MQ=255)
                         cagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgc                                                                                                                                 <  2:1624909/150‑1 (MQ=255)
                          agcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcg                                                                                                                                <  1:5783958/150‑1 (MQ=255)
                                  tGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTccgc                                                                                                                        <  1:6171976/150‑3 (MQ=255)
                                    tgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                   >  1:6677791‑M1/93‑150 (MQ=255)
                                    tgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                   >  1:6958419‑M1/93‑150 (MQ=255)
                                      cGACCTGCCACCAGAAGACCCACGAGCAGGGGGCTGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  2:1419171/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTccgcgcct                                                                                                                    <  1:2403085/150‑7 (MQ=38)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTccgcgcct                                                                                                                    <  2:5022838/150‑7 (MQ=38)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:7196684/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  2:1889272/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:2778523/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  1:1518460/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 <  2:2343095/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGGGAATCCGCGCCGCCCCGCGCGGGCCCGGCGCGGGcc‑‑‑ccccgg                                                                                                                 >  1:3965269/1‑147 (MQ=255)
                                      aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                 >  2:4784842‑M1/91‑150 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                <  1:6388555/150‑1 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                <  1:6248274/150‑1 (MQ=255)
                                        aCCTGCCACAAGAAGACCCACGAGAAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATc                                                                                                                                                          <  1:4165155/111‑1 (MQ=255)
                                           tgttggggcagattgccagcgactctaacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGc                                                                                                            >  2:4741909‑M1/86‑150 (MQ=255)
                                           tgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGc                                                                                                            >  1:1136661‑M1/86‑150 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAATACGCGCCTAAAGGCGCGGACCGGGGGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          >  1:5420674/1‑150 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                              <  2:5840855/146‑1 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  1:5087532/150‑1 (MQ=255)
                                              tggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                         <  2:4350947‑M1/68‑1 (MQ=255)
                                              cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                         <  1:4999977/150‑1 (MQ=255)
                                              cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgga                                                                                                         <  2:261217/150‑2 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  1:21667‑M1/82‑150 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  2:3576001‑M1/82‑150 (MQ=255)
                                               acaAGAAGACCCCCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  2:3043144‑M2/1‑149 (MQ=255)
                                               acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  1:716055‑M2/1‑149 (MQ=255)
                                               acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCCCGTAGCCGGTAAATCCGCACCGACACGCGCGGCCCCGGGGCGGTCC‑‑‑CCCCGAGGAGCgggg                                                                                                        >  1:746982/1‑149 (MQ=255)
                                                gggcagattgccagcgacgataacgggtgggtttatacgggttaaagaacaaaaattgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGagcggttt                                                                                                       >  1:6883933‑M1/81‑147 (MQ=255)
                                                 ggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                      >  2:2856017‑M1/80‑150 (MQ=255)
                                                 ggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                      >  1:4656456‑M1/80‑150 (MQ=255)
                                                  gcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTa                                                                                                     >  2:2072086‑M1/79‑150 (MQ=255)
                                                     gattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGGGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAccc                                                                                                  >  2:6515994‑M1/76‑150 (MQ=255)
                                                     gagtgccagcgactagaacgggggggtttagacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAccc                                                                                                  >  2:6034855‑M1/76‑150 (MQ=255)
                                                      gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccg                                                                                                 >  2:6470108‑M2/1‑142 (MQ=255)
                                                      gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccg                                                                                                 >  1:6023407‑M2/1‑142 (MQ=255)
                                                       ttgccagctactataacgggtggttttataagggttaaagaacaataaatgctttggttatgcttctgcctgtGCCGGCCTCTTCTCGTTTAAATCCGCTCCTACAGGCGCGGACCGTTCGCTGTCC‑‑‑CCCGTAGGAGCGGGTTTACCcgc                                                                                                >  1:560330‑M1/74‑150 (MQ=255)
                                                       aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                >  2:6374643‑M2/1‑141 (MQ=255)
                                                          ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGTGGGTTTACCCGCGaa                                                                                             >  1:4627301‑M1/71‑150 (MQ=255)
                                                          ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                             >  2:605512‑M1/71‑150 (MQ=255)
                                                          ccagcgactataacgggtgggtttagacgggttaaagaacaaaaaatgcttgggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                             >  2:5168595‑M1/71‑150 (MQ=255)
                                                          cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  2:2800850‑M2/1‑138 (MQ=255)
                                                          cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  2:2756432‑M2/1‑138 (MQ=255)
                                                            cgcgactataccgggtggggttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAaga                                                                                           >  1:7176389‑M1/69‑150 (MQ=255)
                                                            cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaaga                                                                                           >  1:2342356‑M2/1‑136 (MQ=255)
                                                             gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaagag                                                                                          >  2:4929532‑M2/1‑135 (MQ=255)
                                                                 aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                      >  2:1082784‑M2/1‑131 (MQ=255)
                                                                      gcggcgTTCTGTTCTCGCTGTGGGATTGGGTTGTTTCCCAACGCGTGATCGCGCCTGCGCGGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGCCCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca                                                                                 <  2:7012349/150‑1 (MQ=38)
                                                                        ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                               >  2:5437963/1‑150 (MQ=255)
                                                                              gtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                             >  1:5319907‑M1/51‑130 (MQ=255)
                                                                              ctgttcgCGCGGTGGGAAGGGGTTGTTTCCCAACGCGTGATCGCGCCGGCGCCGGCCTCTGCGCGGGTAACTCCGCGCCTCCCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCactat                                                                         <  2:3965269/150‑5 (MQ=37)
                                                                                       cTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAgg                                                                                                               <  1:2672038/103‑1 (MQ=255)
                                                                                        tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaagag                                                                                          >  2:5528561‑M2/1‑111 (MQ=255)
                                                                                        tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccccgaagag                                                                                          <  1:5528561‑M2/126‑16 (MQ=255)
                                                                                                 gggTTGTTCCCCAACGCGTGATTGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGGGCGGACCGTGCGGGGTCC‑‑‑CCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                      <  1:3036090/150‑1 (MQ=38)
                                                                                                 caacccattCCCACAGCGTGATAGCGCCTGCGCCGGCCTCATCGCGGGTAAAAACCCACCAACAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCTAGGAGCCGGTTTAACCCCCAAGAGGCCCGCACACGCAGAAGCAAAAACAAAACAattttt                                                      >  1:7266087/8‑150 (MQ=255)
                                                                                                  ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                     <  2:4827378/150‑1 (MQ=255)
                                                                                                        ttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCTCCGACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCGTAACCAAAGCATTTTTTGTTCttt                                               >  1:6739018/1‑150 (MQ=39)
                                                                                                           ccAATGGGTGATCGCGCCGGCGCCGGCCTCTTGGGGGGTAACTCCGCTCCCACCGGCGCGGACCGGGCGCGGCCC‑‑‑CCCGTCGGTGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                            <  1:4802619/150‑1 (MQ=255)
                                                                                                              ccgCGAGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTCCAGGGGCGGACCGGGCGCGGACC‑‑‑CCGGTTGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCg                                         <  1:3691563/148‑1 (MQ=38)
                                                                                                                gCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCACCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAACCGtg                                       >  2:5680763/1‑149 (MQ=255)
                                                                                                                            cTGGGCCGGCCTCTGGGGGTTTAATCCGGCTCCAACGGGCGGGCACGGTGGGGGGGCC‑‑‑CCCGTTGTTGGGGTTTAACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCg                           <  1:1931559/150‑1 (MQ=255)
                                                                                                                                  CGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAg                     <  2:4859445‑M1/150‑1 (MQ=255)
                                                                                                                                   GGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                    <  1:4694978‑M1/54‑1 (MQ=255)
                                                                                                                                      ctctTGGGGGGTAACCCCGCCCCTCCCGGCGCGGCCGGGGCGGGGTCC‑‑‑CCCGTTGGGGGGGGTTTACCCGCGATGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 <  1:1425109/150‑1 (MQ=255)
                                                                                                                                       TCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                <  2:5228260‑M1/150‑1 (MQ=255)
                                                                                                                                          TCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGc             <  1:226756‑M1/150‑1 (MQ=255)
                                                                                                                                              gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtggggct         <  2:4158487‑M2/150‑97 (MQ=255)
                                                                                                                                              GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTCTAAACCCACCCGTTATAGTCGCTGGCAATc         <  2:2124533‑M1/150‑1 (MQ=255)
                                                                                                                                              GGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc         <  2:1454423‑M1/150‑1 (MQ=255)
                                                                                                                                               GGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGACCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt        <  1:4313275‑M1/150‑1 (MQ=255)
                                                                                                                                                GTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg       <  2:5376607‑M1/150‑1 (MQ=255)
                                                                                                                                                GTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg       <  2:4198775‑M1/150‑1 (MQ=255)
                                                                                                                                                GTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg       <  1:4350947‑M1/150‑1 (MQ=255)
                                                                                                                                                  AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc     <  1:820177‑M1/150‑1 (MQ=255)
                                                                                                                                                  AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc     <  1:1247264‑M1/150‑1 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt   <  2:4029423‑M2/150‑103 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt   <  2:4029345‑M2/150‑103 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt   <  1:3002604‑M2/150‑103 (MQ=255)
                                                                                                                                                    ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc   <  2:7105224‑M1/150‑1 (MQ=255)
                                                                                                                                                    ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc   <  2:4682320‑M1/150‑1 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:2233286/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:7238308/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:1585959/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:1393389/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGATGGCAATCTGCCCCa  >  1:1237156/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
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: