Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,805 Δ1 bp 100% intergenic (+42/+711) 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.91.2% 124.6 / 4.6 34intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (19/12);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740661‑4740948
                                                                                                                                                |                                                                                                                                                
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                                                                                                                                      >  4:27703/1‑142 (MQ=255)
        aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                      >  3:259107/1‑142 (MQ=255)
                tccagcccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                              >  8:19732‑M2/9‑133 (MQ=255)
                tccagcccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                              <  7:19732‑M2/125‑1 (MQ=255)
                     cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                       >  4:150675‑M2/4‑149 (MQ=255)
                               aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                             >  5:272339/1‑149 (MQ=255)
                                        ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           <  1:67042/112‑3 (MQ=255)
                                        ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           >  2:67042/1‑110 (MQ=255)
                                        ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                    <  3:5260/149‑1 (MQ=255)
                                         cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgcg                                                                                                     >  1:162942/1‑147 (MQ=255)
                                          gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggcaaggc                                                                                                         <  7:142200/142‑8 (MQ=255)
                                          gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                  <  5:276343/149‑1 (MQ=255)
                                          gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACAGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggcaaggc                                                                                                         >  8:142200/1‑135 (MQ=255)
                                            gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                >  5:161064/1‑149 (MQ=255)
                                               aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                             <  2:16883/149‑1 (MQ=255)
                                                     acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                                <  8:44293/124‑1 (MQ=255)
                                                      cacCGCCGACCAGACCCGCGACATCGTTCCCCTTAAAGCCCCCAAGTACGCTGCATTAGCCGACCGCAAGAAGGCTGCCCTCGGGCTGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGcaca                                                                                       <  4:196404/148‑1 (MQ=255)
                                                         cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACAAACCGGCCATGCCCGCCAAGAGGCCGCcacacaca                                                                                   >  1:244794/1‑149 (MQ=255)
                                                                  taaGCGTGACTTCGTTGCCCTTAAATACCGAATTTATGCTGCAGTAGCGTAACGAAATATGGGTGCCCTGGGTAAGCTTTTTTTTTTTGTCGTGTGCGGGCCTTTGGGGGGCTTGGCCCGCTATGTGGCCGGCACAGACACCACAAAAc                                                                          <  2:191120/148‑1 (MQ=255)
                                                                              cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAACCCTCAGgttttat                                                              >  8:132394/1‑147 (MQ=255)
                                                                               gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCTAACGCAAGAAGGCTGCCGTCGGGATGCGTTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                              <  7:132394/148‑1 (MQ=255)
                                                                                   cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                 <  4:132404/109‑1 (MQ=255)
                                                                                   cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                 >  3:132404/1‑109 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTGCCGCGAAGAGGCCGGCACAGAAACCACAAAAACACAGTTTTtctagatag                                                        >  7:8265/1‑142 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                          >  2:41262/1‑146 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                          <  1:41262/146‑1 (MQ=255)
                                                                                      tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGCCACCACAAAACCTCAGGTTTTCTccgcccgg                                                      >  6:215386/1‑148 (MQ=255)
                                                                                          aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCAGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                                  >  1:115554/1‑149 (MQ=255)
                                                                                               gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                             <  6:272339/149‑1 (MQ=255)
                                                                                                       gCTGCAGTAGCCTAACGCAAGAAGGCTGCCCGCGGGCAGCCTTTTTTATTTTGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                     <  5:91902/149‑1 (MQ=255)
                                                                                                        cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                    <  7:48804/149‑1 (MQ=255)
                                                                                                                    aaCCCAAGACGGCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTCCCGGCCTGTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                        <  3:27703/149‑1 (MQ=255)
                                                                                                                         aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACAACAGTTTTTATAAGCCCGCTCACGCCGCCCTTTAAg                                     >  2:170877/1‑131 (MQ=37)
                                                                                                                         aaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTAAggg                                   >  7:13920/4‑133 (MQ=37)
                                                                                                                          agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCGCAGCAACCACAAAACCGCAGTTTTTATACGCCCGCGCCCGCAGCCCAGTCAGGGGATTGTAAGgcgccgcc                  >  4:237043/1‑145 (MQ=255)
                                                                                                                           gaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCAGCGAAGAGGCCGGCACACACACCACAAAACCGAAGTTTTTATCCGCCCCCTCACGCAACCCATTCAGGGTATTGAAAGGCGCATcca                 >  7:37394/1‑149 (MQ=255)
                                                                                                                                 tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGCCACCACAAAACCTCAGTGTTTCTCCGCCCGCGCCCGCCCCCCGTTCAGGGGTATGACCGGCGCATCcaccacta           >  8:246562/1‑147 (MQ=255)
                                                                                                                                   cccgagggCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGCAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATCCGCCCGCTCACGCAGCCCTTTAAGGGTATTGAACGGCGCATCCCCCACGAAc         >  3:43366/6‑149 (MQ=255)
                                                                                                                                   cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc         <  2:162942/149‑1 (MQ=255)
                                                                                                                                         gggcatgtttttttATTTTTCCGTGTGCTGCCCTTTGCGGGGCTTGGCCGGCTATGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt       <  2:244794/140‑1 (MQ=255)
                                                                                                                                          gCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTgg  >  3:95010/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740661‑4740948

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: