Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,478,802 C→T *602* (TAG→TAA)  kefB ← potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,478,8020CT97.8% 143.5 / NA 46*602* (TAG→TAA) kefBpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/22);  minor base A (1/0);  total (24/22)
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

AACACATTCAACAATATCAATATTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTCTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCT  >  NC_000913/3478685‑3478935
                                                                                                                     |                                                                                                                                     
aaCACATTCAACAATATCAATATTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCg                                                                                                                  <  2:471399/139‑1 (MQ=255)
       tCAACAATATCAATATTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCt                                                                                                           >  1:94044/1‑139 (MQ=255)
            aatatcaatatTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgt                                                                                                      <  1:499484/139‑1 (MQ=255)
              tatcaatatTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgtcg                                                                                                    >  1:359731/1‑139 (MQ=255)
               atcaatatTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgtcgt                                                                                                   >  2:1162097/1‑139 (MQ=255)
               atcaatatTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgtcgt                                                                                                   <  1:729719/139‑1 (MQ=255)
                   atatTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCtt                                                                                               >  1:934876/1‑139 (MQ=255)
                    tattatTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCttg                                                                                              <  2:613115/139‑1 (MQ=255)
                         ttttCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCa                                                                                         >  1:1016888/1‑139 (MQ=255)
                              gcgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAAtt                                                                                                                             <  2:832045/98‑1 (MQ=255)
                              gcgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAAtt                                                                                                                             >  1:832045/1‑98 (MQ=255)
                              gcgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATccc                                                                                                                       >  1:184447/1‑104 (MQ=255)
                              gcgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATccc                                                                                                                       <  2:184447/104‑1 (MQ=255)
                              gcgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCa                                                                                    >  1:1482778/1‑139 (MQ=255)
                               cgcCACATCGCGGTTGAATCCTCCGCCTGACAGGCCGGGCATTTTTCGCCCGCAATAAACCGTTTACGGATTGCATTCTTTACCCTATACTCAACTTCATCCCAGCCGTCCAGCTGGCGTCGTTCTCGTTGCATTTCAc                                                                                   >  2:1101800/1‑139 (MQ=255)
                               cgcCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCAc                                                                                   <  1:1162097/139‑1 (MQ=255)
                                       cGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGaaa                                                                           <  2:934876/139‑1 (MQ=255)
                                            ttGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTc                                                                                                 <  1:534287/112‑1 (MQ=255)
                                            ttGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTc                                                                                                 >  2:534287/1‑112 (MQ=255)
                                               aaTCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGTCGTCCAGCt                                                                                                           <  2:5103/99‑1 (MQ=255)
                                               aaTCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGTCGTCCAGCt                                                                                                           >  1:5103/1‑99 (MQ=255)
                                                 tCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGa                                                                             <  1:1228053/127‑1 (MQ=255)
                                                 tCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGa                                                                             >  2:1228053/1‑127 (MQ=255)
                                                 tCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTACCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCa                                                                 >  2:823229/1‑139 (MQ=255)
                                                         cTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTACCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGt                                                         <  1:823229/139‑1 (MQ=255)
                                                             cAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCggg                                                     <  2:94044/139‑1 (MQ=255)
                                                                gCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCtt                                                  >  2:130100/1‑139 (MQ=255)
                                                                      cATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAAt                                                                         <  2:446543/110‑1 (MQ=255)
                                                                      cATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAAt                                                                         >  1:446543/1‑110 (MQ=255)
                                                                                   cAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgt                                                                                                      >  2:1186245/1‑68 (MQ=255)
                                                                                   cAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGcgt                                                                                                      <  1:1186245/68‑1 (MQ=255)
                                                                                   cAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCgg                                                      >  1:559577/1‑116 (MQ=255)
                                                                                   cAATAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCgg                                                      <  2:559577/116‑1 (MQ=255)
                                                                                     aTAAAACGTTTTCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTa                             <  2:1482778/139‑1 (MQ=255)
                                                                                              tttCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCa                                                                                                                      <  1:1083461/41‑1 (MQ=255)
                                                                                              tttCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCa                                                                                                                      >  2:1083461/1‑41 (MQ=255)
                                                                                              tttCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCttg                                                                                              <  2:384243/65‑1 (MQ=255)
                                                                                              tttCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCttg                                                                                              >  1:384243/1‑65 (MQ=255)
                                                                                                tCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAg                                                                                                             >  1:810323/1‑48 (MQ=255)
                                                                                                tCGGATTGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAg                                                                                                             <  2:810323/48‑1 (MQ=255)
                                                                                                      tGCCATCTTTACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCAtt            >  1:746099/1‑139 (MQ=255)
                                                                                                              ttACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAg    <  2:359731/139‑1 (MQ=255)
                                                                                                               tACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGgc                   <  1:338256/123‑1 (MQ=255)
                                                                                                               tACCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGgc                   >  2:338256/1‑123 (MQ=255)
                                                                                                                aCCCTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGct  >  1:972045/1‑139 (MQ=255)
                                                                                                                  ccTTTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGt                         >  1:1003946/1‑114 (MQ=255)
                                                                                                                     |                                                                                                                                     
AACACATTCAACAATATCAATATTATTTTCGCGCCACATCGCCATTGAATCCTGCGCCTGACAGGCCGGGCATTTTGCGCCCGCAATAAAACGTTTTCGGATTGCCATCTTTACCCTCTACTCAAATTCATCCCAGCCGTCCAGCTGGCGTCGTTCTTGTTGCATTTCACGCTGGAAAATCTCTTCCAGTTCGCGTCGGGCTTCCCTGGCGCGAGAAATTTGTACGGTATCGGCATGCATTGGGATGAGCT  >  NC_000913/3478685‑3478935

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