Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,447,744 G→A 8.9% A232V (GCA→GTA)  prmB ← N5‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,7440GA 8.9% 119.9 / 2.8 47A232V (GCA→GTA) prmBN5‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/24)
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.51e-01

CCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGAC  >  NC_000913/2447599‑2447886
                                                                                                                                                 |                                                                                                                                              
ccagGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGACAGTGCCAGATGCCa                                                                                                                                             <  2:575449/149‑1 (MQ=255)
ccagGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAAGGCGACGCGGCAGTTTCAGGCCGTCAGGGCCAGCTGCCa                                                                                                                                             >  2:510052/1‑149 (MQ=255)
           gacactCCGGTTGTTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACCCGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTATCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTcc                                                                                                                                  <  2:548721/144‑1 (MQ=255)
           gAACATCGGGGTATTGTTCCATAAGATTGCCCATGCTGTTGCCGACTTCACAAATCATTACGCTATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTcc                                                                                                                                  <  1:10171/149‑1 (MQ=255)
            aaCATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATACCAGGCCCAGTTCCg                                                                                                                                 >  1:820807/1‑149 (MQ=255)
            aaCATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATACCAGGCCCAGTTCCg                                                                                                                                 >  1:820777/1‑149 (MQ=255)
                           ttCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCGCGTGGCGGTAtt                                                                                                                  >  2:93442/1‑149 (MQ=255)
                                    aTGTACCATTCTGTTGCCGACTTCACAAATCAACACGCCATCACCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCa                                                                                                         <  2:580291/149‑1 (MQ=255)
                                        aCCATGCTGTTGCCTACTTCACAAATCACCACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTc                                                                                                     <  2:372579/149‑1 (MQ=255)
                                        aCCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTc                                                                                                     <  1:348577/149‑1 (MQ=255)
                                           aTGCCGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGa                                                                                                  <  1:190814/149‑1 (MQ=255)
                                               tGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCCGTTTCAGGCCGTCAGTGACAGCTGCCAGGCCCAGTTCCGGCCCGTGGCGGCATTCGTTTGGCAggacggaccaa                                                                                              >  2:692885/1‑146 (MQ=255)
                                                 ttGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCCGTGCCAGATGCCAGGCCCAGTTCCGCCTCGTGGCGGTATCCGCTTGCCCGGTCGGACAAATc                                                                                            >  1:101996/1‑149 (MQ=255)
                                                    ccGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTc                                                                                         >  2:612701/1‑149 (MQ=255)
                                                     cGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCGGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTcc                                                                                        <  1:631174/149‑1 (MQ=255)
                                                      gACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCg                                                                                       <  1:380983/149‑1 (MQ=255)
                                                            acaAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa                                                                                 >  1:861841/1‑149 (MQ=255)
                                                            acaAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGGCAGTTTCCGGCCGTCAGTGCCAGATGCCAGGCCCAGTGCCGGCTCGTGGCGGTATTCGTTTGGCGGGTCGGACGTAACGTCCGCCTCGa                                                                                 >  2:397382/1‑149 (MQ=255)
                                                             caAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGGATTCGTTTGGCAGGTCGGACAAATCTTCCGCCTCGAc                                                                                >  1:159299/1‑149 (MQ=255)
                                                                     cccGCCTTCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATAcggcg                                                                        <  2:835377/147‑1 (MQ=255)
                                                                       cGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGGCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGCATTCGTTTGGCAGGTCGGACATATCTTCCGCGTCGGCATACGGCggg                                                                      >  2:675146/1‑149 (MQ=255)
                                                                            tcatcaGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGt                                                                 <  2:452223/149‑1 (MQ=255)
                                                                            tcatcaGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTAGTCGTTTGGCAGGTCGGACATATCCTCCGCAGCGGCGTACGTCGGGTTAgg                                                                 >  2:167866/1‑148 (MQ=255)
                                                                              atcaGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATGCGGTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGa                                                               >  2:884840/1‑149 (MQ=255)
                                                                                caGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATGTTCCGCATCGACATGCGGCGGGTTAGTGACa                                                             >  2:293614/1‑149 (MQ=255)
                                                                                 aGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGCTGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACaa                                                            >  1:894288/1‑149 (MQ=255)
                                                                                 aGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACaa                                                            <  1:857774/149‑1 (MQ=255)
                                                                                  gCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATACCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAAt                                                           <  2:820777/149‑1 (MQ=255)
                                                                                  gCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATACCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAAt                                                           <  2:820807/149‑1 (MQ=255)
                                                                                    catgTACTCTGCCGCGTTACCGATATGGCGCCGCGTCAGTTTCAGGCCGTCAGTGCCGGATGCCAGGCCCAGTTCCGGCTCGCTGCGGTATTCGTTTGGCAGGTCGTACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCa                                                         <  2:227706/146‑1 (MQ=255)
                                                                                    aaGGCAATCGCCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGGGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCTACATACGGCGGGTTAGTGACAATCa                                                         <  1:460858/149‑1 (MQ=255)
                                                                                       gTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGt                                                      >  1:138095/1‑149 (MQ=255)
                                                                                           tCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTa                                                  <  1:383316/149‑1 (MQ=255)
                                                                                           tCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTa                                                  <  1:578833/149‑1 (MQ=255)
                                                                                                cccattACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGCCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCa                                             <  1:614286/145‑1 (MQ=255)
                                                                                                              tGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGCTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTcg                               >  2:747523/1‑149 (MQ=255)
                                                                                                               gccacccgTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTcgc                              <  2:865014/143‑1 (MQ=255)
                                                                                                                   tgcgTCCGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCAGCGGCATACGGCGGGTTAGTGACAATCCGGTCGTACTGCACTTTCGGCAAGTCGCGGaa                          <  2:118521/148‑1 (MQ=255)
                                                                                                                      gTCAGTTTCAGGCCGTCAGTGCCAGATGCAAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATAACTTCCGCATCTACATACGGCGGGGTAGTGACAATCAGGTAGGACTGCCCATACGGCAAGGCGCGGAACac                       >  2:69444/1‑148 (MQ=255)
                                                                                                                         aGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCTGACATATCTTCCGCATCGACATACGGCAGGTTAGTGACAATCAGGTCGTACTGCACTTTCCGCAAGTCGCGGAACAGATc                    >  1:674992/1‑149 (MQ=255)
                                                                                                                         aGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATc                    >  1:487778/1‑149 (MQ=255)
                                                                                                                           tttCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTAGTGCACTTTCGGACGGTCGCGGAACAGATCgg                  >  2:179420/1‑149 (MQ=255)
                                                                                                                              ccggCCGTCAGGTCCAGATGTCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGGGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAAc               <  2:809701/147‑1 (MQ=255)
                                                                                                                                   cGTCAGTGCCAGATCCCGGGGCCAGTTCCGGCGCGTGGCGGTATTCGTTTGGCAGGTCGCACATATCTTCCCCATCGACAGACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATc          <  2:328323/149‑1 (MQ=255)
                                                                                                                                    gTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCg         <  1:124685/149‑1 (MQ=255)
                                                                                                                                         tGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCACGTCGCGGCACAGATCGGAACGAATCGGAATg    >  2:842446/1‑149 (MQ=255)
                                                                                                                                           ccAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGGCATGTCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGAc  >  2:39567/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGACATACGGCGGGTTAGTGACAATCAGGTCGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGAC  >  NC_000913/2447599‑2447886

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