Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,655,360 (C)5→4 coding (314/327 nt) ynfA ← UPF0060 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,655,3560C.96.4% 106.9 / ‑2.3 28coding (318/327 nt)ynfAUPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (15/12);  minor base G (1/0);  total (16/12)
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.76e-01

TGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACG  >  NC_000913/1655256‑1655479
                                                                                                    |                                                                                                                           
tGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                      <  1:184631/139‑1 (MQ=255)
       tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                                               >  1:512647/1‑139 (MQ=255)
                 aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                     >  2:486604/1‑139 (MQ=255)
                 aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                     >  2:382154/1‑139 (MQ=255)
                                gAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                      >  2:654761/1‑107 (MQ=255)
                                gAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                      <  1:654761/107‑1 (MQ=255)
                                 aaTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                                                     >  2:591360/1‑139 (MQ=255)
                                                   acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaaca                                                                                                              >  1:780387/1‑64 (MQ=255)
                                                   acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaaca                                                                                                              <  2:780387/64‑1 (MQ=255)
                                                        tCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAg                              >  1:567283/1‑139 (MQ=255)
                                                         ccTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                             <  2:567283/139‑1 (MQ=255)
                                                                   tatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgc                                  <  1:1114400/124‑1 (MQ=255)
                                                                   tatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgc                                  >  2:1114400/1‑124 (MQ=255)
                                                                   tatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg                   <  1:486604/139‑1 (MQ=255)
                                                                    atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  <  1:204597/107‑1 (MQ=255)
                                                                    atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  >  2:204597/1‑107 (MQ=255)
                                                                     ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                            >  1:288447/1‑80 (MQ=255)
                                                                     ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                            <  2:288447/80‑1 (MQ=255)
                                                                     ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  >  1:405472/1‑106 (MQ=255)
                                                                     ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  <  2:405472/106‑1 (MQ=255)
                                                                          aaaaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa            >  1:1713/1‑139 (MQ=255)
                                                                          aaaaTGGATTTATGTTCTTACGTG‑GGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa            >  1:693489/1‑139 (MQ=255)
                                                                                  tttATGTTCTCACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATaaa    >  1:1074086/1‑139 (MQ=255)
                                                                                    tATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACg  <  2:693489/139‑1 (MQ=255)
                                                                                             aCGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                       <  2:164972/61‑1 (MQ=255)
                                                                                             aCGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                       >  1:164972/1‑61 (MQ=255)
                                                                                             aCGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGc             <  2:1024011/119‑1 (MQ=255)
                                                                                             aCGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGc             >  1:1024011/1‑119 (MQ=255)
                                                                                                    |                                                                                                                           
TGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACG  >  NC_000913/1655256‑1655479

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