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.94.4% 141.5 / 1.2 36intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (19/15);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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.

CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740657‑4740948
                                                                                                                                                    |                                                                                                                                               
cGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                 >  8:169602/1‑149 (MQ=255)
cGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                 >  5:158385/1‑149 (MQ=255)
 gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCCt                                                                                                                                                 <  2:19333/148‑1 (MQ=255)
   aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                              <  5:57897/149‑1 (MQ=255)
    cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGCATGCCCTCGGGCAGCCttttt                                                                                                                                             >  6:155687/1‑149 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                          >  2:268181/1‑147 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                          <  3:31445/149‑3 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           <  4:329570/148‑2 (MQ=255)
        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                         <  6:158385/149‑4 (MQ=255)
         aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                        >  5:296305/1‑145 (MQ=255)
         aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                        <  4:345429/149‑5 (MQ=255)
            aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                     <  7:169602/149‑8 (MQ=255)
                      cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGcc                                                                                                                             <  6:296305/146‑1 (MQ=255)
                       aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                              <  7:263125/129‑1 (MQ=255)
                         cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                       >  3:350610/1‑149 (MQ=255)
                                   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggc                                                                                                              >  4:18973/1‑146 (MQ=255)
                                    tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                            <  2:349122/149‑1 (MQ=255)
                                      ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                           <  1:268181/148‑1 (MQ=255)
                                         cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                       >  3:90952/1‑149 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                   <  3:180302/149‑1 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                    <  1:91866/148‑1 (MQ=255)
                                               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                 <  5:155687/149‑1 (MQ=255)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                <  3:18973/149‑1 (MQ=38)
                                                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                             <  4:90952/149‑1 (MQ=255)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                >  5:202456/1‑138 (MQ=255)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                <  6:202456/138‑1 (MQ=255)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCACGCGAAGAGGCCGGCACAgc                                                                                     >  8:80706/1‑148 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                 <  7:153787/127‑1 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                 >  8:153787/1‑127 (MQ=255)
                                                                          cgcgACAGCGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCC‑TTTTTAATTTGACCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                      <  3:226276/149‑1 (MQ=255)
                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                 <  3:163833/106‑1 (MQ=25)
                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                 >  4:163833/1‑106 (MQ=25)
                                                                                            aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTATACGCCCGCt                                                     >  8:301340/1‑148 (MQ=255)
                                                                                               caCCGAAAGTAAGCGGCAGTAGCCTAACGAAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                 <  8:116639/148‑1 (MQ=255)
                                                                                                       gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGACCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                         >  2:54581/1‑149 (MQ=255)
                                                                                                        tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCttt                                        >  7:198557/1‑149 (MQ=255)
                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCCCAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                    >  5:352618/1‑149 (MQ=255)
                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                    <  8:198557/149‑1 (MQ=255)
                                                                                                               cAGTAGCCTACCGCATCAAGGCGGCCCTCGGGCAG‑CTTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTa                                 <  4:350610/149‑1 (MQ=255)
                                                                                                                             aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCCCCACAAAACCTCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTCAGGGGATTGCACGGCGCATc                   >  7:56001/1‑149 (MQ=255)
                                                                                                                               gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTACGCCCGCACCCGCAGCCCTTTCAGGGTATTGAACGGCGCCTcca                 >  4:255224/1‑149 (MQ=255)
                                                                                                                                     tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCAAGACACCACAAAACCTCAGTTTTTCTCAGCCCGCTAACGCAGCCCTTTAAGGGTATTGAACGGCGCATCCACCACg            >  7:249603/1‑148 (MQ=255)
                                                                                                                                     tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCACTTTTTATCCCCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg            >  1:198702/1‑148 (MQ=255)
                                                                                                                                       cccTCGGTCAGGC‑TTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCAGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc         <  8:241129/149‑1 (MQ=255)
                                                                                                                                         cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCATGCCCGCGAAGAGGCCGGCACAAAAACCACAAAACCACACATTTTCTCCGCCCGCACCCGCAGCCCATTAAGGGTATTGAACCGCGCATCCACCACGAACta       >  2:285735/1‑148 (MQ=255)
                                                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg      >  8:149691/1‑149 (MQ=255)
                                                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGAAACCACAAAACAACAGTTTTTATACCCCCCCTCACGCACCCCGTTCAGTGTATTTAACGGCGCATc                   >  8:221383/1‑136 (MQ=255)
                                                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTATGCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                   <  7:221383/136‑1 (MQ=255)
                                                                                                                                              gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTCGCGGACATGCCCCCGAAGCGGACGGCACACCCACCACAAAACCACACTTTTTCTCCCCCCGCTAACGCAGCCCTTTCAGGGTAATTAACGGCGCATCCACCACGAACTTGATgg  >  2:100556/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740657‑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: ‑