Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,491,756 G→A 13.3% G27D (GGT→GAT)  ydcA → putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,491,7560GA13.3% 71.6 / 4.8 30G27D (GGT→GAT) ydcAputative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base A (2/2);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

ACGGTATCAAAATAAACCAGCGTATAACGTTCAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCG  >  NC_000913/1491608‑1491899
                                                                                                                                                    |                                                                                                                                               
aCGGTATCAAAATAAACCAGCGTATAACGTTCAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTgg                                                                                                                                                 <  3:316892/149‑1 (MQ=255)
        aaaaTAAACCAGCGTATAACGTTCAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAg                                                                                                                                         >  1:347848/1‑149 (MQ=255)
                  aGCGTATAACGTTCAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGAtctc                                                                                                                               <  4:329868/149‑1 (MQ=255)
                               cAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCa                                                                                                                  <  4:180851/149‑1 (MQ=255)
                                   tCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAggcg                                                                                                              <  4:450175/149‑1 (MQ=255)
                                             gATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGt                                                                                                    <  3:84204/149‑1 (MQ=255)
                                                tCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTg                                                                                                 <  1:525084/149‑1 (MQ=255)
                                                  gATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTa                                                                                               <  1:329293/149‑1 (MQ=255)
                                                   aTAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTaa                                                                                              <  3:477359/149‑1 (MQ=255)
                                                        ccATGGTGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAatgatg                                                                                         >  2:442342/1‑149 (MQ=255)
                                                                 ttATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTa                                                                                <  2:354079/149‑1 (MQ=255)
                                                                           aaaCTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGAtctc                                                                                                                               <  2:538821/92‑1 (MQ=255)
                                                                           aaaCTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGAtctc                                                                                                                               >  1:538821/1‑92 (MQ=255)
                                                                                ggCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGaaaaaaa                                                                 <  2:347848/149‑1 (MQ=255)
                                                                                    cTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCGTACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATg                                                             >  3:409459/1‑149 (MQ=255)
                                                                                                aTGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAg                                                 >  3:108129/1‑149 (MQ=255)
                                                                                                  gggAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGATTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAg                                                                               >  2:50678/1‑117 (MQ=255)
                                                                                                  gggAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGATTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAg                                                                               <  1:50678/117‑1 (MQ=255)
                                                                                                  gggAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGATTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGACAAATTTGTCTGTAATGATGGTTCTATTAg                                                                               <  3:138595/117‑1 (MQ=255)
                                                                                                  gggAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGATTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGACAAATTTGTCTGTAATGATGGTTCTATTAg                                                                               >  4:138595/1‑117 (MQ=255)
                                                                                                       cGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAGACCCTGTTCTGGTTCGAAAGGGGGGATCTCATACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTaaaa                                          <  4:409459/149‑1 (MQ=255)
                                                                                                        gCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAg                                         >  2:168943/1‑149 (MQ=255)
                                                                                                          ttGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAggg                                       <  3:66256/149‑1 (MQ=255)
                                                                                                                  ttttATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATg                               >  1:140043/1‑149 (MQ=255)
                                                                                                                         ccgatgccgGGCGCAAACCCTGTTCTGGTTGGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAAt                        <  4:163623/149‑1 (MQ=255)
                                                                                                                                cgGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACa                                                               <  2:417963/103‑1 (MQ=255)
                                                                                                                                cgGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACa                                                               >  1:417963/1‑103 (MQ=255)
                                                                                                                                   gcgcAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTa              <  1:495673/149‑1 (MQ=255)
                                                                                                                                   gcgcAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTa              <  1:168943/149‑1 (MQ=255)
                                                                                                                                               ttCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGcgcg  >  3:258419/1‑149 (MQ=255)
                                                                                                                                               ttCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGcgcg  >  2:332461/1‑149 (MQ=255)
                                                                                                                                                    gTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAgg                                        >  4:111026/1‑106 (MQ=255)
                                                                                                                                                    gTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAgg                                        <  3:111026/106‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACGGTATCAAAATAAACCAGCGTATAACGTTCAAATCGAACATAAGATTCGATAAACCATGGAGGTTATATGAAAAAACTGGCACTTATTTTGTTTATGGGAACGCTTGTTTCCTTTTATGCCGATGCCGGGCGCAAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCG  >  NC_000913/1491608‑1491899

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