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.7% 100.2 / ‑3.2 30coding (1034/1245 nt)sstTsodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (13/16);  total (13/17)
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.G86.7% 88.1 / 7.3 30coding (1034/1245 nt)sstTsodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (12/14);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑GGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAA  >  NC_000913/3240842‑3241096
                                                                                                                                        ||                                                                                                                          
tGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑ggg                                                                                                                        <  2:701802/139‑1 (MQ=255)
  gcgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGAACTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑ggggg                                                                                                                      >  1:797836/1‑139 (MQ=255)
   cgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggggg                                                                                                                     >  2:1085795/1‑139 (MQ=255)
     ggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGACGCATCCGGCGTGGC‑‑‑ggggg                                                                                                                      <  2:722998/136‑1 (MQ=255)
     ggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGACGCATCCGGCGTGGC‑‑‑ggggg                                                                                                                      >  1:722998/1‑136 (MQ=255)
      gcgcAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCATCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggggg‑‑ggg                                                                                                                >  1:818804/1‑137 (MQ=255)
                    aTTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGG‑‑GTCTCTGCTGCTGa                                                                                                     >  1:992516/1‑139 (MQ=255)
                      tACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggggg‑‑gg                                                                                                                 >  1:870223/1‑121 (MQ=255)
                      tACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑gggggg‑‑gg                                                                                                                 <  2:870223/122‑2 (MQ=255)
                        ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  2:973597/139‑1 (MQ=255)
                        ccGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  2:797836/139‑1 (MQ=255)
                         cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  2:992516/139‑1 (MQ=255)
                         cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  2:990245/138‑1 (MQ=255)
                         cGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                >  1:990245/1‑138 (MQ=255)
                          gtgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACt                                                                                              <  2:1261875/139‑1 (MQ=255)
                          gtgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACt                                                                                              <  1:907559/139‑1 (MQ=255)
                            gtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCG‑‑GGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           <  1:1085795/139‑1 (MQ=255)
                             ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           <  1:1352866/139‑1 (MQ=255)
                                      gctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAatat                                                                                  >  2:703176/1‑139 (MQ=255)
                                                     cTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGaa                                                                   <  2:1176411/139‑1 (MQ=255)
                                                       gggTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACg                                                                 >  2:1299796/1‑139 (MQ=255)
                                                           aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTc                                                                      <  2:1255571/130‑1 (MQ=255)
                                                           aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTc                                                                      >  1:1255571/1‑130 (MQ=255)
                                                                  tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTg                                                                                             <  1:875047/100‑1 (MQ=255)
                                                                  tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTg                                                                                             >  2:875047/1‑100 (MQ=255)
                                                                  tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           <  2:1222976/102‑1 (MQ=255)
                                                                  tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           >  1:1222976/1‑102 (MQ=255)
                                                                    gATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                >  1:770955/1‑95 (MQ=255)
                                                                    gATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  2:770955/95‑1 (MQ=255)
                                                                                    tgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           <  2:1283466/84‑1 (MQ=255)
                                                                                    tgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGc                                                                                           >  1:1283466/1‑84 (MQ=255)
                                                                                    tgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGc                                    >  2:723538/1‑139 (MQ=255)
                                                                                           gAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCCTGCAGGTGGTTGCCGTCGGCTTTatca                             >  1:1071278/1‑139 (MQ=255)
                                                                                           gAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatca                             >  2:1247356/1‑139 (MQ=255)
                                                                                                            tgtgtgCCTGTGGCGCATCCGGCGTGGC‑‑‑GGGGGGGTGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGAc            <  2:153922/139‑1 (MQ=255)
                                                                                                             gtgtgCCTGTGGCGCATCCGGCGTGGCG‑‑GGGGGG‑‑GTCTCTGCTGCTGATCCCa                                                                                                <  1:387932/53‑1 (MQ=255)
                                                                                                             gtgtgCCTGTGGCGCATCCGGCGTGGCG‑‑GGGGGG‑‑GTCCCTGCTGCTTATcccc                                                                                                >  2:387932/1‑52 (MQ=255)
                                                                                                               gtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGActct         >  1:1309227/1‑139 (MQ=255)
                                                                                                                      tgGCGCATCCGGCGTGGCGGGGGGGG‑‑GTCGCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGaaa  <  1:723538/139‑1 (MQ=255)
                                                                                                                                        ||                                                                                                                          
TGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑‑GGGGGG‑‑GTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAA  >  NC_000913/3240842‑3241096

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