Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,795 T→G 10.9% intergenic (+32/+721) 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,7950TG10.9% 104.1 / 4.5 45intergenic (+32/+721)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base G (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAC  >  NC_002947/4740652‑4740941
                                                                                                                                               |                                                                                                                                                  
GATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCggg                                                                                                                                                >  6:193569‑M1/1‑148 (MQ=255)
         cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGCTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                      <  1:143061/149‑1 (MQ=255)
            gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                   >  8:48804/1‑147 (MQ=255)
            gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCCGCCCTCGGGCAGCCttttttat                                                                                                                                   <  7:41632/149‑3 (MQ=255)
             tccagtgcgagatccagcccGTACAAGAGCGGCAGTCCGCGCCTGCGGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCATGTAAGCTGCAGTAGCCGAACGCAAGCAGGCTGCCCTCGGGCAGccttttttat                                                                                                                                   <  4:47223‑M2/128‑3 (MQ=255)
             aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGTCTGCCCTCGGGCAGCCttttttatt                                                                                                                                  <  5:164356/149‑4 (MQ=255)
                gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgt                                                                                                                               <  5:230714/149‑7 (MQ=255)
                 aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                              >  3:259107/1‑142 (MQ=255)
                 aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                              >  4:27703/1‑142 (MQ=255)
                         tccagcccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                      >  8:19732‑M2/9‑133 (MQ=255)
                         tccagcccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                      <  7:19732‑M2/125‑1 (MQ=255)
                              cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                >  4:150675‑M2/4‑149 (MQ=255)
                                        aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                      >  5:272339/1‑149 (MQ=255)
                                                 ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                   >  2:67042/1‑110 (MQ=255)
                                                 ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                   <  1:67042/112‑3 (MQ=255)
                                                 ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                             <  3:5260/149‑1 (MQ=255)
                                                  cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgcg                                                                                              >  1:162942/1‑147 (MQ=255)
                                                   gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggcaaggc                                                                                                  <  7:142200/142‑8 (MQ=255)
                                                   gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                           <  5:276343/149‑1 (MQ=255)
                                                   gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACAGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggcaaggc                                                                                                  >  8:142200/1‑135 (MQ=255)
                                                     gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                         >  5:161064/1‑149 (MQ=255)
                                                        aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                      <  2:16883/149‑1 (MQ=255)
                                                              acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                         <  8:44293/124‑1 (MQ=255)
                                                               cacCGCCGACCAGACCCGCGACATCGTTCCCCTTAAAGCCCCCAAGTACGCTGCATTAGCCGACCGCAAGAAGGCTGCCCTCGGGCTGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGcaca                                                                                <  4:196404/148‑1 (MQ=255)
                                                                  cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACAAACCGGCCATGCCCGCCAAGAGGCCGCcacacaca                                                                            >  1:244794/1‑149 (MQ=255)
                                                                           taaGCGTGACTTCGTTGCCCTTAAATACCGAATTTATGCTGCAGTAGCGTAACGAAATATGGGTGCCCT‑GGGTAAGCTTTTTTTTTTTGTCGTGTGCGGGCCTTTGGGGGGCTTGGCCCGCTATGTGGCCGGCACAGACACCACAAAAc                                                                   <  2:191120/148‑1 (MQ=255)
                                                                                       cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAACCCTCAGgttttat                                                       >  8:132394/1‑147 (MQ=255)
                                                                                        gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCTAACGCAAGAAGGCTGCCGTCGGG‑ATGCGTTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                       <  7:132394/148‑1 (MQ=255)
                                                                                            cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                          <  4:132404/109‑1 (MQ=255)
                                                                                            cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                          >  3:132404/1‑109 (MQ=255)
                                                                                              cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTGCCGCGAAGAGGCCGGCACAGAAACCACAAAAACACAGTTTTtctagatag                                                 >  7:8265/1‑142 (MQ=255)
                                                                                              cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                   <  1:41262/146‑1 (MQ=255)
                                                                                              cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                   >  2:41262/1‑146 (MQ=255)
                                                                                               tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGCCACCACAAAACCTCAGGTTTTCTccgcccgg                                               >  6:215386/1‑148 (MQ=255)
                                                                                                   aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCAGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                           >  1:115554/1‑149 (MQ=255)
                                                                                                        gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                      <  6:272339/149‑1 (MQ=255)
                                                                                                                gCTGCAGTAGCCTAACGCAAGAAGGCTGCCCGCGGGCAGCC‑TTTTTTATTTTGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                              <  5:91902/149‑1 (MQ=255)
                                                                                                                 cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                             <  7:48804/149‑1 (MQ=255)
                                                                                                                             aaCCCAAGACGGCTGCCCTCGGGCAGCCTTTTTTT‑TTTGTCCTGTCCCGGCCTGTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                 <  3:27703/149‑1 (MQ=255)
                                                                                                                                  aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACAACAGTTTTTATAAGCCCGCTCACGCCGCCCTTTAAg                              >  2:170877/1‑131 (MQ=37)
                                                                                                                                  aaaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTAAggg                            >  7:13920/4‑133 (MQ=37)
                                                                                                                                   agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCGCAGCAACCACAAAACCGCAGTTTTTATACGCCCGCGCCCGCAGCCCAGTCAGGGGATTGTAAGgcgccgcc           >  4:237043/1‑145 (MQ=255)
                                                                                                                                    gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCAGCGAAGAGGCCGGCACACACACCACAAAACCGAAGTTTTTATCCGCCCCCTCACGCAACCCATTCAGGGTATTGAAAGGCGCATcca          >  7:37394/1‑149 (MQ=255)
                                                                                                                                          tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGCCACCACAAAACCTCAGTGTTTCTCCGCCCGCGCCCGCCCCCCGTTCAGGGGTATGACCGGCGCATCcaccacta    >  8:246562/1‑147 (MQ=255)
                                                                                                                                            cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc  <  2:162942/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAC  >  NC_002947/4740652‑4740941

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: ‑