Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,263,834 G→A 27.7% intergenic (‑339/‑29) fruB ← / → setB fructose‑specific PTS multiphosphoryl transfer protein FruB/sugar exporter SetB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,263,8340GA27.7% 36.8 / 15.3 30intergenic (‑339/‑29)fruB/setBfructose‑specific PTS multiphosphoryl transfer protein FruB/sugar exporter SetB
Reads supporting (aligned to +/- strand):  ref base G (7/14);  new base A (4/4);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCG  >  NC_000913/2263709‑2263970
                                                                                                                              |                                                                                                                                        
cGTTATCAGCAATCTTATTTCAGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTTGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGt                                                                                                                         <  2:30075/143‑1 (MQ=255)
   tATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGAt                                                                                                                      <  1:172243/143‑1 (MQ=255)
   tATCAGCAATTTCATTTCAGCAGCGGCAGCAGCCAAGCAGCAACAGTGAAATGGCATGTTGACTATTCGTGGCGCTGTTGAGTTCGCTGGCTGTTTATTTTCTGACAAGCAGCGTGAACTCCGCGTCTTCCTCTTCCAGTGAt                                                                                                                      <  2:26024/143‑1 (MQ=255)
          aaTTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAg                                                                                                               <  2:429913/143‑1 (MQ=255)
           gttttatttcaGCAGCGTCAGCAGCCAATCAGCAACAGTGGAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAgc                                                                                                              <  2:517668/142‑1 (MQ=255)
            tttcattttaGTAGCGTAAGCAGCCAATGAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGTAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAgca                                                                                                             <  1:23683/143‑1 (MQ=255)
            tttcatttcaGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAgca                                                                                                             >  1:497360/1‑143 (MQ=255)
                 tttcaGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAgcatgcat                                                                                                        <  1:309818/143‑1 (MQ=255)
                 tttcaGCAGCGGGAGGAGCGAAGCGGCAACAGTGGAATGGCATGTTGACGATTCGTGGCGCTTTTGATTTCGCTGTCGGTTTATTTTCTGACAAGCAGCGTACCAATGCCTTCTTCCTCTTCCAGTGATCGTCCAgcatgcat                                                                                                        <  2:118254/143‑1 (MQ=255)
                              agcagcCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAg                                                                                           >  1:392549/1‑143 (MQ=255)
                              agcagcCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAg                                                                                           >  1:390781/1‑143 (MQ=255)
                              agcagcCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAg                                                                                           >  1:390778/1‑143 (MQ=255)
                              agcagcCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAg                                                                                           >  1:390774/1‑143 (MQ=255)
                                   ccAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTcc                                                                                      >  1:195529/1‑143 (MQ=255)
                                       tCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAgcg                                                                                  <  1:531167/143‑1 (MQ=255)
                                       tCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAgcg                                                                                  <  2:289251/143‑1 (MQ=255)
                                               aGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCg                                                                          <  1:339100/143‑1 (MQ=255)
                                                gggAAGTGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGt                                                                         <  2:371596/141‑1 (MQ=255)
                                                 tGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGGCTGCGGCGCGAACTAGtt                                                                        >  1:72248/1‑143 (MQ=255)
                                                      tGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCCCAGTCCCCCGCCCGCAATCGttccacc                                                                   >  2:535632/1‑138 (MQ=255)
                                                              tGACGATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCg                                                           <  1:518815/143‑1 (MQ=255)
                                                                aCTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGAc                                                         <  2:390774/143‑1 (MQ=255)
                                                                aCTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGAc                                                         <  2:390778/143‑1 (MQ=255)
                                                                aCTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGAc                                                         <  2:390781/143‑1 (MQ=255)
                                                                aCTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCACGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGAc                                                         <  2:392549/143‑1 (MQ=255)
                                                                            gcTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAAt                                             <  1:239602/143‑1 (MQ=255)
                                                                                ttGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACATGACCTCGACGGCGTGTTTAATCGtt                                         >  1:500653/1‑143 (MQ=255)
                                                                                          tgtctgtTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                               <  2:497360/143‑1 (MQ=255)
                                                                                                        cTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTCTCTCACCGATATGGCGGgc                 >  2:2235/1‑143 (MQ=255)
                                                                                                                        aaCTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTCTTAATCGTTGCCTCTCTCACCGGTATGACGGGCGCTCCGCCAACCCCg  >  1:471089/1‑143 (MQ=255)
                                                                                                                              |                                                                                                                                        
CGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACAGTGAAATGGCATTTTGACTATTCGTGGCGCTTTTGATTTCGCTGTCTGTTTATTTTCTGACAAGCAGCGTAAACTCCGCGTCTTCCTCTTCCAGTGATCGACCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCG  >  NC_000913/2263709‑2263970

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