Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,299,142 G→A 12.8% G182D (GGC→GAC)  ychE → UPF0056 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,299,1420GA12.8% 74.7 / 2.2 32G182D (GGC→GAC) ychEUPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base A (2/2);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TTCTACCATCGTCTGGGGTACGCGTTATCACAGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGT  >  NC_000913/1299002‑1299289
                                                                                                                                            |                                                                                                                                                   
ttCTACCATCGTCTGGGGTACGCGTTATCACAGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACg                                                                                                                                             >  2:139535/1‑149 (MQ=255)
                   aCGCGTTATTGCAGCATTCGCTATCTGTTTGGTTTTTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATggg                                                                                                                          <  2:200217/149‑1 (MQ=255)
                     gcgTTATCACAGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGc                                                                                                                        >  1:129246/1‑149 (MQ=255)
                      cgTTATCACAGCGTTAGATATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTGTATGTTGTTGGGGATTGTTCCGCGTTGCACCGTGGCTGGTACGGGGTTTACGCCATACCGCCTTCAACGTGCTTCCGCGTATCATGGAGCt                                                                                                                       >  2:207662/1‑149 (MQ=255)
                      cgTTATCACAGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGGGTATTATGGGGCt                                                                                                                       >  2:15436/1‑149 (MQ=255)
                      cgTTATCACAGCATTAGCTATCTGTATGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTGTTATGGGGCt                                                                                                                       >  1:404551/1‑149 (MQ=255)
                             acaGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGACATCAACGTGATTACGCGTATTATGGGGCTATTGCTg                                                                                                                >  2:311230/1‑149 (MQ=255)
                             acaGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGACATCAACGTGATTACGCGTATTATGGGGCTATTGCTg                                                                                                                >  2:311198/1‑149 (MQ=255)
                                      gCTATCTGTTTGTTTTCTTTGTGGCTATTGCATTTTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTATACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGACGGCATTg                                                                                                       <  2:577046/149‑1 (MQ=255)
                                         aTCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTgggg                                                                                                    >  2:773722/1‑149 (MQ=255)
                                             gTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTg                                                                                                >  1:6976/1‑149 (MQ=255)
                                                           tgGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTAATCGTTACt                                                                                  >  2:282398/1‑149 (MQ=255)
                                                             gCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTgg                                                                                >  2:773355/1‑149 (MQ=255)
                                                                 ttGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTAtt                                                                            <  1:47503/149‑1 (MQ=255)
                                                                      ttGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGACATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAggg                                                                       <  1:311198/149‑1 (MQ=255)
                                                                      ttGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGACATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAggg                                                                       <  1:311230/149‑1 (MQ=255)
                                                                            gCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGGTTACGCGTATTATGGGGCTATTGCTGATGGCATTTGGGATTGAAATTATCGTTACTGGTCTTAAGGGGCtttt                                                                 >  2:184460/1‑149 (MQ=255)
                                                                            gCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGATTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGGTGGCATTGGGCGTTGAAATTATCATTCCTGTTATTTAGGGGAtttt                                                                 >  2:184404/1‑149 (MQ=255)
                                                                                tATGTTGTTGTGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTcccc                                                             <  1:458209/149‑1 (MQ=255)
                                                                                               tGTTCCGCATGGCACCGTGGCTGGTACTGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCttaattaa                                              <  1:619571/149‑1 (MQ=255)
                                                                                                  tCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTGTCGTTACTGGTTTTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTc                                           >  2:348302/1‑149 (MQ=255)
                                                                                                  tCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTGCCCGGCCTGCTTAATTAATTc                                           >  2:348591/1‑149 (MQ=255)
                                                                                                      cATGCCACCGTGGCTGGTACGGGTTTTCCGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTCTCGTTCCTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCttt                                       <  1:751696/149‑1 (MQ=255)
                                                                                                       aTGGCACCTGGGCTGGTACGGGTTTTACGCCAGACCGGCATCATCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTc                                      <  2:405514/149‑1 (MQ=255)
                                                                                                       aTGACACCGTGGCTGGGACGGGTTTTCCGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTc                                      <  2:19470/149‑1 (MQ=255)
                                                                                                                              tttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCt               <  2:157834/149‑1 (MQ=255)
                                                                                                                               ttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTc              <  1:324597/149‑1 (MQ=255)
                                                                                                                                tACGCCAGGCCGGCATCAACGTGATTACGCGTATTATGGGGCTATGGCTGGTGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCAAGCCTACTTAATTAGTTCCTTTCAAATGAAACGGAGCTGCCATGCTcc             >  2:357965/1‑149 (MQ=255)
                                                                                                                                tACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGCTTTTCACCGGCCTGGGAAATATAGTTCTTTCAAATGAAACGGAGCTGCCATGCTcc             >  2:357911/1‑149 (MQ=255)
                                                                                                                                  cGCAAGACCGGCAGCAACGTGATTACGCTTATTATGGGGCTATTGCTGATGGCATTTGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGTCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGt           <  1:140409/149‑1 (MQ=255)
                                                                                                                                    ccAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTCATCGCTGCCTGCTTAATTATTGCCTTTCAAATGAAACGGAGCTGCCATGCTCCGttt         >  1:569524/1‑149 (MQ=255)
                                                                                                                                        aCCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACtt     >  1:593260/1‑149 (MQ=255)
                                                                                                                                           ggCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGt  >  1:634079/1‑149 (MQ=255)
                                                                                                                                           ggCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGt  <  1:139535/149‑1 (MQ=255)
                                                                                                                                           ggCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGAATGATGTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGt  >  1:737932/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TTCTACCATCGTCTGGGGTACGCGTTATCACAGCATTAGCTATCTGTTTGGTTTCTTTGTGGCTATTGCATTGTTCGCTTTATGTTGTTGGGGATTGTTCCGCATGGCACCGTGGCTGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGT  >  NC_000913/1299002‑1299289

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: ‑