Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,796,782 G→A T268T (ACG→ACA waaC → ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,7820GA92.3% 32.9 / ‑5.7 13T268T (ACG→ACAwaaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (4/8);  total (4/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGA  >  NC_000913/3796539‑3796881
                                                                                                                                                                                                                                                   |                                                                                                    
gacCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCAc                                                                                               <  1:388341/251‑1 (MQ=255)
          gatAAACACTGGTCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAAGACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTAGTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGCGGCTAAATTGGTAGTGTCAGTGGATACGGGGTTAAGCCATTGATCAGCGGCACGGTATAGAc                                                                                      <  2:51684/250‑1 (MQ=255)
             aaaCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCa                                                                                   >  1:79314/1‑250 (MQ=255)
                  cTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATc                                                                              >  2:283719/1‑250 (MQ=255)
                    ggCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACg                                                                           <  1:283797/251‑1 (MQ=255)
                                                         ggCTGAGTTAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGGCTGGCGGAAGGATTTGCTTATGTTGAAGTATGGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATtggtgg                                       <  1:288267/249‑1 (MQ=255)
                                                              aTTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATggg                                >  1:109764/1‑251 (MQ=255)
                                                                tCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGga                               <  2:109764/250‑1 (MQ=255)
                                                                         cGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGAt                     <  1:283719/251‑1 (MQ=255)
                                                                          ggATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGAGTTGCTTATGTTGAAGTAGTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAAGTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTGTATGGACCTACCGATCCGGGATTAATTGGGGGTTATGGGAAGAATCAGATg                    <  2:174758/251‑1 (MQ=255)
                                                                             ttAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGtt                 >  1:280539/1‑250 (MQ=255)
                                                                                aaCTTCCGTGGGGCGCGCCGCATGAGGAAGAACTGGCGATACGACTGGCGGAAGGATTTGCGTATGTTGAAGTATTGGCTAAGATGAGTCTGGAAGTCGTTGCTCGCGTGCTGTCCGGGGCTAAATTTGGAGTTTCAGTGGATACGGGGTTAGGCCCTTTAACGGCGGCAATGGATGGTGTCAATATCTCGGTTTATTGTACTAGGATTCCGGGATTAATTGATGGGTATGGGAAGAATCAGATGGTTTGt              <  1:175282/251‑1 (MQ=255)
                                                                                            gcgcgaCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTAGGGCTccgggg  <  2:149413/251‑5 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                    
GACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGA  >  NC_000913/3796539‑3796881

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