Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,021 G→A 100% Q113* (CAA→TAA)  fleQ ← transcriptional regulator FleQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,965,0210GA98.1% 163.8 / ‑4.4 51Q113* (CAA→TAA) fleQtranscriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/27);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTA  >  NC_002947/4964876‑4965148
                                                                                                                                                  |                                                                                                                               
tgcATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCt                                                                                                                              >  7:97754/1‑149 (MQ=255)
  catcatCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCt                                                                                                                              <  6:145038/147‑1 (MQ=255)
             gCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTggcgg                                                                                                                  >  5:263165/1‑148 (MQ=255)
                    cTGAATGGCACTACTGGTGCCAACCAGGCTGCGGAACAGATTGGTTTCGCGCTGGCTACCCCGCTCGCGCGACTGGTCATACATTTCGCGGTA‑GACCTGGGCTCGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGTCGTCATTTcc                                                                                                           <  2:120500/148‑1 (MQ=255)
                     tGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAg                                                                                                         >  4:246849/1‑149 (MQ=255)
                     tGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAg                                                                                                         >  7:272275/1‑149 (MQ=255)
                     tGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAg                                                                                                         >  7:271316/1‑149 (MQ=255)
                       aaTGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGt                                                                                                       <  2:153872/149‑1 (MQ=255)
                        aTGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGtt                                                                                                      >  3:97170/1‑149 (MQ=255)
                         tGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTg                                                                                                     <  3:240583/149‑1 (MQ=255)
                             aCGTCTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                 <  6:191977/149‑1 (MQ=255)
                             aCGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATTCTTTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCATTTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                 <  1:100722/149‑1 (MQ=255)
                                  tGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCAcc                                                                                            >  6:291152/1‑149 (MQ=255)
                                    gTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCcg                                                                                          >  8:182382/1‑149 (MQ=255)
                                       ccAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCcggcg                                                                                       >  3:184525/1‑149 (MQ=255)
                                        cAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCcggcgg                                                                                      <  7:182382/149‑1 (MQ=255)
                                        cAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCcggcgg                                                                                      <  8:97754/149‑1 (MQ=255)
                                            cAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGaa                                                                                  <  8:271316/149‑1 (MQ=255)
                                            cAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGaa                                                                                  <  4:97170/149‑1 (MQ=255)
                                             aGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAg                                                                                 <  4:104972/149‑1 (MQ=255)
                                                    ggAACAGCTTGGGGTCGCGCTGGCGACCCCGCTCGCGCGCCGGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTcc                                                                          <  8:272275/149‑1 (MQ=255)
                                                      aaCAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                 >  8:61036/1‑124 (MQ=255)
                                                      aaCAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                 <  7:61036/124‑1 (MQ=255)
                                                           aTTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGc                                                                    <  2:212379/148‑1 (MQ=255)
                                                           aTTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGc                                                                    >  1:212379/1‑148 (MQ=255)
                                                              gggTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTcg                                                                 <  7:10059/148‑1 (MQ=255)
                                                               ggTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTcggc                                                               >  1:280530/1‑149 (MQ=255)
                                                               ggTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTcggc                                                               <  4:184525/149‑1 (MQ=255)
                                                                 ttCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCgg                                                                        <  5:26365/138‑1 (MQ=255)
                                                                 ttCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCgg                                                                        >  6:26365/1‑138 (MQ=255)
                                                                       cTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGaaa                                                       >  3:36351/1‑149 (MQ=255)
                                                                        tGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAAt                                                      >  5:23121/1‑149 (MQ=255)
                                                                          gCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAAttt                                                    <  6:23121/149‑1 (MQ=255)
                                                                               cccGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGc                                                                       >  7:217447/1‑125 (MQ=255)
                                                                               cccGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGc                                                                       <  8:217447/125‑1 (MQ=255)
                                                                               cccGCTCGCGCGCCTCGACATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGAACCCGGCGGCGAAGCTCTTGCGGCAGCGCGGCGGAAGAAATTTCAcc                                                >  6:175903/1‑148 (MQ=255)
                                                                                ccGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACGCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCt                                               >  5:92114/1‑148 (MQ=255)
                                                                                  gCTCGCGCGCCTGGTCATACATTTCGCGGTAGGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTa                                              >  1:37588/1‑148 (MQ=255)
                                                                                   cTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTaaaa                                           <  4:36351/149‑1 (MQ=255)
                                                                                      gcgcgcCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGc                                         <  3:246849/148‑1 (MQ=255)
                                                                                             tGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCg                                  >  6:189543/1‑148 (MQ=255)
                                                                                               gTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGa                                                                                   >  3:114786/1‑97 (MQ=255)
                                                                                               gTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGa                                                                                   <  4:114786/97‑1 (MQ=255)
                                                                                                   tACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCggaaggaa                           <  2:37588/149‑1 (MQ=255)
                                                                                                   tACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCggaaggaa                           <  5:189543/149‑1 (MQ=255)
                                                                                                      aTTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACt                         >  1:187506/1‑148 (MQ=255)
                                                                                                       tttCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCa                       <  6:263165/149‑1 (MQ=255)
                                                                                                                   aCCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGc            >  8:281285/1‑149 (MQ=255)
                                                                                                                     cTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCACAGCCGGAAGGAACTCATCCCAGGTAGccc          >  2:189017/1‑148 (MQ=255)
                                                                                                                        ggCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAATCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACt        <  2:280530/148‑1 (MQ=255)
                                                                                                                             ggTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTTGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTa  <  5:291152/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                               
TGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTA‑GACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTA  >  NC_002947/4964876‑4965148

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