Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,072,962 G→C 11.8% S450R (AGC→AGG) ‡ PP_5327 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,072,9620GC11.8% 89.4 / 5.8 34S450R (AGC→AGG) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (16/14);  new base C (2/2);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGG  >  NC_002947/6072815‑6073106
                                                                                                                                                   |                                                                                                                                                
cATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCt                                                                                                                                                >  2:1059164/1‑150 (MQ=255)
  tGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTgg                                                                                                                                              >  2:609416/1‑150 (MQ=255)
      cttctGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCAt                                                                                                                                          >  2:768628/1‑150 (MQ=255)
              ccAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGa                                                                                                                                  <  1:1044137/150‑1 (MQ=255)
                  gTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGACAGTTCTGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGAt                                                                                                                              >  1:10912/1‑150 (MQ=255)
                                        tccaCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTg                                                                                                        >  1:415392/1‑150 (MQ=255)
                                             cGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTgg                                                                                                                                              >  2:5460/1‑107 (MQ=255)
                                             cGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTgg                                                                                                                                              <  1:5460/107‑1 (MQ=255)
                                                            cAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                 >  1:94195/1‑137 (MQ=255)
                                                            cAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                 <  2:94195/137‑1 (MQ=255)
                                                                     cacaCCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCtt                                                                                                         >  2:881784/1‑120 (MQ=255)
                                                                     cacaCCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCtt                                                                                                         <  1:881784/120‑1 (MQ=255)
                                                                         ccTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcacgcacg                                                                       <  1:609416/150‑1 (MQ=255)
                                                                          cTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACgg                                                                      >  1:5194/1‑150 (MQ=255)
                                                                          cTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCAGCCCGGCCATga                                                                                                                                        <  2:890541/84‑2 (MQ=255)
                                                                          cTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCAGCCCGGCCATga                                                                                                                                        >  1:890541/1‑83 (MQ=255)
                                                                             ccAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTgc                                                                   <  2:384310/150‑1 (MQ=255)
                                                                                      ccggcgCCACTTCGCCGGTGGCACGGGCTACAGCGAGGATCATGACGGTCATCATGGCCGGGCTGGCCATCGGCAGGACGATTTTCCCGAGCGTATCGGCCTTGGTTGCGCCCCAAGCCAGCTAGCCTTCTCGCACGGTGTGCGGAATTc                                                          <  2:100409/149‑1 (MQ=255)
                                                                                              aCTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGGCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAgg                                                  >  1:651684/1‑150 (MQ=255)
                                                                                                ttCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                 >  2:1072034/1‑101 (MQ=255)
                                                                                                ttCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                 <  1:1072034/101‑1 (MQ=255)
                                                                                                       gcggcACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCg                                         >  1:222675/1‑150 (MQ=255)
                                                                                                            aCGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGc                                    >  2:265576/1‑150 (MQ=255)
                                                                                                              aggCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCa                                  >  1:876900/2‑150 (MQ=255)
                                                                                                                  cacaGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGAt                              <  1:768628/150‑1 (MQ=255)
                                                                                                                   acaGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATc                             <  1:172728/150‑1 (MQ=255)
                                                                                                                                     ctgcccatGGCCAGCCCGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGc           >  1:832650/6‑150 (MQ=255)
                                                                                                                                      tgcccatGGCCAGCCCGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAg                                                   <  1:977283/105‑1 (MQ=255)
                                                                                                                                         tcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATg       >  2:1032523/1‑150 (MQ=255)
                                                                                                                                         tcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATg       <  1:1059164/150‑1 (MQ=255)
                                                                                                                                            tGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCa    >  1:547189/1‑150 (MQ=255)
                                                                                                                                            tGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCCGGACCGCCAGCAATGCCa    <  2:651684/150‑1 (MQ=255)
                                                                                                                                             ggCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAg   <  2:876900/150‑1 (MQ=255)
                                                                                                                                              gCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAgg  <  2:5194/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGG  >  NC_002947/6072815‑6073106

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