Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,610 C→G 13.9% intergenic (+65/‑334) alx → / → sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6100CG13.9% 126.2 / 11.1 43intergenic (+65/‑334)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base G (3/3);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTG  >  NC_000913/3239493‑3239736
                                                                                                                     |                                                                                                                              
ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                           <  1:529658/123‑1 (MQ=255)
ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                           >  2:529658/1‑123 (MQ=255)
 tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgaga                                                                                                                        >  3:431859/1‑125 (MQ=255)
     atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                         <  3:200763/120‑1 (MQ=255)
     atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                         >  4:200763/1‑120 (MQ=255)
         aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                 >  1:48770/1‑124 (MQ=255)
         aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                 <  2:48770/124‑1 (MQ=255)
                       tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTc                                                                                                                                >  4:139215/1‑95 (MQ=255)
                       tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTc                                                                                                                                <  3:139215/95‑1 (MQ=255)
                       tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                   <  4:25937/108‑1 (MQ=255)
                       tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                   >  3:25937/1‑108 (MQ=255)
                         tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                               >  2:979683/1‑94 (MQ=255)
                         tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                               <  1:979683/94‑1 (MQ=255)
                            gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                       <  3:57485/99‑1 (MQ=255)
                            gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                       >  4:57485/1‑99 (MQ=255)
                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                        >  1:990253/1‑65 (MQ=255)
                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                        <  2:990253/65‑1 (MQ=255)
                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                   >  3:322657/1‑70 (MQ=255)
                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                   <  4:322657/70‑1 (MQ=255)
                                                                   aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                <  2:950831/83‑1 (MQ=255)
                                                                   aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                >  1:950831/1‑83 (MQ=255)
                                                                     gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                              >  4:252187/1‑51 (MQ=255)
                                                                     gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                              <  3:252187/51‑1 (MQ=255)
                                                                         tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                 <  1:453659/124‑1 (MQ=255)
                                                                         tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                 >  2:453659/1‑124 (MQ=255)
                                                                               cTGTCATTAACATGAGCCCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                         <  2:582081/94‑1 (MQ=255)
                                                                               cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                         >  1:582081/1‑94 (MQ=255)
                                                                                     ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                 <  1:10480/48‑1 (MQ=255)
                                                                                     ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                 >  2:10480/1‑48 (MQ=255)
                                                                                      tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTcag                                                                                                                 >  3:384171/1‑45 (MQ=255)
                                                                                      tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                 <  4:384171/47‑1 (MQ=255)
                                                                                      tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                  <  2:348898/142‑1 (MQ=255)
                                                                                      tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                  >  1:348898/1‑142 (MQ=255)
                                                                                         cATGAGCACCGTTTTCTCCCTCTCCCTTCCGGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                    <  4:58004/137‑1 (MQ=255)
                                                                                         cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                    >  3:58004/1‑137 (MQ=255)
                                                                                              gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                    <  3:185368/84‑1 (MQ=255)
                                                                                              gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                    >  4:185368/1‑84 (MQ=255)
                                                                                                aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg  >  1:310786/1‑148 (MQ=255)
                                                                                                 ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                              <  3:788292/119‑1 (MQ=255)
                                                                                                 ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                              >  4:788292/1‑119 (MQ=255)
                                                                                                          ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGc                                                                                      <  4:22262/54‑1 (MQ=255)
                                                                                                          ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGc                                                                                      >  3:22262/1‑54 (MQ=255)
                                                                                                            ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                               <  4:1051634/75‑1 (MQ=21)
                                                                                                            ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                               >  3:1051634/1‑75 (MQ=21)
                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                          >  2:615712/1‑79 (MQ=255)
                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                          <  1:615712/79‑1 (MQ=255)
                                                                                                             tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                              <  2:84828/107‑1 (MQ=255)
                                                                                                                     |                                                                                                                              
TTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTG  >  NC_000913/3239493‑3239736

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