Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,591,925 C→T 17.4% G134D (GGC→GAC)  ugpB ← sn‑glycerol‑3‑phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,591,9250CT17.4% 53.2 / 5.5 23G134D (GGC→GAC) ugpBsn‑glycerol‑3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (5/14);  new base T (2/2);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGCCGGAGGCTTTCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGTCATACACCGGTTTAATGGCTTTCGACGCCATCATGGTGGCGGTGCCAACTTCATAAAC  >  NC_000913/3591776‑3592067
                                                                                                                                                     |                                                                                                                                              
tGCCGGAGGCTTTCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTATTCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGc                                                                                                                                                <  1:69661/150‑1 (MQ=255)
     gAGGCTTTCAGTTTAGCGGCATAGTCCGCCACATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATATAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGtt                                                                                                                                           >  1:117199/1‑150 (MQ=255)
       ggCTTTGAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGtttt                                                                                                                                         <  1:175595/150‑1 (MQ=255)
       ggCTTTCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGtttt                                                                                                                                         <  1:242086/150‑1 (MQ=255)
          tttCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCt                                                                                                                                      <  2:165193/150‑1 (MQ=255)
          tttCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCt                                                                                                                                      <  2:131861/150‑1 (MQ=255)
                  tCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGagt                                                                                                                              >  2:120734/1‑150 (MQ=255)
                                  cAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGt                                                                                                              <  1:40829/150‑1 (MQ=255)
                                    gATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCg                                                                                                            <  1:57915/150‑1 (MQ=255)
                                    gATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCg                                                                                                            <  2:77822/150‑1 (MQ=255)
                                    gATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCg                                                                                                            <  2:77820/150‑1 (MQ=255)
                                     aTCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCgg                                                                                                           <  2:24608/150‑1 (MQ=255)
                                                 tttCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGa                                                                                               <  2:47717/150‑1 (MQ=255)
                                                       ctgcTGTACCTGGTCTAATCCTTCTTTCTTGAAGGCGTCTTTGTTGTAATTGAGTACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTATCCTGAAACCGTCGGCACAAACTTCGACTCATc                                                                                         <  2:208376/148‑1 (MQ=255)
                                                       cggcTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATc                                                                                         <  2:208377/150‑1 (MQ=255)
                                                          cTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGaa                                                                                      <  1:8355/150‑1 (MQ=255)
                                                                               tttCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTaaa                                                                 >  1:202484/1‑150 (MQ=255)
                                                                               tttCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTaaa                                                                 >  1:202481/1‑150 (MQ=255)
                                                                                                                               gAATGGCTGGGAGAGTAAGTGGTCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGt                                                            <  2:86952/107‑1 (MQ=255)
                                                                                                                               gAATGGCTGGGAGAGTAAGTGGTCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGt                                                            <  2:86978/107‑1 (MQ=255)
                                                                                                                               gAATGGCTGGGAGAGTAAGTGGTCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGt                                                            >  1:86978/1‑107 (MQ=255)
                                                                                                                               gAATGGCTGGGAGAGTAAGTGGTCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGt                                                            >  1:86952/1‑107 (MQ=255)
                                                                                                                                         gagaGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGTCATACACCGGTTTAATGGCTTTCGACGCCATCATGGTGGCGGTGCCAACTTCa       <  1:79162/150‑1 (MQ=255)
                                                                                                                                              tAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGTCATACACCAGTTTAATGGCTTTCGACGCCATCATGGTGGCGGTGCCAACTTCATAAAc  >  1:228855/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGCCGGAGGCTTTCAGTTTCGCGGCATAGTCCGCCAGATCCTGCCAGGTTTTCGGCGGCTGTTCCGGGTCTAATCCTGCTTTCTTGAAGGCGTCTTTGTTGTAATAGAGAACGGGGGTCGAGCTGTTGAATGGCTGGGAGAGTAAGTGGCCCGTTTTGCTGTCGGAGTAGTAACCTGAAACCGTCGGCACAAACTGCGACTCATCGAACTGAATCCCTGCCTCTTTAAACACGTCATACACCGGTTTAATGGCTTTCGACGCCATCATGGTGGCGGTGCCAACTTCATAAAC  >  NC_000913/3591776‑3592067

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