Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,106,743 A→C 46.4% L186R (CTG→CGG)  PP_1878 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,106,7430AC46.4% ‑6.4 / 3.2 13L186R (CTG→CGG) PP_1878hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (2/4);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  NC_002947/2106617‑2106871
                                                                                                                              |                                                                                                                                
ggggCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGccc                                                                                                                         >  1:221046/1‑135 (MQ=255)
      tggCGGCCCAGTACGTGGGATTGGGCCGGGGCCAACAGTTTCTCAAGAAACGTCAGCTGCTCTACATACCCCTTCGGCGCCGGCGGCTGGCGAAAATACTCCACGCCGCGCTCGACGGCCCGGCGCgccgggccgg                                                                                                                   <  2:164988/134‑1 (MQ=255)
                 tACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGGCCGGCggggcggtt                                                                                                        >  2:233763/1‑130 (MQ=255)
                         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGCCGTTTTTCAgggg                                                                                                >  1:453898/1‑135 (MQ=255)
                         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCCGGGTGGGCGTTTTTAAgggg                                                                                                >  1:106737/1‑135 (MQ=255)
                         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCCCCTGGGCGCGCTGGGGGTTTTGCTGGGc                                                                                                >  1:34842/1‑136 (MQ=255)
                                    cccAACATTTTCTCATGAGACGTCTGCTGGTCTTCAAACCGCTTCGGCTCGGGCTGGTGGCTCAAATCCGGCCCGCCGCGCTCGTCGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGa                                                                                     <  2:427752/135‑1 (MQ=255)
                                           ttttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                              <  1:161682/135‑1 (MQ=255)
                                                caaaGATCCGTCAGCTGCTCTTCATCCACCTGCGGCGCGGGCGGCGTGCTCAAATACGCCCCGCCGCGCTCGCAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                         <  1:427819/134‑1 (MQ=255)
                                                                gctCTACATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                         <  2:336477/136‑1 (MQ=255)
                                                                                                              gccgcGCTCGAAGGCCAGGCGCGCCTGTCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTTGCCTCGAcc           >  1:254234/1‑135 (MQ=255)
                                                                                                                     tCGAAGGCCAGGCGCGCCTTGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCa    >  2:386647/1‑136 (MQ=255)
                                                                                                                       taaGGCCCGGCGCGGCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  2:453898/135‑1 (MQ=255)
                                                                                                                              |                                                                                                                                
GGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  NC_002947/2106617‑2106871

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