Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,510,111 G→A G161D (GGC→GAC)  yfeO → putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,510,1110GA92.1% 96.7 / 1.4 38G161D (GGC→GAC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (19/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

CCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCG  >  NC_000913/2509977‑2510230
                                                                                                                                      |                                                                                                                       
ccTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACt                                                                                                                     >  2:1321210/1‑139 (MQ=255)
                         aCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATAtttt                                                                                            >  1:1521721/1‑139 (MQ=255)
                             tCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCa                                                                                        >  1:838302/1‑139 (MQ=255)
                                           gCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCa                                                                          >  2:829887/1‑139 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTgcggcg                                                                                                         <  2:656984/102‑1 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTgcggcg                                                                                                         >  1:656984/1‑102 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCaaa                                                                                      >  2:429984/1‑121 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCaaa                                                                                      <  1:429984/121‑1 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTg                                                                    >  2:1279983/1‑139 (MQ=255)
                                                 gCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTg                                                                    >  1:933102/1‑139 (MQ=255)
                                                      tGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGtt                                                               <  2:1707242/139‑1 (MQ=255)
                                                       ggCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTc                                                              <  1:1279983/139‑1 (MQ=255)
                                                       ggCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTc                                                              <  2:933102/139‑1 (MQ=255)
                                                                        gcgcgTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTc                                                              >  2:408804/1‑122 (MQ=255)
                                                                        gcgcgTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTc                                                              <  1:408804/122‑1 (MQ=255)
                                                                        gcgcgTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGtctc                                             >  1:1224729/1‑139 (MQ=255)
                                                                           cgTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAAt                                                                              <  1:1676751/103‑1 (MQ=255)
                                                                           cgTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAAt                                                                              >  2:1676751/1‑103 (MQ=255)
                                                                            gTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTg                                         <  1:1321210/139‑1 (MQ=255)
                                                                               aaCCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTgcgc                                      <  1:829887/139‑1 (MQ=255)
                                                                                      tGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATgg                                                                            <  1:715835/94‑1 (MQ=255)
                                                                                      tGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATgg                                                                            >  2:715835/1‑94 (MQ=255)
                                                                                                   tttAGCCTCTGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCt                  >  1:1279695/1‑139 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACtt         >  1:273387/1‑139 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCa                                                                                        >  2:375/1‑60 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCa                                                                                        <  1:375/60‑1 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTgc                                        <  2:1316523/108‑1 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTgc                                        >  1:1316523/1‑108 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATggcgg                          <  1:998906/122‑1 (MQ=255)
                                                                                                            tGCCGGAACCATCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATggcgg                          >  2:998906/1‑122 (MQ=255)
                                                                                                                      aTCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTgcg                                       >  2:261058/1‑99 (MQ=255)
                                                                                                                      aTCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTgcg                                       <  1:261058/99‑1 (MQ=255)
                                                                                                                      aTCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTg              <  1:915002/124‑1 (MQ=255)
                                                                                                                      aTCGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTg              >  2:915002/1‑124 (MQ=255)
                                                                                                                       tCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCg  <  2:279836/135‑1 (MQ=255)
                                                                                                                       tCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCg  >  1:279836/1‑135 (MQ=255)
                                                                                                                        cGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGc                                                                                         <  2:1317447/47‑1 (MQ=255)
                                                                                                                        cGGCGCACTGTTTGACACTCCTGTTGCGGCGGCGTTGATATTTTCGc                                                                                         >  1:1317447/1‑47 (MQ=255)
                                                                                                                                      |                                                                                                                       
CCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCGTTGATATTTTCGCAAACCTTAAATGGCAGTAGTGAAGTTCCGCTATGGGATCGTCTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCG  >  NC_000913/2509977‑2510230

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