Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,811 (T)7→6 intergenic (+48/+705) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8050T.87.5% 125.5 / 12.5 40intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base . (19/16);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740660‑4740947
                                                                                                                                                 |                                                                                                                                              
aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                             <  2:52876/149‑1 (MQ=255)
 cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                            >  5:204395/1‑149 (MQ=255)
 cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                            <  2:297499/149‑1 (MQ=255)
    gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                         >  4:87967/1‑147 (MQ=255)
    gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                         <  3:87967/149‑3 (MQ=255)
     aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                        >  7:168198/1‑146 (MQ=255)
         aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                    >  3:225480/1‑142 (MQ=255)
         aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                    >  5:71134/1‑142 (MQ=255)
                            cAAGATCGGCCGTGCGCGCCAGCTGTTCACCGGCGAGCAGAAGCGCTACATCGTTGCCCTGAAATACCGCAAGTAAGCTGGAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTGTCCTGTGCCGGCCTCTTcg                                                                                                                  <  1:249875/148‑1 (MQ=255)
                             aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAACCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                               <  5:365905/149‑1 (MQ=255)
                                          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                 <  8:303234/134‑1 (MQ=255)
                                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                <  7:384555/149‑1 (MQ=38)
                                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                <  7:262939/149‑1 (MQ=38)
                                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                <  6:204395/149‑1 (MQ=255)
                                             gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                               <  4:217861/149‑1 (MQ=38)
                                             gcCAGCTGTACACCGGCGAGCAGAAGCGAGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               >  5:428696/1‑149 (MQ=255)
                                               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                             <  2:106794/133‑1 (MQ=255)
                                               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                             >  1:106794/1‑133 (MQ=255)
                                                    gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                 <  7:223352/140‑1 (MQ=38)
                                                    gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                 >  8:223352/1‑140 (MQ=38)
                                                     tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTATTAGCGGCCATGCCCGCGAAGAGGGCGGcac                                                                                       >  1:327653/1‑149 (MQ=255)
                                                       cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                        >  1:88527/1‑146 (MQ=255)
                                                                  ataaGCGCTACATGTTTGCCTTGAATGACCTCAAGTATGGTGCATTAGCCTAACGCAATATGGCTGCCCTGGGGCTGCCTTTTTTTTTTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCGGCTATGTGGCCGGCACAGACACCACaaa                                                                           <  2:88527/146‑1 (MQ=255)
                                                                  aGAAGCCCTACTGCGTGGCCCTTAACCACCGAATGTAAGCCGCAGTCGCCCAACGCACGATGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGGGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACaaaa                                                                          <  4:225480/149‑1 (MQ=255)
                                                                      gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCCCGGCCGTTCCCGCGAAGCGGCCCGCACAGGCACCACAAAAACTc                                                                      >  5:88874/1‑149 (MQ=255)
                                                                        gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                               <  4:103919/106‑1 (MQ=255)
                                                                        gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                               >  3:103919/1‑106 (MQ=255)
                                                                                gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                          >  5:346863/1‑135 (MQ=39)
                                                                                gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                          <  6:346863/135‑1 (MQ=255)
                                                                                 ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                           <  6:428696/149‑1 (MQ=255)
                                                                                         aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTAGGGCAGCCTTTTTTATTTGTCCTGTGCAGGCCACTTCGCCCCCTTTACCGCGAAGCGGCCGGCACAGACACCACAAAACCTCACATTTTATAAGCCCGCt                                                    >  2:4310/1‑148 (MQ=255)
                                                                                              ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                  <  8:8268/97‑1 (MQ=25)
                                                                                              ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                  >  7:8268/1‑97 (MQ=25)
                                                                                                    gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCCGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCCGcccgt                                        >  6:232870/1‑147 (MQ=255)
                                                                                                        gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                    >  8:250588/1‑149 (MQ=255)
                                                                                                        gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                    >  1:159703/1‑149 (MQ=255)
                                                                                                             aGTAGCCGAACGATAGAAGGCTGCCCTGGGGCTGCTTTTTTTTTTTGTCCTGCGCCGGCCTCTTCGCGGCCTTGCCGGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTa                                <  7:344420/148‑1 (MQ=255)
                                                                                                                  ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                              >  7:85645/1‑113 (MQ=38)
                                                                                                                  ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                              <  8:85645/113‑1 (MQ=39)
                                                                                                                    gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgg                        >  1:81623/1‑149 (MQ=255)
                                                                                                                      ctatctgaagGCGGCCCTCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                      <  1:355543/143‑1 (MQ=255)
                                                                                                                          aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                  >  7:292946/1‑149 (MQ=255)
                                                                                                                          aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCATTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACACTTTTTATCCGCCCCCTCACCCACCCCATTAAGGGTATTTAAAGGCGCaaa                  >  2:9834/1‑147 (MQ=255)
                                                                                                                             aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac               >  1:256956/1‑149 (MQ=255)
                                                                                                                                  tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGa          <  1:9834/149‑1 (MQ=255)
                                                                                                                                     ccTGGGGCAGCTTTTTTTTTTTGCCGTGTGCCGGCCTCTTCGCGGCGTTGCCCGCTATGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                  <  4:354239/122‑1 (MQ=255)
                                                                                                                                       tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg     >  1:319/1‑149 (MQ=255)
                                                                                                                                          ggCAGCCTTTTTTTTTTGCCCGGTGCCGGCCTCTTCGGGGCCTTCCCCGCGAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg  <  6:88874/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740660‑4740947

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