Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,674,892 A→C 20.0% intergenic (+135/‑26) PP_3236 → / → PP_3237 lipoprotein OprI/universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,674,8920AC20.0% 34.2 / 3.8 20intergenic (+135/‑26)PP_3236/PP_3237lipoprotein OprI/universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base C (2/2);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

TGCAGGCGTTCGCACGCTTGATAAAAACCAAGCGCCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCATCGGCCAG  >  NC_002947/3674765‑3675018
                                                                                                                               |                                                                                                                              
tGCAGGCGTTCGCACGCTTGATAAAAACCAAGCGCCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAggccggcc                                                                                                                        >  2:437799/1‑136 (MQ=255)
                                 gcCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAAc                                                                                       >  1:248266/1‑136 (MQ=255)
                                 gcCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAAc                                                                                       >  1:336762/1‑136 (MQ=255)
                                 gcCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAAc                                                                                       >  1:523510/1‑136 (MQ=255)
                                 gcCAAGGGCGGCATGAACAGGATGCTGGCAGCAGCAGTCAAACACTGCCCGCCGAAAGCCCCCTCTCGAAGCAACCCGAGCGACCTGAGTGAACAGGCCGGCCCACCGTTGGAGGTAATCAAGAGCCACTATCAAc                                                                                       >  1:336763/1‑136 (MQ=255)
                                                    ggATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCCGGCCGGCCCACCGTTGGAGGTAAt                                                                                                        <  1:23260/100‑1 (MQ=255)
                                                    ggATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCCGGCCGGCCCACCGTTGGAGGTAAt                                                                                                        <  1:25990/100‑1 (MQ=255)
                                                    ggATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCCGGCCGGCCCACCGTTGGAGGTAAt                                                                                                        >  2:23260/1‑100 (MQ=255)
                                                    ggATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCCGGCCGGCCCACCGTTGGAGGTAAt                                                                                                        >  2:25990/1‑100 (MQ=255)
                                                                      tCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGc                                                  <  1:385292/136‑1 (MQ=255)
                                                                      tCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGc                                                  <  1:468228/136‑1 (MQ=255)
                                                                             cTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGcc                                           >  2:342602/1‑136 (MQ=255)
                                                                                            ccccTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCAGt                            <  1:342602/136‑1 (MQ=255)
                                                                                             cccTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTg                           >  1:251561/1‑136 (MQ=255)
                                                                                                               aGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCAt         <  1:408233/136‑1 (MQ=255)
                                                                                                               aGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCAt         <  1:189200/136‑1 (MQ=255)
                                                                                                                gCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCATc        <  1:206295/136‑1 (MQ=255)
                                                                                                                    ccTGAGTGATCAGGCCGGCGCACCGTTGGAGGTAATCATGAGCCACTATCTACGCGTGTTGCTGTTCATGGGCGCATTGCTGCACCAGGCGGTGGCGCTGGAGCTGGCCGGTGGTCTGGCCCGATCTGCATCGGcc    <  2:336763/136‑1 (MQ=255)
                                                                                                                    ccTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCATCGGcc    <  2:336762/136‑1 (MQ=255)
                                                                                                                      tGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGCTCAGCATCGGCCAg  >  2:27757/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                              
TGCAGGCGTTCGCACGCTTGATAAAAACCAAGCGCCAAGGGCGGCATGAACTGGATGCTGGCAGCAGCAGTCAATCACTGCCCGCCGAAAGCCCCCTCTCGATGCAACCCGAGCGACCTGAGTGATCAGGCCGGCCCACCGTTGGAGGTAATCATGAGCCACTATCAACGCCTGTTGCTGATCATGCGCCCAATGCTGCACCAGGCGGTCGCCCTGGAGCAGGCCGGTGCACTGGCCCGATCAGCATCGGCCAG  >  NC_002947/3674765‑3675018

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: