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.81.2% 80.1 / 20.1 32intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base . (13/13);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 5.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740658‑4740951
                                                                                                                                                   |                                                                                                                                                  
gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                    >  3:28735/1‑148 (MQ=255)
   cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACGGCAAGTAAGCTGCAGTAGGCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                                 <  1:190713/148‑1 (MQ=255)
    tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                                <  7:1854/148‑1 (MQ=255)
       aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                             <  7:161285/148‑3 (MQ=255)
       aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                               <  1:18045/146‑1 (MQ=255)
       aaagaaaTGCTCTCCAGCCCGTACAAGATAGGCCGTCCGCGCCAGCTGTACAACGGCAAGCAGAAGCGCGACATCGTTGCAATGAAAGAACGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCATCCttttttat                                                                                                                                             >  4:82556/1‑146 (MQ=255)
                      aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                               <  3:73756/146‑1 (MQ=255)
                                             gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                        >  2:84274/1‑130 (MQ=255)
                                             gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                      <  4:77585/148‑1 (MQ=255)
                                              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                    <  2:18662/148‑1 (MQ=38)
                                              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                       <  1:84274/130‑1 (MQ=255)
                                               gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                    <  1:4595/148‑1 (MQ=38)
                                               gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                    <  8:50153/148‑1 (MQ=255)
                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGCgggc                                                                                                >  5:158765/1‑146 (MQ=255)
                                                                     gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGGCTTTGCCGCGAAGCGGGCGGCACAGACACAACaaaa                                                                              >  1:67700/1‑148 (MQ=255)
                                                                          gcgACATGGTTGCCCTTAAACACCGAATTTACGCTGCAGTAGCCTAACCCAATACGGGTGCCCGGGGGCTGCTTTTTTTATTTGTCGTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCa                                                                         <  6:158765/148‑1 (MQ=255)
                                                                          gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCGCTTTGCGGCCGTGGCGGCGGAGGGGGCGGCCGAGCCCCCCCaaaaaccac                                                                         >  1:163488/1‑143 (MQ=255)
                                                                                 cGTTGCCCTGAAAGACCGCTAGTAAGCTGCAGTCGCGGAACGCAATAAGGCTGCCCTCGGGCTGCCTTTTTTTTTTGGGCTGTGCGGGCCTCTTCGCGGCGTTGCCCGGTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                  <  4:55377/148‑1 (MQ=255)
                                                                                        cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCAGGCCTCTTCGCGGCCTTGACCGCGAAGAGGCCGGCACAGACACCACAAAACCTCACTTATTCTacgccc                                                           >  6:172653/1‑148 (MQ=255)
                                                                                           aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                        >  4:40246/1‑148 (MQ=255)
                                                                                              gACCGCAAGTAAGCTGCAGTCGCCGAACGAATGAAGGCTGCCCTCGGGCAGCCTTTTTTTTTTGTCGTGTGCCGGCCTCTTGGCGGCGTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                                     <  6:178821/148‑1 (MQ=255)
                                                                                              gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGAAGAGGCCGGCACAGACACCACAAAACCTCAGGGTTTGTCCGCCCCCTccc                                                     >  5:352/1‑146 (MQ=255)
                                                                                                          gctgctTTAGCCGAACGAATGATGGCTGCCCTCGGGCAGCTTTTTTTTTTTTTCGTGTGCCGGCCTCTTCGGGGCCTTGCCCGCTAATTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  6:352/148‑1 (MQ=255)
                                                                                                          gCTGGAGTAGCCCAACGAAAGACGGGGGCCGTGGGGCTGGTTTTTTTTTTTGTCGTGGGCCGGCCTCTGCGCGGCCTTGCCGGGTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  4:102344/148‑1 (MQ=255)
                                                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACaaa                                                                               >  6:130986/1‑110 (MQ=38)
                                                                                                           cggCAGTCGCCCACCCCACGACGGCGGCCCTCGGCCAGCCTTTTTTATTTGGCCTGGCCCGGCCCCTCCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                              <  5:130986/108‑1 (MQ=16)
                                                                                                            tGCAGTACCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                       <  4:112458/148‑1 (MQ=255)
                                                                                                              cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                      >  1:14669/1‑147 (MQ=255)
                                                                                                               aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGt                                     <  2:14669/147‑1 (MQ=255)
                                                                                                               aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggt                                      <  3:40246/146‑2 (MQ=255)
                                                                                                                 tAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                  <  2:54900/148‑1 (MQ=255)
                                                                                                                          gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                          >  2:116779/1‑147 (MQ=255)
                                                                                                                            aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                       >  4:45311/1‑148 (MQ=255)
                                                                                                                               aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCCGCCCTTTCAGGGTATTGAACGGCGCATccc                    >  7:146082/1‑147 (MQ=255)
                                                                                                                                       ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTATTCGAGGCCTTGCCCCCGAAGAGGCCGGCACAGACACCACAAAACCTAAGTTTTTATCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAACCACCACGAAc            >  2:127506/1‑148 (MQ=255)
                                                                                                                                             gCAGCCTTTTTAATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  4:72530/148‑1 (MQ=255)
                                                                                                                                             gCAGCCTTTTTAATTTGCCTGGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  6:179600/148‑1 (MQ=255)
                                                                                                                                                  cTTTTTTTTTGGGTCTGGGGCCGCCGGCTGCGCGGCGTGCCCCGGCACGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca  <  2:73828/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740658‑4740951

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: