Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,131,524 A→G 18.6% S295S (AGT→AGC) ‡ yghQ ← putative inner membrane polysaccharide flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,131,5240AG18.6% 95.3 / 18.2 43S295S (AGT→AGC) ‡yghQputative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base G (4/4);  total (22/21)
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.67e-02

GCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCA  >  NC_000913/3131385‑3131621
                                                                                                                                           |                                                                                                 
gCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                           >  4:448235/1‑148 (MQ=255)
                       aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                    >  2:550059/1‑148 (MQ=255)
                       aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                          <  1:110599/126‑1 (MQ=255)
                       aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                          >  2:110599/1‑126 (MQ=255)
                       aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                          >  2:522911/1‑126 (MQ=255)
                       aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                          <  1:522911/126‑1 (MQ=255)
                               cgcgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             <  2:504468/115‑1 (MQ=255)
                               cgcgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             >  1:504468/1‑115 (MQ=255)
                                    gcgcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                       <  3:448235/148‑1 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                              >  4:365927/1‑103 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                              <  3:365927/103‑1 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                              <  1:32900/103‑1 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                              >  2:32900/1‑103 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             >  3:628024/1‑104 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             <  4:628024/104‑1 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                           <  1:124169/106‑1 (MQ=255)
                                          ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                           >  2:124169/1‑106 (MQ=255)
                                                    gATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             >  3:979352/1‑94 (MQ=255)
                                                    gATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                             <  4:979352/94‑1 (MQ=255)
                                                             cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                       <  1:873866/91‑2 (MQ=255)
                                                             cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                       >  2:873866/1‑90 (MQ=255)
                                                                   tAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTACCGGCGGCGTCgaa                                                          <  2:513107/114‑1 (MQ=255)
                                                                   tAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGAGCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                          >  1:513107/1‑114 (MQ=255)
                                                                      aCCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                          <  2:515109/111‑1 (MQ=255)
                                                                      aCCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                          >  1:515109/1‑111 (MQ=255)
                                                                              tCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGt                                                    <  1:337991/109‑1 (MQ=255)
                                                                              tCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGt                                                    >  2:337991/1‑109 (MQ=255)
                                                                                 cacCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                    >  4:1074093/1‑90 (MQ=255)
                                                                                 cacCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                    <  3:1074093/90‑1 (MQ=255)
                                                                                     tAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                       >  2:478280/1‑66 (MQ=38)
                                                                                     tAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                       <  1:478280/67‑2 (MQ=38)
                                                                                        caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaaga                                                        >  1:1078410/1‑95 (MQ=255)
                                                                                        caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaaga                                                        <  2:1078410/95‑1 (MQ=255)
                                                                                        caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaa                                      <  1:53627/113‑1 (MQ=255)
                                                                                        caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaa                                      >  2:53627/1‑113 (MQ=255)
                                                                                             aCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                           <  3:928439/103‑1 (MQ=255)
                                                                                             aCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                           >  4:928439/1‑103 (MQ=255)
                                                                                                                      aCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCcca  >  3:72878/1‑119 (MQ=255)
                                                                                                                      aCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCcca  <  4:72878/119‑1 (MQ=255)
                                                                                                                              ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                >  1:469583/1‑97 (MQ=255)
                                                                                                                              ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                <  2:469583/97‑1 (MQ=255)
                                                                                                                                tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                   <  1:403268/92‑1 (MQ=255)
                                                                                                                                tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                   >  2:403268/1‑92 (MQ=255)
                                                                                                                                           |                                                                                                 
GCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCA  >  NC_000913/3131385‑3131621

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