Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,183,190 C→T A214V (GCG→GTG)  zupT → heavy metal divalent cation transporter ZupT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,1900CT96.3% 85.4 / NA 27A214V (GCG→GTG) zupTheavy metal divalent cation transporter ZupT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (1/25);  minor base G (1/0);  total (2/25)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGG  >  NC_000913/3183055‑3183313
                                                                                                                                       |                                                                                                                            
ccGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGGTGCAGGGATCATGGGGGGGctc                                                                                                                          >  2:818380/1‑140 (MQ=255)
  cagATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCtctc                                                                                                                        <  1:607317/138‑1 (MQ=255)
   cgATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGGCGGGATTATTGTGGGGctcttg                                                                                                                       >  2:584834/1‑138 (MQ=255)
    gATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGGTTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCgg                                                                                                                      <  2:743706/140‑1 (MQ=255)
                ggggATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGGAGGAAATATGGGGGGGGCCTCGGGTGGTGAAtttt                                                                                                         >  2:607317/1‑138 (MQ=255)
                  ggATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCCGTTTCAGGGAATTTGGGGGGGCTCTCGGGTGATGAAATTATg                                                                                                       >  2:916096/1‑140 (MQ=255)
                      ttCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGc                                                                                                   <  1:51939/140‑1 (MQ=255)
                        cTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTc                                                                                                 <  2:657902/140‑1 (MQ=255)
                              tAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCaaa                                                                                           >  1:173128/1‑140 (MQ=255)
                              tAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCaaa                                                                                           <  1:402868/140‑1 (MQ=255)
                                gCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGGTGCAGGGATTATGGGGGGGGTCTCGGGTGATGAAATTAATCCGCCCGCCcaaaa                                                                                         >  2:402868/1‑138 (MQ=255)
                                    aaaTTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGAtt                                                                                     <  2:917260/140‑1 (MQ=255)
                                    aaaTTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGAtt                                                                                     <  1:749107/140‑1 (MQ=255)
                                    aaaTTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGAtt                                                                                     <  1:231585/140‑1 (MQ=255)
                                      aTTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGa                                                                                   <  1:710073/140‑1 (MQ=255)
                                        tCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGAcc                                                                                 <  1:551675/140‑1 (MQ=255)
                                        tCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGAcc                                                                                 <  1:472390/140‑1 (MQ=255)
                                         cTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGGGATTGAccc                                                                                <  1:798532/140‑1 (MQ=255)
                                         cTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGAccc                                                                                <  1:221792/140‑1 (MQ=255)
                                              tggtgTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGGTGCAGGGAATATGGTGGTGCTCTCGGGTGATGAAATAAAGCCGCCCGCCAAAAAGAATTAACCcaata                                                                           >  2:710073/1‑140 (MQ=255)
                                                 tgtgCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTTCAGGGATTATGGGGGTGGTCTCGGGTGGTGAATTTATGCCGCTCGCCCAAGAGATTGGCCCttataata                                                                        >  2:420060/1‑140 (MQ=255)
                                                  gtgCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGGGGTGCTCTCGGGTTATGAATTAATGCCGCTCGCCAAAAAGATTTACCCtaataaaaa                                                                       >  2:789392/1‑140 (MQ=255)
                                                          gTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCt                                                               <  2:140576/140‑1 (MQ=255)
                                                          gTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCt                                                               <  2:250991/140‑1 (MQ=255)
                                                          gTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCt                                                               <  1:420060/140‑1 (MQ=255)
                                                          gTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCt                                                               <  2:387513/140‑1 (MQ=255)
                                                          gTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCt                                                               <  1:611451/140‑1 (MQ=255)
                                                             gTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGGTGCAGGGATTTTTGGGGGGCTCTCGGGTGATGAATTTATGCCGCTCGCCAAAGAGAATGGCCCTAATAATAAACCCCGCTATg                                                            >  2:102804/1‑140 (MQ=255)
                                                                aaTCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCg                                                         <  1:456521/140‑1 (MQ=255)
                                                                aaTCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCg                                                         <  2:590803/140‑1 (MQ=255)
                                                                    cTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACt                                                     <  2:173128/140‑1 (MQ=255)
                                                                               gATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGt                                          <  1:789392/140‑1 (MQ=255)
                                                                                          tggtCATGGCGGCAATCATGGCGGCGGTTGCCGGGATTATGGGGGGGTTCTCGGGTGATGAATTAATGCCGCCCGCCAAAAAGAATTACCCCCAAAATAAACCCACCCCTGGCGGACTATGTGGGATCGCAATGAAgggg                               >  2:116690/1‑139 (MQ=255)
                                                                                          tggtCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATTGGGGGGCTCTCGGGTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAAACCCAGCTATGGCGTACTATGTGGAAAGTCAGTGATgggg                               >  2:636944/1‑139 (MQ=255)
                                                                                            gtCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGAtt                             <  2:361947/140‑1 (MQ=255)
                                                                                                 ggcggcAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGtt                        <  1:636944/140‑1 (MQ=255)
                                                                                                   cggcAATCATGGCGGCGGTTGCAGGAATTATGGTGGTGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTa                      <  2:616325/140‑1 (MQ=255)
                                                                                                    ggcAATCATGGCGGCGGTTGCAGGAATTATGGGGGGGCTCTCGGGTGATGAATTAATGCCGCTCGCCAAAGAGATTGAACCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAg                     >  1:616325/1‑140 (MQ=255)
                                                                                                                       gCAGGAATTATGGGGGGGCTCTCGGGTGATGAATTAATGCCGCTCGCCAAAAAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCggg  >  2:791028/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                            
CCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTCCCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGG  >  NC_000913/3183055‑3183313

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