Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,453,565 C→A 23.5% P64P (CCC→CCA PP_1272 → multidrug MFS transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,453,5650CA23.5% 29.2 / 4.5 17P64P (CCC→CCAPP_1272multidrug MFS transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (8/5);  new base A (2/2);  total (10/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.64e-01

CCAACAACCGCAAAACCCTTTTCATCGGCTCGGTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTGCGCAGGCGCAGTTGCTGGTAGCCC  >  NC_002947/1453420‑1453701
                                                                                                                                                 |                                                                                                                                        
ccAACAACCGCAAAACCCTTTTCATCGGCTCGGTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCAAAAg                                                                                                                                      >  2:330246/1‑150 (MQ=255)
ccAACAACCGCAAAACCCTTTTCATCGGCTCGGTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCAAAAg                                                                                                                                      >  2:330242/1‑150 (MQ=255)
                ccTTTTCATCGGCTCGGTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATcaa                                                                                                                      >  1:136249/1‑150 (MQ=255)
                            cTCGGTGCTGGCCGTGGCCCTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGCCGCCTATGTGCCCGCCGACTACACCGTGGTCGCGCCCAACGTCGCCGGTTTCATCAAGCAAGTGCtggt                                                                                                          <  1:210569/150‑1 (MQ=255)
                                gTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGaa                                                                                                      <  1:461309/150‑1 (MQ=255)
                                   cTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCAAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGAc                                                                                                   <  1:330242/150‑1 (MQ=255)
                                   cTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCAAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGAc                                                                                                   <  1:330246/150‑1 (MQ=255)
                                       ccgtggccgtCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATc                                                                                               >  1:161476/1‑150 (MQ=255)
                                                                    cccTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCGCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCAAGTTGCTGGca                                                                  >  2:88339/1‑149 (MQ=255)
                                                                      cTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGAAAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGAc                                                                >  2:119136/1‑150 (MQ=255)
                                                                                                     aCCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTg                                 >  1:257961/1‑150 (MQ=255)
                                                                                                        aaTGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGAt                              <  2:161476/150‑1 (MQ=255)
                                                                                                        aaTGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTGCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGAt                              >  1:257784/1‑150 (MQ=255)
                                                                                                                   tgtgACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCTTCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTgcgcaggc                   <  1:437361/150‑1 (MQ=255)
                                                                                                                           ccgACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTGCGCAGGCGCAGTTGc           >  1:333150/1‑150 (MQ=255)
                                                                                                                              aCTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTGCGCAGGCGACGTTGATgg        <  2:136249/150‑1 (MQ=255)
                                                                                                                                    ccGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGCCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTGCGCAGGCGCAGTTGCTGGTAGccc  >  1:139778/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
CCAACAACCGCAAAACCCTTTTCATCGGCTCGGTGCTGGCCGTGGCCGTCCTGGCCGGCATCGTCGGCCCCTGGATGTTTGGCAGCGATCACCGCCAGAGCACCAATGACGCCTATGTGACCGCCGACTACACCGTGGTCGCGCCCAAAGTCGCCGGTTTCATCAAGCAAGTGCTGGTGGAAGACAATCAGCAGGTCACAGCTGGCCAGTTGCTGGCGACCATCGACGCGCGTGACTACCAGGCCGCGCTGGATGCTGCGCAGGCGCAGTTGCTGGTAGCCC  >  NC_002947/1453420‑1453701

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