Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,366 T→C 13.9% 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,898,3660TC13.9% 311.1 / 31.8 125intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/66);  new base C (10/7);  total (50/75)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 139 total aligned reads displayed.

CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4898217‑4898515
                                                                                                                                                     |                                                                                                                                                         
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGTAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                           <  2:2298955/150‑1 (MQ=255)
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                           <  1:839262/150‑1 (MQ=255)
  atgcgatgcgTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATcc                                                                                                                                                         <  1:3350579/150‑1 (MQ=255)
    tgaagagttttttcacgttgcgtccttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                       >  2:1160287‑M1/125‑150 (MQ=255)
    gCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                       <  1:448599/150‑1 (MQ=255)
     ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCt                                                                                                                                                      <  1:5164783/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                     <  1:107666/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                     <  1:420951/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATACGCTc                                                                                                                                                     >  2:4318477/1‑150 (MQ=255)
       agagttttttcacgttgcgtccttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCATCGCGGGTAAAAccgcacc                                                                                                                                                    <  1:4318477‑M1/29‑4 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                                  <  1:1553821/150‑1 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                                  <  1:655708/150‑1 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                                  <  1:1553529/150‑1 (MQ=255)
           ttttttcacgttgcgtccttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACa                                                                                                                                                >  2:1254383‑M1/118‑150 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                               <  1:1701447/150‑1 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                               <  1:4106149/150‑1 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                               <  2:4954897/150‑1 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                               <  2:5217595/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                           <  1:4220759/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                           <  1:2626987/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                           <  1:2244189/150‑1 (MQ=255)
                 agcagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGcgcg                                                                                                                                          <  1:1778163/150‑1 (MQ=255)
                       tGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCGCCTGTAGGAGCGGattta                                                                                                                                    <  1:3312892/150‑5 (MQ=255)
                        gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGt                                                                                                                                    <  2:1855211/149‑1 (MQ=255)
                        gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGt                                                                                                                                    >  1:1855173/1‑149 (MQ=255)
                        gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGt                                                                                                                                    >  1:1855155/1‑149 (MQ=255)
                        gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGt                                                                                                                                    <  2:1855173/149‑1 (MQ=255)
                            caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                               <  1:469937/150‑1 (MQ=255)
                                  tGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGTCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑cc                                                                                                                      <  2:2042915/150‑1 (MQ=255)
                                   ttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                     >  1:3093394‑M1/94‑150 (MQ=255)
                                    tgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                    >  1:1369318‑M1/93‑150 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                    <  1:1902723/150‑1 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                    <  2:1509902/150‑1 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                    <  2:1509870/150‑1 (MQ=255)
                                    cccGACCAGCCACAAGTAGACCCACGAGCAGGGGGCGTCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccgc                                                                                                                    <  2:2361768/150‑3 (MQ=255)
                                     gaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                   >  2:532783‑M1/92‑150 (MQ=255)
                                     gaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                   >  2:532967‑M1/92‑150 (MQ=255)
                                     ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGCTCTGTTCTCGCTGCGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                   <  1:125040/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                  <  2:4691087/150‑1 (MQ=255)
                                      cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                  <  1:3618952/150‑1 (MQ=255)
                                       gACCTGCCACACGAAGCCCCACGAGCAGGGGGCGGCGGTCTGTTCTCGCTGTGGGAAGGGGTTGTGTCCCCCAGCGGGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                 <  2:4744572/150‑1 (MQ=255)
                                       gACCCGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                 <  2:3520329/150‑1 (MQ=255)
                                        aCCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAgg                                                                                                                >  2:845726/1‑150 (MQ=255)
                                         ggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCg                                                                                                           >  1:1047075‑M1/88‑150 (MQ=255)
                                         ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                               <  1:5138037/150‑1 (MQ=255)
                                           gggtggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgctactgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGCG‑‑GACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                         >  1:4803599‑M1/86‑150 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          >  1:5267800‑M1/84‑150 (MQ=255)
                                             ccACAAGCCGACCCACGCGCAGGGGGCGGCGGTCTGTTCTCGCTGTGGGAATGGGGTGTTTCCCACAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                          <  2:4748283/150‑1 (MQ=38)
                                              cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                         <  2:328063/150‑1 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                        >  2:443935‑M1/82‑150 (MQ=255)
                                                gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                     >  1:2717522‑M1/81‑137 (MQ=255)
                                                caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                       <  2:2812169‑M2/150‑3 (MQ=255)
                                                       aCCCCCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                >  1:1885819‑M2/1‑141 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                            >  1:1478131/1‑124 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcga                                                                                              <  2:2546244‑M2/150‑12 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcga                                                                                              >  1:2615930‑M2/1‑139 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTAGCGGGTAAAACCGCTCCGACAGGCGCGGGCCCGGCGCGGGCC‑‑‑CCCCGTGGGGCGGGtttacccccgc                                                                                              <  1:1194294‑M2/150‑12 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCACGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACGTCGCGGGGAAATACGCCCCCAACGGCGCCGCCCGGGCCCGGGCC‑‑‑CCCCGAGGCGCCGGGTTTcccccca                                                                                               >  1:2880675/1‑147 (MQ=255)
                                                          cccgaGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  1:4077862‑M2/3‑138 (MQ=255)
                                                          cACGCGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  1:4282123‑M2/1‑138 (MQ=255)
                                                          cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                             >  2:2690981‑M2/1‑138 (MQ=255)
                                                                        ggcgTTCTGTTCTCGCTGTGGGATTGGGTTGTTTCCCAAAGCGTGATCGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGACC‑‑‑CCCGTCGCAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                               <  2:4244561/150‑1 (MQ=39)
                                                                        ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCACTGGGTGATAGCGCCTGCGCGGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACGGTGCGCGGTCC‑‑‑CCGGTCGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAgg                                                                               <  2:3189597/150‑1 (MQ=255)
                                                                                gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCCCACCTACAGGCGCGGGCCGGGCGCGGTCC‑‑‑CCCATAGGAGCGGGGTTACCCGCGCAGAGGCCGGCCCAGGCACAAGCa                                                                       >  1:1745526/1‑150 (MQ=255)
                                                                                  tctcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATa                                                                     <  2:4991122/150‑1 (MQ=255)
                                                                                    tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTAGCGGGTAAATGCCCACCAACACGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGATGAGCGGGGTTACCCCCGAATAGGCCGGCACAGGCAGAAGCAAAAc                                                                   >  1:1217192/1‑150 (MQ=255)
                                                                                            ggAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCg                                                                                                                                     <  2:1409697/80‑1 (MQ=255)
                                                                                            ggAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCg                                                                                                                                     >  1:1409697/1‑80 (MQ=255)
                                                                                              aaTGGGTTTTTTCCCAAAGCGTTATCGCGCCTGCGCCGGCCTCTTGGCGGGTAATTCCGCTCCCACAGGGGCGGACCGTGCGCGGTCC‑‑‑CCCGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAttt                                                         <  1:1753877/150‑1 (MQ=38)
                                                                                               aacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTa                                                                                                     <  2:4325891‑M1/72‑1 (MQ=255)
                                                                                               aacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTa                                                                                                     >  1:4325891‑M1/34‑105 (MQ=255)
                                                                                               aTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCGCCTACACGCGCGGACCCGGCGCGGGCC‑‑‑CCCCGAGGAGCGGGTTTAACCCCGAAGAGGGCGGCACCGGCAGAAGCAAAACCAAAGCaaaat                                                        >  2:4119711/1‑146 (MQ=255)
                                                                                                   gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCAACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTACCCCCGAAGACGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGt                                                    >  2:2782166/1‑150 (MQ=255)
                                                                                                   gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATACGCTCCGACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGCGGGGTTACCCGCGAAGAGGCCGGCACAGGCAGACGCCTAACCAAAGCATTTTTTGt                                                    >  1:4893257/1‑150 (MQ=255)
                                                                                                        tccccAACGGGTGATGGCGCCGGCGCCGGCCTTTTCGGGGGTAAATCCGCACCTCCCGGCGCGGCGGGTGGGGGGTCC‑‑‑CCCTGGGGTGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCttt                                               <  1:1029841/148‑1 (MQ=255)
                                                                                                           ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                               >  2:147618/1‑84 (MQ=255)
                                                                                                           ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                               <  1:147618/84‑1 (MQ=255)
                                                                                                            cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                           >  2:1761434/1‑150 (MQ=255)
                                                                                                             aacgCGTTATAGCGCCTGCGCCGGCCTCTTGGCGGGTAAATCCGCTCCTCCCGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                                          <  1:1790737/147‑1 (MQ=38)
                                                                                                                                 CCGGCCTCGTCGCGGGTAATTCCGCTCCTACTGGCGCGGACCGTGGGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                      >  2:2483159‑M1/1‑150 (MQ=255)
                                                                                                                                  CGGCCTGTTCGGGGGTAAATCCGCTCCTACCGGCGCGGACGGTGCGGGGGCC‑‑‑CCCGTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAg                     >  2:1820372‑M1/1‑150 (MQ=255)
                                                                                                                                  CGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAg                     <  2:2697496‑M1/150‑1 (MQ=255)
                                                                                                                                    GCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                   <  1:1854095‑M1/55‑1 (MQ=255)
                                                                                                                                    GCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                   >  2:1854095‑M1/1‑55 (MQ=255)
                                                                                                                                     CCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACTGGCAGTAGCATAACCAAAGCATTTTTTGTTCTTTTACCGGTGTAATGCCTCGCGGTATAGTCg                  <  2:3030427‑M1/150‑1 (MQ=255)
                                                                                                                                       TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                <  1:5113035‑M1/150‑1 (MQ=255)
                                                                                                                                         TTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg              <  1:2860458‑M1/150‑1 (MQ=255)
                                                                                                                                              GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTCTAGTCGCTGGCAATc         <  2:4548723‑M1/150‑1 (MQ=255)
                                                                                                                                               GGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt        <  1:3233007‑M1/150‑1 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg  <  2:1130364‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttc       >  2:1509335‑M2/1‑47 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttc       <  1:1509335‑M2/145‑99 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcga                                        >  2:3703762‑M2/1‑47 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcga                                        <  1:3703762‑M2/112‑66 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc                                                                                    <  2:3056948‑M2/68‑22 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc                                                                                    >  1:3056948‑M2/1‑47 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGTATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:710213/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:3948766/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:3376953/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg  <  2:2169307‑M2/150‑104 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4898217‑4898515

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: