Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,762,606 C→G 25.0% L38L (CTG→CTC PP_5055 ← Sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,762,6060CG25.0% 25.4 / 7.0 16L38L (CTG→CTCPP_5055Sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base G (2/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCGCGCCGACGATATGACCGGCGGAGTATTCCTTGGCCGGGCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTCGCGAATTGAATCAGGTGAGCAACCATCAGCGGTGT  >  NC_002947/5762506‑5762728
                                                                                                    |                                                                                                                          
ccGCGCCGACGATATGACCGGCGGAGTATTCCTTGGCCGGGCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGag                                                                                         <  1:334290/136‑1 (MQ=255)
 cgcgCCGACGATATGACCGGCGGAGTATTCCTTGGCCGGGCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGagc                                                                                        <  1:109035/136‑1 (MQ=255)
 cgcgCCGACGATATGACCGGCGGAGTATTCCTTGGCCGGGCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGagc                                                                                        <  1:183655/136‑1 (MQ=255)
                                       ggCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGaa                                                                 <  1:170950/121‑1 (MQ=255)
                                       ggCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGaa                                                                 >  2:170950/1‑121 (MQ=255)
                                          gAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTa                                               >  2:498272/1‑136 (MQ=255)
                                                     aCCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTcg                                    >  2:459179/1‑136 (MQ=255)
                                                          ggCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCg                                                            <  1:383708/107‑1 (MQ=255)
                                                          ggCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCg                                                            >  2:383708/1‑107 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTGAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAAcc                                                       <  1:354450/93‑1 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTGAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAAcc                                                       <  1:357170/93‑1 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTGAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAAcc                                                       >  2:354450/1‑93 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTGAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAAcc                                                       >  2:357170/1‑93 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTCGCGAATTGAATCAGGTGAGCAACCa            >  1:109447/1‑136 (MQ=255)
                                                                             aGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTCGCGAATTGAATCAGGTGAGCAACCa            >  2:237019/1‑136 (MQ=255)
                                                                                       gTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTCGCGAATTGAATCAGGTGAGCAACCATCAGCGgtgt  >  1:102151/1‑136 (MQ=255)
                                                                                                    |                                                                                                                          
CCGCGCCGACGATATGACCGGCGGAGTATTCCTTGGCCGGGCGAATGTCGATCACCAGGGCCTTCTCGGCATTGACCAGGGCTGTCAGTTGGCCGTTGGTCAGGCTTTGGCCACCACGACGGATTTCATTGATGAGCAGCAGGACCAGCAGAATGAGGAAGATCGCAACCAGGATGTAGTGATCTGTCGCGAATTGAATCAGGTGAGCAACCATCAGCGGTGT  >  NC_002947/5762506‑5762728

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