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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,240,9771.G96.9% 109.7 / ‑3.6 32coding (1034/1245 nt)sstTsodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (14/17);  total (14/18)
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
*NC_0009133,240,9772.G90.6% 101.6 / 3.4 32coding (1034/1245 nt)sstTsodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (14/15);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAA  >  NC_000913/3240843‑3241106
                                                                                                                                         ||                                                                                                                                   
ggcgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggg                                                                                                                                >  2:601108/1‑139 (MQ=255)
 gcgggcgCAGCAATCACTATTACCGTGTTGACCCTCGCCGCGGTTTATACGCTGGGTATTTCGGTCGATCTGCCCACGGCGCTGCCGTTTGGCGTGGTGGGTTCTCTGTGTGCCTTTGGCGCATCCGGCGTGGC‑‑‑ggggg                                                                                                                               >  1:300963/1‑139 (MQ=255)
    ggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggggggg                                                                                                                           <  2:410793/139‑2 (MQ=255)
     gcgcAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑ggggggggg                                                                                                                          >  1:592209/1‑137 (MQ=255)
               cACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCG‑‑GGGGGGGTCTctgctgct                                                                                                                 <  2:430266/139‑1 (MQ=255)
                  tATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGa                                                                                                               <  2:384515/139‑1 (MQ=255)
                       ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  1:703825/139‑1 (MQ=255)
                       ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  2:396064/139‑1 (MQ=255)
                       ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  2:319809/139‑1 (MQ=255)
                       ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCa                                                                                                          <  1:9305/139‑1 (MQ=255)
                          tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTctgct                                                                                                                    <  2:433947/126‑1 (MQ=255)
                          tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTctgct                                                                                                                    >  1:433947/1‑126 (MQ=255)
                           gtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCG‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     <  1:738286/139‑1 (MQ=255)
                           gtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCAGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTgg                                                                                                      <  2:300963/139‑1 (MQ=255)
                            ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGTGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     >  2:520381/1‑139 (MQ=255)
                            ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     >  1:601423/1‑139 (MQ=255)
                            ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     <  1:527011/139‑1 (MQ=255)
                            ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     <  2:601423/139‑1 (MQ=255)
                            ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     >  2:654317/1‑139 (MQ=255)
                             tGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                                                    >  2:334648/1‑139 (MQ=255)
                              gACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGTGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCt                                                                                                   >  1:742890/1‑139 (MQ=255)
                                              aaTACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTa                                                                                    >  1:801924/1‑139 (MQ=255)
                                                     tGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAAc                                                                            >  1:106099/1‑139 (MQ=255)
                                                          aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                       >  2:681274/1‑139 (MQ=255)
                                                            tCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCg                                                                     >  2:201581/1‑139 (MQ=255)
                                                                               gcgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCg                                                                               <  1:611000/110‑1 (MQ=255)
                                                                               gcgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCg                                                                               >  2:611000/1‑110 (MQ=255)
                                                                               gcgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGt                                                  <  2:659783/139‑1 (MQ=255)
                                                                                   tgctgTTTAGCGTGGTGGCTTCTCCGGGTGCCTGTGGCCCATCCGGCTTTGGTG‑‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGc                                              <  1:591447/139‑1 (MQ=255)
                                                                                    gctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCt                                             <  2:580932/139‑1 (MQ=255)
                                                                                                ggtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCg                                 <  1:207585/139‑1 (MQ=255)
                                                                                                  tggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     <  2:148255/69‑1 (MQ=255)
                                                                                                  tggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                                                     >  1:148255/1‑69 (MQ=255)
                                                                                                                     tgGCGCATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGaaa            >  1:690529/1‑139 (MQ=255)
                                                                                                                          cATCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACtgcg       <  2:472476/139‑1 (MQ=255)
                                                                                                                            tCCGGCGTGGCGGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCt     >  2:476269/1‑139 (MQ=255)
                                                                                                                              cGGTGTGGC‑‑‑GGGGGGTCTCTTTGCTGCTCAGCCCACTCGCCTGTAATATTTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTGATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGaa  <  1:449466/139‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                                   
GGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑GGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAA  >  NC_000913/3240843‑3241106

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