Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,030,403 T→G 9.6% intergenic (+32/+35) PP_5729 → / ← rpmG hypothetical protein/50S ribosomal protein L33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,030,4030TG 9.6% 113.5 / 4.9 43intergenic (+32/+35)PP_5729/rpmGhypothetical protein/50S ribosomal protein L33
Reads supporting (aligned to +/- strand):  ref base T (16/22);  new base G (2/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ACCGGCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGGTGTAGAA  >  NC_002947/6030255‑6030548
                                                                                                                                                     |                                                                                                                                                 
aCCGGCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCtt                                                                                                                                                  >  2:642735/1‑150 (MQ=255)
  cGGCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTc                                                                                                                                                <  1:226796/150‑1 (MQ=255)
     ccGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGa                                                                                                                                             >  1:162967/1‑150 (MQ=255)
        cACCCGGGTGGTCGGCCCCAAAGAAATGTGCACCCTCTAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAACACCCGCCTCCTTTCGGATGc                                                                                                                                          <  1:311497/150‑1 (MQ=255)
         aCCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCg                                                                                                                                         <  2:1053217/150‑1 (MQ=255)
           ccGAGTGGACGGCCCCAACGAAACGCGCACCCGCGAGTGCGAGGTCGAGCGTAACCAGAACAATGCCGGTGCCGACCTGCTGATCCAGGGCTTCACCTTCAACTGCCACAGAAAAACAA‑AAAAAAAAACCCGCCTCCTTTCGGATGCggg                                                                                                                                       <  2:676731/150‑1 (MQ=255)
             gggTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGtt                                                                                                                                     <  2:877662/150‑1 (MQ=255)
                   tCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTc                                                                                                                               >  1:1122234/1‑150 (MQ=255)
                    cGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTCCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCa                                                                                                                              <  2:402410/150‑1 (MQ=255)
                     ggCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAg                                                                                                                             >  1:1133045/1‑150 (MQ=255)
                      gCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAgg                                                                                                                            <  2:1133045/150‑1 (MQ=255)
                               aaaTGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCATCCGAAAGGAGGCGGGTTTTTTTCAGGTCGCCGATc                                                                                                                   <  2:1144242/150‑1 (MQ=39)
                                  tgtgCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAAt                                                                                                                <  1:923968/150‑1 (MQ=255)
                                       aCCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACtt                                                                                                           <  2:1109458/150‑1 (MQ=255)
                                       aCCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACtt                                                                                                           >  2:1051209/1‑150 (MQ=255)
                                            cgAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCt                                                                                                      <  2:1097311/150‑1 (MQ=255)
                                                          aGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTac                                                                                        <  1:344064/150‑1 (MQ=255)
                                                                tCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCATCCGAAAGGAGGCGGGtttttttt                                                                                                                               <  1:1194815/105‑2 (MQ=25)
                                                                tCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCATCCGAAAGGAGGCGGGtttttttt                                                                                                                               >  2:1194815/1‑104 (MQ=25)
                                                                       cAAGGCCGCTGCCGCCCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTACAACAGCAAAAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTg                                                                            <  2:1030009/150‑1 (MQ=255)
                                                                          ggCCGGTGCCGACCTGCTGATCCAGGGCCTCGACTCCAACTGACCCTGCGACGGCG‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGaa                                                                        >  1:1056653/1‑150 (MQ=255)
                                                                             cGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAACCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGaacaa                                                                     <  1:1051209/150‑1 (MQ=255)
                                                                                   cGACATACTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑ACAAAATGTACAACTTGCCCGTGGAGGGTGTTCCTTTTCAGTGCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGAt                                                               <  1:274558/150‑1 (MQ=255)
                                                                                             acccAGGGCGTCAACTTCAACTGACCCTGCAACAGCT‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTTGTTGATCCCTTACTTCATCATTGCTTACTTATTCACCACGTTCGTTCTAACTACCGGAGCATATGTCtt                                                     <  2:12488/148‑1 (MQ=255)
                                                                                                     tttCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGa                                             <  2:1061712/149‑1 (MQ=255)
                                                                                                      ttCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTAc                                                                                                             <  1:137889/85‑1 (MQ=255)
                                                                                                      ttCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTAc                                                                                                             >  2:137889/1‑85 (MQ=255)
                                                                                                          aCTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTg                                        >  1:305128/1‑150 (MQ=255)
                                                                                                            ttCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCtt                                                                                            <  2:778838/96‑1 (MQ=255)
                                                                                                            ttCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCtt                                                                                            >  1:778838/1‑96 (MQ=255)
                                                                                                               aaCTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCgg                                   >  1:1019412/1‑150 (MQ=255)
                                                                                                                   gACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGggtg                               <  2:162967/150‑1 (MQ=255)
                                                                                                                      ccTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGgtggtg                            >  1:1062915/1‑150 (MQ=255)
                                                                                                                        tGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTgcg                          >  2:1061387/1‑150 (MQ=255)
                                                                                                                                 c‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCtt                 <  2:1019412/150‑1 (MQ=255)
                                                                                                                                   tgaaaaaaaCCCGCATCCGAAAGGAGGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTACGGGGGGGTGCGCTTGTTCttt                >  2:160811/3‑149 (MQ=38)
                                                                                                                                   aaaaaaaaaCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCCATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGGTCATATTTCTTGATCTCGATTTTGTCCGGCGTGGTGCGCTTGTTCTTg                >  1:657214/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaaaCCCGCCTCCTTTCGGAAGCGGGGTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTg                >  1:1125468/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCgg            >  1:332594/1‑150 (MQ=255)
                                                                                                                                        aaaaCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCggt           <  2:1122234/150‑1 (MQ=255)
                                                                                                                                           aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCggtggt        <  2:1062915/150‑1 (MQ=255)
                                                                                                                                             ccGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGggtg                               >  1:1150321/1‑125 (MQ=255)
                                                                                                                                             ccGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGggtg                               <  2:1150321/125‑1 (MQ=255)
                                                                                                                                             ccGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTAATCTCGTTTTTCTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGgtgt      >  1:1172291/1‑150 (MQ=255)
                                                                                                                                                 ctcctTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCatat                                                          <  2:375018/94‑1 (MQ=255)
                                                                                                                                                 ctcctTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCatat                                                          >  1:375018/1‑94 (MQ=255)
                                                                                                                                                 ctcctTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGGTGTAGaa  <  2:305128/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCGGCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGGTGTAGAA  >  NC_002947/6030255‑6030548

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