Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,401,467 G→A 9.1% A265A (GCG→GCA)  hflX → GTPase, stimulated by 50S subunit binding; Mn(2+) homeostasis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,401,4670GA 9.1% 106.9 / 3.6 44A265A (GCG→GCA) hflXGTPase, stimulated by 50S subunit binding; Mn(2+) homeostasis regulator
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base A (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

TCACCGAAGCGCGGGTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACC  >  NC_000913_3_bme_pgi/4401322‑4401608
                                                                                                                                                 |                                                                                                                                             
tCACCGAAGCGCGGGTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCACCACCTGCCGCACGATCTGGTGGCGGCa                                                                                                                                            >  2:421051/1‑149 (MQ=255)
      aaGCGCGGGTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTaaa                                                                                                                                      <  1:435769/149‑1 (MQ=255)
            gggTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACg                                                                                                                                <  1:514716/149‑1 (MQ=255)
             tgTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGt                                                                                                                               >  1:482173/2‑149 (MQ=255)
                  aggcgGCAGACCAGTTTTTTCCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGCGACCCCGTGCCTACCGGTGACGGTGGTTTTTTTCGCAACCTGCTTCACGATCTGGTGCCGGGATTTTAAGCCACGTTACaa                                                                                                                          <  2:639916/147‑1 (MQ=255)
                    gcgGTAGACTAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGAAGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAaga                                                                                                                        <  1:120967/149‑1 (MQ=255)
                       gCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCCCGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGAc                                                                                                                     >  1:870630/1‑149 (MQ=255)
                          gACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTCTTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGCGAcgcg                                                                                                                  >  1:233147/1‑149 (MQ=255)
                              ggtttttTGCCACCCTCGACCCGACGTGGCGGTGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCCCCGCCGGCCGGACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCaa                                                                                                              <  2:370664/148‑1 (MQ=255)
                              aGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCAGCATTTAAAGCCACGTTACAAGAGACGCGGCaa                                                                                                              >  1:39250/1‑149 (MQ=255)
                              aGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCAGCATTTAAAGCCACGTTACAAGAGACGCGGCaa                                                                                                              >  1:40168/1‑149 (MQ=255)
                                ttgttTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACGTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGc                                                                                                            <  1:849220/149‑1 (MQ=255)
                                       ccACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTa                                                                                                     >  1:719210/1‑149 (MQ=255)
                                           ccTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGACCCCCCGCAAGCCCCATTActgc                                                                                                 >  2:578482/1‑149 (MQ=255)
                                            cTCGACCCGACGTTGCGGCGTAGTGACGTTGCAGATGTCGGTGCAACCGTACTTGCAGATACCGTGGGTGGTATTCGCCGCCGGGCGCACGATCCGATCGTGGCATTTGAAGTCAACGTGCAAACGACACGGCACGCCAGCTTAcccgt                                                                                                >  2:677466/1‑145 (MQ=255)
                                             tCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCCTTTAAAGCCACACTACAAGACACGCCGCACGCCACATTActgctg                                                                                               >  2:476975/1‑149 (MQ=255)
                                                   caaccgtGCGGCGTATTAACGTTGCAGAGGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTc                                                                                         <  1:444431/143‑1 (MQ=255)
                                                      cGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCAtt                                                                                      <  1:679608/149‑1 (MQ=255)
                                                        tcgcggcgTATTGTCGTTGCAGATGTCGGGGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGa                                                                                    <  2:389212/147‑1 (MQ=255)
                                                                              aggTCGATGAAACCGTCCTTTCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGTTCTGGTGGCGGCATTCAAAGCCACGTTACAAGCGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTaaaa                                                              <  2:235364/147‑1 (MQ=255)
                                                                              aTGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCAGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACaa                                                              <  2:40168/149‑1 (MQ=255)
                                                                              aTGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCAGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACaa                                                              <  2:39250/149‑1 (MQ=255)
                                                                               tGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATATGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAg                                                             >  2:201297/1‑149 (MQ=255)
                                                                                 tCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGAAGCGTCGGATGTGCGTGTACAAGaa                                                           >  1:67240/1‑149 (MQ=255)
                                                                                   ggTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGaaaa                                                         <  1:519372/149‑1 (MQ=255)
                                                                                     tGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAATGCCACGTTACAAgaga                                                                                                                      <  1:784614/86‑1 (MQ=255)
                                                                                     tGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAgaga                                                                                                                      >  2:784614/1‑86 (MQ=255)
                                                                                          ccGTACTTGCAGATACCGTAGGGTTTATTCGCCACATGACAAAAGATCTGGTGGCGAAATCGAAAACCACTTTACAAGAAACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCaaa                                                  >  2:324848/1‑149 (MQ=255)
                                                                                             tACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCAGTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGGGCGTGTACAAGAAAACATCGAAGCg                                               >  2:468002/1‑149 (MQ=255)
                                                                                                     gATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATAc                                       >  1:210656/1‑149 (MQ=255)
                                                                                                           gTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAAAACGGTTCt                                 >  2:836485/1‑149 (MQ=255)
                                                                                                           gTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGGACAAGAAAACATCGAAGCGGTGAATACGGGTCt                                 >  1:40893/1‑149 (MQ=255)
                                                                                                            tAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCtt                                <  1:443093/149‑1 (MQ=255)
                                                                                                               ggTTTATTCGCCACCTGCCGCCTGATCTGGTGGCGGCATTTAAAGCCACCTTCCAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGCAAACATCGAAGCGGTGAATACGGTTCTTGaa                             <  2:617668/149‑1 (MQ=255)
                                                                                                               ggTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCATGTTTTTGTCGCGACGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTgga                             >  1:429701/1‑147 (MQ=255)
                                                                                                                gTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGCCGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGAGTACAAGAAAAAATCGAAGCGGTGAATACGGTTCTTCAag                            >  2:383197/1‑149 (MQ=255)
                                                                                                                    aTTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCTACGTTACTAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGAt                        <  1:102467/149‑1 (MQ=255)
                                                                                                                     ttCGCAACCTGCCGCAAGATCTGGTGGAGGCATTTAAAGCCACGTTACAAGAGAAGCGGCAAGCGACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGGCGCGGTGAATACGGTTATTGAagatact                       >  2:388614/1‑145 (MQ=255)
                                                                                                                        gCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGAc                    >  1:271654/1‑149 (MQ=255)
                                                                                                                         ccACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGCCAACGTAGAGAATGAAGCGGTGAATACGGTTCTTGAGAAGATCGACg                   >  1:615997/1‑149 (MQ=255)
                                                                                                                                 cGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCCCGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCAAAGCGCTGAATACGGTTCTTGAAGCGAACGACCCGCACaag           >  2:117078/1‑149 (MQ=255)
                                                                                                                                  gCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGCCGTGCGTGTACAAGAAAACATCCAAGCGATGAACACGGAACTTGAAAAGAACGACGCTCAGgaga          >  1:240764/1‑149 (MQ=255)
                                                                                                                                     cGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCAGATTACTGCTGCACGTGATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGTGGTGTGCGGTCACGAGATcc       >  1:866170/1‑149 (MQ=255)
                                                                                                                                          tggtggCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACTTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGATCTTGAAGAGATCGACGCTCACGAGAACCCAAcc  >  1:418149/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TCACCGAAGCGCGGGTCTACGCGGCAGACCAGTTGTTTGCCACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACC  >  NC_000913_3_bme_pgi/4401322‑4401608

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