Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,608,512 G→A R171C (CGC→TGC)  yneE ← conserved inner membrane protein YneE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,5120GA95.7% 68.6 / NA 24R171C (CGC→TGC) yneEconserved inner membrane protein YneE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCG  >  NC_000913/1608382‑1608638
                                                                                                                                  |                                                                                                                              
aaaaTCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACCGGATATCCGAAAGCGGTCCATCGCGGCACt                                                                                                                              >  1:176543/1‑133 (MQ=255)
   aTCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGaa                                                                                                                           <  2:65151/133‑1 (MQ=255)
         gtgtAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCa                                                                                                                     >  1:264332/1‑133 (MQ=255)
            tAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAAcc                                                                                                                  >  1:455859/1‑133 (MQ=255)
             aGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGATAGCGGTCTATTGCGGCGCTGAACCGCCAACCa                                                                                                                 >  2:145522/1‑133 (MQ=255)
                  aaGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATtctc                                                                                                            >  1:119359/1‑133 (MQ=255)
                            gCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTaa                                                                                                  <  2:455859/133‑1 (MQ=255)
                            gCGTATAGCCAATGCCCTCGCACCCTGCCATGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTaa                                                                                                  <  1:145522/133‑1 (MQ=255)
                                 taGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGa                                                                                             >  2:471942/1‑133 (MQ=255)
                                   ggCAATGCTCTCGCACCCTGCCCGGACCGCAGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATa                                                                                           <  1:142896/133‑1 (MQ=255)
                                      aaTGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACgg                                                                                        >  2:378887/1‑133 (MQ=255)
                                       aTGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGt                                                                                       <  1:378887/133‑1 (MQ=255)
                                          cgcTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTAAACCGCCAACCAATCTCCCATTATTGACAAAATACGGTTAg                                                                                    >  2:206122/1‑133 (MQ=255)
                                                                cTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATAATTAACAAGATAAGGTGAGCTGGAGCGTTCGAAGCCAGTac                                                              >  2:144322/1‑133 (MQ=255)
                                                                 tGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTaca                                                             <  2:446509/133‑1 (MQ=255)
                                                                  gAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCTCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTacac                                                            <  1:471942/133‑1 (MQ=255)
                                                                      tatCATTACGCCGATTGTTGAGGCTAATAAGCAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTg                                                        <  1:206122/133‑1 (MQ=255)
                                                                             aaGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCt                                                 >  2:255308/1‑133 (MQ=255)
                                                                                   aTCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTc                                           >  2:488964/1‑133 (MQ=255)
                                                                                                    aCAGGATATCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAg                          <  2:257149/133‑1 (MQ=255)
                                                                                                         ataaCTGAAACCTGTGCATTGCGGGGCTCAACCGCCAATCAGTCTCCCATTATTAACAAGATCCGGTTAGCTGGAGAGTTCGAACCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACtt                     <  2:176543/133‑1 (MQ=255)
                                                                                                           atCTGAAATCTGTCCATTCCGGTACTGAACTGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTcc                   <  1:159827/133‑1 (MQ=255)
                                                                                                           atCTGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTcc                   <  1:255308/133‑1 (MQ=255)
                                                                                                              tGAAAGCTGTCCATTGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGcc                >  1:373360/1‑133 (MQ=255)
                                                                                                                            tGCGGCACTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCg  <  2:373360/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                              
AAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCG  >  NC_000913/1608382‑1608638

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