Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,463 T→A 8.9% intergenic (‑62/+52) yihR ← / ← yihS putative aldose 1‑epimerase YihR/sulfoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4630TA 8.9% 163.8 / 6.2 56intergenic (‑62/+52)yihR/yihSputative aldose 1‑epimerase YihR/sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base A (3/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

CGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGT  >  NC_000913/4070314‑4070609
                                                                                                                                                     |                                                                                                                                                  
cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                    <  1:1185700/150‑1 (MQ=255)
  tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                  <  1:622533/150‑1 (MQ=255)
    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                >  1:481246/1‑150 (MQ=255)
     tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                                                                               <  1:293/150‑1 (MQ=255)
               cgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                                                                     <  2:1572190/150‑1 (MQ=255)
                ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTACCCTCTCCCCAAAGGGGCGAggggaccg                                                                                                                                    <  1:41081/150‑6 (MQ=255)
                  cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                                                                                                                                  <  2:125718/150‑1 (MQ=255)
                    aGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGgaga                                                                                                                                >  1:297322/1‑150 (MQ=255)
                     gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                                                                               <  1:566170/150‑1 (MQ=255)
                     gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                                                                               >  2:1131276/1‑150 (MQ=255)
                       gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                                                                             >  2:1087656/1‑150 (MQ=255)
                            ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                                                                        <  2:178327/150‑1 (MQ=255)
                            ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                                                                        <  1:1131276/150‑1 (MQ=255)
                             gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                       <  2:1039710/150‑1 (MQ=255)
                                 gTCCACTGCAATGCATATTGGTTATTTGCATACAAGGCACCTTAGTTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                                                                   <  1:1018687/150‑1 (MQ=255)
                                 gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                                                                   <  2:1122741/150‑1 (MQ=255)
                                 gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                                                                   >  1:191760/1‑150 (MQ=255)
                                 ctaaaCTGCAATACATATTGGTGATTTGCATACAAGGAAAATTAATTAACGACATACAGCATACCGTGAGGAAGATGTCAAACACAAAAACTACACCGGACGGTCCCCTCGACCCTTTGTGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                                                                   >  2:1012748/5‑150 (MQ=255)
                                   ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                                                                                 >  2:1579925/1‑150 (MQ=255)
                                    cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                                                                                >  1:87860/1‑150 (MQ=255)
                                          aaTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAACCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAttt                                                                                                          >  2:1553648/1‑150 (MQ=255)
                                           aTGCATATTGGTGATTTGCATACAAGGCACCTTAATTACCGACATACCGCCGCCCGGGCGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                                                                         <  1:1331493/150‑1 (MQ=39)
                                              cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg                                                                                                      <  2:997416/150‑1 (MQ=255)
                                                               tACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                                                                                     >  1:1247327/1‑150 (MQ=255)
                                                                      cACCTTAATTAACGACATACAGCCTCCCGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATcc                                                                              <  2:1456978/150‑1 (MQ=255)
                                                                          ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                                                                          <  2:87860/150‑1 (MQ=255)
                                                                               taaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAccc                                                                     >  1:987473/1‑150 (MQ=255)
                                                                                   gACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg                                                                 <  2:481246/150‑1 (MQ=255)
                                                                                      aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                              >  1:19814/1‑150 (MQ=255)
                                                                                       tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAc                                                             >  1:640131/1‑150 (MQ=255)
                                                                                       tACAGCCTCCCGAGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAACCCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAc                                                             >  2:744515/1‑150 (MQ=255)
                                                                                              tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGa                                                      <  2:19814/150‑1 (MQ=255)
                                                                                                        cTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCgg                                            >  2:1169203/1‑150 (MQ=255)
                                                                                                         tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCggg                                           >  2:1023067/1‑150 (MQ=255)
                                                                                                         tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCggg                                           <  1:1023067/150‑1 (MQ=255)
                                                                                                               cAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTaa                                     >  2:235259/1‑150 (MQ=255)
                                                                                                                  aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                          >  1:884839/1‑46 (MQ=255)
                                                                                                                  aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                          <  2:884839/46‑1 (MQ=255)
                                                                                                                   cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACggg                                 <  2:1215237/150‑1 (MQ=255)
                                                                                                                   cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACggg                                 >  1:380547/1‑150 (MQ=255)
                                                                                                                   cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACggg                                 >  2:364114/1‑150 (MQ=255)
                                                                                                                           cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCgg                         <  2:1247327/150‑1 (MQ=255)
                                                                                                                             gCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGa                       <  2:640131/150‑1 (MQ=255)
                                                                                                                                  ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCAcc                  <  1:1169203/150‑1 (MQ=255)
                                                                                                                                    ccctaacccTCTCCCCAAAGGGGCGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAg                >  1:1255842/7‑150 (MQ=255)
                                                                                                                                    aCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAg                <  1:1087656/150‑1 (MQ=255)
                                                                                                                                     cGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg               >  1:1413629/1‑150 (MQ=255)
                                                                                                                                        aacccTCTCCCCAAAGGGGCGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATc                     >  2:484216/3‑141 (MQ=255)
                                                                                                                                        aacccTCTCCCCAAAGGGGCGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATc                     <  1:484216/139‑1 (MQ=255)
                                                                                                                                         acccTCTCCCCAAAGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAAt           >  2:697133/2‑150 (MQ=255)
                                                                                                                                             tCGCCCCTTTTGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcag       >  2:503133/1‑150 (MQ=255)
                                                                                                                                              cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcagc      >  2:1551455/1‑150 (MQ=255)
                                                                                                                                              cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcagc      <  1:364114/150‑1 (MQ=255)
                                                                                                                                                ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGcagcag    >  1:646898/1‑150 (MQ=255)
                                                                                                                                                  ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATATGCGGGATc                     <  2:169290/131‑1 (MQ=255)
                                                                                                                                                  ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATATGCGGGATc                     >  1:169290/1‑131 (MQ=255)
                                                                                                                                                  ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGt  >  1:32528/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGT  >  NC_000913/4070314‑4070609

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