Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,724,285 G→A R53Q (CGG→CAG)  yfiM → DUF2279 domain‑containing protein YfiM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,724,2850GA92.6% 71.8 / ‑1.1 28R53Q (CGG→CAG) yfiMDUF2279 domain‑containing protein YfiM
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (9/16);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

ATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACA  >  NC_000913/2724047‑2724523
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
aTCCGCATCGACCGATTAATTATCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCTGCCGGAAATGAATATTCACAGCATCAGGTGATGAGCCAGGATCGCAGTg                                                                                                                                                                                                                                     <  1:312785/250‑1 (MQ=255)
        cGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGAAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTACGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTg                                                                                                                                                                                                                            <  1:118991/251‑1 (MQ=255)
        cGACCGAGTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTACGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTg                                                                                                                                                                                                                            <  1:117319/251‑1 (MQ=255)
                taattaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCTTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGATAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTACGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGAt                                                                                                                                                                                                                     <  1:301411/250‑1 (MQ=255)
                  attaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCCTCAGGGGATGAGCCCGGATCGCAGTGCCATGTTTGGATTGATGtt                                                                                                                                                                                                                  >  1:288536/1‑251 (MQ=255)
                     aGCCGCATCCTGTAATCACAACCCCGGCCTGTACGGGGTTTGTTTTTTGGAGGACACGTTTTGCGAATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCCTATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTAACGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCGGTGCAATGTTTGGTTTGagtttctc                                                                                                                                                                                                               >  2:379682/1‑244 (MQ=255)
                                ggaaTCACAACCCCTTCCTGTACGGGGTTTGTTTGTGGGAGGCCACGTTTTGCGTATTCTTGTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTAGCGCCTCGGCGATGCTGTGCGCCGCCGGAAATTAATATTCACAGCAGCAGGGGATGGGTCAGGATCGCAGTGCCATGTGTGGATTGATGTTCGCTGGCAGttt                                                                                                                                                                                                      <  2:173778/247‑1 (MQ=255)
                                   aTCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTggggg                                                                                                                                                                                                 <  1:12737/251‑1 (MQ=255)
                                       caACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTc                                                                                                                                                                                             <  2:279215/251‑1 (MQ=255)
                                                ccTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGc                                                                                                                                                                                     <  1:291401/250‑1 (MQ=255)
                                                  tGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTc                                                                                                                                                                                                           <  1:187372/226‑1 (MQ=255)
                                                  tGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTc                                                                                                                                                                                                           >  2:187372/1‑226 (MQ=255)
                                                       ggggTTAGTTTTTTGGAGGCCACGTTTGGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCCGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGCGCTTTGGGAt                                                                                                                                                                             >  2:269375/1‑251 (MQ=255)
                                                                 ttttGGAGGCCACGGTTTGCGTATTCTTTTTTGCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGCAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGcccc                                                                                                                                                                    <  1:69790/250‑1 (MQ=255)
                                                                          ccACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAACGAGCTTTGGGATAGCCGCCCCGCAGGGAGCg                                                                                                                                                          >  2:175931/1‑251 (MQ=255)
                                                                                                    gTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTg                                                                                                                                 <  2:42913/250‑1 (MQ=255)
                                                                                                             tGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGt                                                                                                                                                                                                        <  1:292958/170‑1 (MQ=255)
                                                                                                             tGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGt                                                                                                                                                                                                        >  2:292958/1‑170 (MQ=255)
                                                                                                                              aGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCt                                                                                                       >  2:89608/1‑250 (MQ=255)
                                                                                                                                          aaCGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGc                                                                                           <  2:108729/250‑1 (MQ=255)
                                                                                                                                                                               ttATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCttt                                                      <  1:89608/250‑1 (MQ=255)
                                                                                                                                                                                            gATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATc                                         >  1:324496/1‑250 (MQ=255)
                                                                                                                                                                                                 tttCCGCCGCCGGAAATGAATATTCACACCATCAGGGGATGAGACAGGATCGCAGTGCAATGTTTGGATTGATGATCACTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGACGCCCCGAAGGGAGCGGCTGAAGCTGGAAGGATTTGGCCTGGGAAGTCGACGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACACGTCACTAAAAACGCATCCCCTTCACCTTGCGATGTGGCATCAAAGa                                    >  1:288787/1‑250 (MQ=255)
                                                                                                                                                                                                   tcggccgccgGAAATGAATATTCACAGAATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGaaa                                  <  2:117995/247‑1 (MQ=255)
                                                                                                                                                                                                       ccgccgGAAATGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGa                             <  1:187720/251‑1 (MQ=255)
                                                                                                                                                                                                                aTGAATATTCACAGCATCAGGGGATGAGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacc                    >  1:38815/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cAGGGGATGGGCCAGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTca    <  2:324496/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aggggaggggCCAGGATCTCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTcaca  <  2:248079/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
ATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCGCCGCCGGAAATGAATATTCACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACA  >  NC_000913/2724047‑2724523

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