Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,283,291 C→T 25.2% C326Y (TGC→TAC)  PP_2885 ← transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,283,2910CT25.2% 15.3 / 5.0 16C326Y (TGC→TAC) PP_2885transporter
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base T (2/2);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

GCCCGATCACCGCGCCCATGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCTGGGCCAGCATCAGTGC  >  NC_002947/3283151‑3283427
                                                                                                                                            |                                                                                                                                        
gCCCGATCACCGCGCCCATGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGTAAGCCcaca                                                                                                                                 >  1:285387/1‑150 (MQ=255)
gCCCGATCACCGCGCCCATGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGTAAGCCcaca                                                                                                                                 >  1:285419/1‑150 (MQ=255)
               ccATGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCCTTCGCCGCCAATCCCGCAAACGCCGCGCTGACTAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCgtgg                                                                                                                  <  1:492864/150‑1 (MQ=255)
                 aTGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCgtggtg                                                                                                                >  2:3558/1‑150 (MQ=255)
                 aTGCCTGGGGCGGCGCCAATGGAACCGGCCTTGGCTCCTTGTTCGGAGGCCAGCATGGCATTGGCCGCCAGTGCCATAAACTTCGTTAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCgtggtg                                                                                                                <  2:250198/150‑1 (MQ=255)
                     cTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCCAAAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCa                                                                                                            >  1:149008/1‑150 (MQ=255)
                                 cAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCCCAgg                                                                                                >  2:48752/1‑150 (MQ=255)
                                       aaCCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATCTTCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGc                                                                                          <  1:148738/150‑1 (MQ=255)
                                                       ttGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGt                                                                          <  2:270473/150‑1 (MQ=255)
                                                         gTTCGGGGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTgg                                                                        <  2:521253/150‑1 (MQ=255)
                                                                                gccgccAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCGTCatc                                                 >  2:409259/1‑148 (MQ=255)
                                                                                                     gtcaaCACTAAGCCCATGCCGTCCGCCGCAACAAAATAGCAAGCCCACATCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGgcgc                            <  1:115889/147‑1 (MQ=255)
                                                                                                               aGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCt                  >  1:15962/1‑150 (MQ=255)
                                                                                                                 cccATGCCGGCCGCCGCAACAAAATAGTAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCATAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCTgg                <  2:285387/150‑1 (MQ=255)
                                                                                                                 cccATGCCGGCCGCCGCAACAAAATAGTAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCATAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCTgg                <  2:285419/150‑1 (MQ=255)
                                                                                                                               cgcAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCTGGGCCAGCATCAGTGc  <  1:146695/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
GCCCGATCACCGCGCCCATGCCTTGGGCGGCGCCAATGGAACCGGCCGTGGCGCCTTGTTCGGAGGCCTGCATGGCATTGGCCGCCAGTGCCGAAAACGCCGGGAACACGAAGCCCATGCCGGCCGCCGCAACAAAATAGCAAGCCCACAGCCACGGCGCCGTGGTGGCCAGCGAGCCACAGGCAAAGCCCAGTGCCGATACCGTGGCACCCACCCGGATCATCTTCAGCGGTGGCCACTCCAGTTGGCGCAGTATCACCTGGGCCAGCATCAGTGC  >  NC_002947/3283151‑3283427

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: