Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,240,984 (G)7→8 coding (1041/1245 nt) sstT → serine/threonine:Na(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,240,9771.G92.9% 93.8 / 0.8 28coding (1034/1245 nt)sstTserine/threonine:Na(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (1/25);  total (2/26)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAA  >  NC_000913/3240844‑3241106
                                                                                                                                      |                                                                                                                                  
gctggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑gggggg                                                                                                                             >  2:913453/4‑140 (MQ=255)
 cgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGGTCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑ggggggg                                                                                                                            <  2:46262/140‑1 (MQ=255)
    gcgcAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑gggggggggt                                                                                                                         <  1:124901/140‑4 (MQ=255)
      gcagcaATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑ggggggggtctc                                                                                                                       <  1:787865/140‑6 (MQ=255)
           aaTCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTctgctg                                                                                                                   <  1:983396/140‑1 (MQ=255)
                cTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGa                                                                                                               <  1:1172063/140‑1 (MQ=255)
                cTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGa                                                                                                               <  1:850210/140‑1 (MQ=255)
                  aTTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATcc                                                                                                            >  1:632062/1‑140 (MQ=255)
                    tACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  2:1224133/140‑1 (MQ=255)
                    tACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  1:873416/140‑1 (MQ=255)
                      ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACt                                                                                                        <  1:366033/140‑1 (MQ=255)
                      ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACt                                                                                                        <  1:936956/140‑1 (MQ=255)
                       cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTg                                                                                                       >  2:1010177/1‑140 (MQ=255)
                       cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTg                                                                                                       <  1:283129/140‑1 (MQ=255)
                       cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCAccg                                                                                                       >  1:40441/1‑138 (MQ=255)
                         tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     >  2:342324/1‑140 (MQ=255)
                         tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGGTCCCCCCGGc                                                                                                     >  2:663752/1‑140 (MQ=255)
                           ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCCCCGGccc                                                                                                   >  2:334443/1‑139 (MQ=255)
                                   ggctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATg                                                                                           >  1:914194/1‑140 (MQ=255)
                                     ctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCCCTTGCCCGTAATATGtt                                                                                         >  1:1184276/1‑140 (MQ=255)
                                                 cGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGaa                                                                             <  2:632062/140‑1 (MQ=255)
                                                       gTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                       <  1:630032/140‑1 (MQ=255)
                                                       gTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                       <  2:558515/140‑1 (MQ=255)
                                                       gTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                       <  1:1034657/140‑1 (MQ=255)
                                                         aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCg                                                                     <  2:710779/140‑1 (MQ=255)
                                                          ttCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCg                                                                     >  1:237890/1‑140 (MQ=255)
                                                              ggTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATg                                                                <  1:663752/140‑1 (MQ=255)
                                                                           aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCg                                                   <  2:685420/140‑1 (MQ=255)
                                                                             ggCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTTCTTCTGATCCCCCTGGCCCGTAATATTTTTTGGATTTTGAACGATATCGCCCTGCAGGGGGGTTCCGtt                                                 >  2:280441/1‑139 (MQ=255)
                                                                              gcgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTTCTGCTGATCCCACCGGCCCGTAAAATGTTCCGTATTTTGGACGGTATCTCCCTGGAGGGGGGTGCCggcc                                                >  2:917366/1‑137 (MQ=255)
                                                                                   gctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCttt                                           <  1:866602/140‑1 (MQ=255)
                                                                                       ttGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑GGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcat                                      >  2:443593/1‑140 (MQ=255)
                                                                                        tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcat                                      <  1:808698/140‑1 (MQ=255)
                                                                                        tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcat                                      <  1:342324/140‑1 (MQ=255)
                                                                                           gCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCg                                    <  1:474512/140‑1 (MQ=255)
                                                                                            cGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGc                                  <  1:1010177/140‑1 (MQ=255)
                                                                                            cGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCg                                 <  1:443593/140‑1 (MQ=255)
                                                                                              tggtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGGCTCTTCTGCTGATCCCACCGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCATGCAGGGGGGTGCCGGCGGCTTTTTTATCgggg                                 >  2:474512/1‑138 (MQ=255)
                                                                                                        tctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGAc                      <  1:334443/140‑1 (MQ=255)
                                                                                                         ctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGAc                      <  1:890234/140‑1 (MQ=255)
                                                                                                          tgtgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACGGTATCGCCCTGCCGGGGGGTGCCGTCGGCTTTTTCATCGGCGTATTGCAGGAcc                     >  2:648421/1‑139 (MQ=255)
                                                                                                                  tgtgGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGGGGGTGCCGGCGGCTTTTTCATCGGCGTATTGCAGGACTCTTGCGaaa            >  2:764971/1‑140 (MQ=255)
                                                                                                                   gtgGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGaaa            <  1:648421/140‑1 (MQ=255)
                                                                                                                    tgGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCCTGCAGGGGGGTGCCGTCGGCTTTTTTATCGGCGTATTGCAGGACTCTTTCGAAAcc          >  2:779342/1‑140 (MQ=255)
                                                                                                                           tCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTg    <  1:306387/140‑1 (MQ=255)
                                                                                                                            ccGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGaa  <  1:913453/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
GCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAA  >  NC_000913/3240844‑3241106

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