Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,464 T→A 6.8% intergenic (‑63/+51) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4640TA 6.8% 170.4 / 3.4 53intergenic (‑63/+51)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (33/13);  new base A (2/2);  total (36/17)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGA  >  NC_000913/4070315‑4070611
                                                                                                                                                     |                                                                                                                                                    
gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                                                                      <  1:1252164/150‑1 (MQ=255)
   aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                   <  1:231894/150‑1 (MQ=255)
          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGTCACCTTAATTAACGACATCCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGG‑CGAggg                                                                                                                                           <  1:111265/150‑1 (MQ=255)
          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGG‑CGAggg                                                                                                                                           <  1:111318/150‑1 (MQ=255)
          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAggg                                                                                                                                           >  1:615624/1‑150 (MQ=255)
            cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGtt                                                                                                                                         >  1:361791/1‑150 (MQ=255)
            cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGtt                                                                                                                                         >  1:361784/1‑150 (MQ=255)
             gcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTa                                                                                                                                        <  1:1158715/150‑1 (MQ=255)
                 cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGga                                                                                                                                    <  2:1083183/150‑1 (MQ=255)
                   aGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGgaga                                                                                                                                  <  2:1345794/150‑1 (MQ=255)
                     tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAGACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAgg                                                                                                                                <  2:254997/150‑1 (MQ=255)
                      gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAggg                                                                                                                               <  2:723668/150‑1 (MQ=255)
                      gAAACGGTCTGTCCACGGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAggg                                                                                                                               <  2:887260/150‑1 (MQ=255)
                        aaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCCAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGGGGGGG‑AGAGGGGGAGGGAGAggggc                                                                                                                             >  2:740184/1‑149 (MQ=255)
                             tCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAAcc                                                                                                                        >  1:1283344/1‑150 (MQ=255)
                             tCTGACAACTGCAATGCGTATTGGTGATTTGCATACAAGGCGCCTGGTTTGAGGACTGGCAGCCTCCCGGGTGGCTGTAGTCAGACCCGCGGGCTGCAGGGGGCGGTCCGGTCGCGCGTGTGGGG‑TGTGGGGTTGGGGGGGGGGacctcg                                                                                                                        >  2:942341/1‑145 (MQ=255)
                                 tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGACAACCCGCa                                                                                                                    >  2:1069581/1‑150 (MQ=255)
                                 tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCa                                                                                                                    >  2:1342023/1‑150 (MQ=255)
                                   cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                                                                                  >  1:713300/1‑150 (MQ=255)
                                   cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                                                                                  >  1:841520/1‑150 (MQ=255)
                                        cAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt                                                                                                             >  1:535361/1‑150 (MQ=255)
                                                atTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATc                                                                                                     >  1:204309/1‑150 (MQ=255)
                                                   ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTCGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                                                                                  >  1:113868/1‑150 (MQ=255)
                                                    gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGtt                                                                                                 >  1:595492/1‑150 (MQ=255)
                                                      gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAt                                                                                               >  2:1395718/1‑150 (MQ=255)
                                                           gCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                                                                                          >  2:182207/1‑150 (MQ=255)
                                                                 aaGGCACCTTAATTAACCACATACATCCACCCGTGAGGCAGATTTCAAACACAAAAACTGCACCGGACGGTACCCTCGCACATTTGCGG‑AGAGGGTTAGGGAGAGGGGAAAACCAGCACCCTATTGTAAGATAAGTTATTTCGCATTAat                                                                                    >  2:624112/1‑150 (MQ=255)
                                                                        cTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagc                                                                             >  2:872171/1‑150 (MQ=255)
                                                                         ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                                                                            >  2:913424/1‑150 (MQ=255)
                                                                         ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                                                                            >  2:1338543/1‑150 (MQ=255)
                                                                           aattacCGCCTTACCGCCTCCCGTGTGGCTGTTGTCAAACCCAAAAACGGCCCCGGACGGACCCCCCGCCCCTTTGGGG‑TGCGGGTTAGGGAGTGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTATTATCCAGCAGa                                                                          <  2:702829/150‑1 (MQ=255)
                                                                                   aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCTTTACTATCCAGCAGACCCGCTGc                                                                  >  1:1068960/1‑150 (MQ=255)
                                                                                     aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                                <  2:916231/150‑1 (MQ=255)
                                                                                     aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                                >  1:1043445/1‑150 (MQ=255)
                                                                                             tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGa                                                        >  2:928973/1‑150 (MQ=255)
                                                                                             tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGa                                                        >  2:928993/1‑150 (MQ=255)
                                                                                                cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGCCGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTTGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGcc                                                     <  2:1288470/150‑1 (MQ=255)
                                                                                                cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGTGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGcc                                                     <  2:1288469/150‑1 (MQ=255)
                                                                                                        tgatgaAAAACCCAAAACCCCCACCGGCCGGCCCCCCCGCCCCGTGGGGG‑TGAGGTTTCGGGTGTGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCggg                                             <  1:1375867/150‑1 (MQ=255)
                                                                                                        tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCggg                                             >  2:896808/1‑150 (MQ=255)
                                                                                                                  ccacaaaaCCGCCCCGGTCGGGCCCCCCGCCCCGTTGGGGTTTTTGTTTTGGTTGCGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACgg                                    <  2:13078/146‑1 (MQ=255)
                                                                                                                     aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                >  2:83109/1‑150 (MQ=255)
                                                                                                                     aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                >  2:1021452/1‑150 (MQ=255)
                                                                                                                             cACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGAt                        >  1:642300/1‑150 (MQ=255)
                                                                                                                                    cctaacccTCTCCCCAAAGGGG‑CGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATACAGCAGACCCGCAGCAACCGCAGGAGCCATGCCGGGCGCTAACGGGAAACGCGGGATCACCAgg                 >  1:926842/6‑150 (MQ=255)
                                                                                                                                     gggccccTCGCCCTTTGGGGG‑TGTGGGTTCGGGAGTGGGTAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGc                <  1:866445/147‑1 (MQ=255)
                                                                                                                                     ggTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGc                >  2:1078606/1‑150 (MQ=255)
                                                                                                                                      gTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCa               >  1:232743/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCGCCCGTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCaa              <  2:341379/150‑1 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCaa              >  1:1230563/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCaa              >  1:529048/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCaa              >  1:222847/1‑150 (MQ=255)
                                                                                                                                        acccTCTCCCCAAAGGGG‑CGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAAt             >  2:563397/2‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgca          >  2:311583/1‑150 (MQ=255)
                                                                                                                                                   tttGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGa  >  1:1164030/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGA  >  NC_000913/4070315‑4070611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑