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.92.9% 99.0 / 2.5 28coding (318/327 nt)ynfAUPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (14/12);  total (15/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 = 3.95e-01

GCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAAGC  >  NC_000913/1655225‑1655490
                                                                                                                                   |                                                                                                                                      
gCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                                 <  1:169958/139‑5 (MQ=255)
     aGTAATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                           <  2:388052/139‑1 (MQ=255)
       tGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                                         <  2:605656/139‑1 (MQ=255)
                                gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCaa                                                                                                >  1:714301/1‑139 (MQ=255)
                                 cgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                                 <  2:106544/106‑5 (MQ=255)
                                 cgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                                 >  1:106544/1‑102 (MQ=255)
                                         aTGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                                       <  2:34696/139‑1 (MQ=255)
                                              tGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                                  <  2:773802/139‑1 (MQ=255)
                                                aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                                    <  2:533615/103‑1 (MQ=255)
                                                aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                                <  2:478326/139‑1 (MQ=255)
                                                aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                                >  1:720587/1‑139 (MQ=255)
                                                aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                                >  1:778233/1‑139 (MQ=255)
                                                aCAAAAGATGAAATAAAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                                    >  1:533615/1‑103 (MQ=255)
                                                      gATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATa                                                                          <  1:302473/139‑1 (MQ=255)
                                                       aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaa                                                                         >  2:513098/1‑139 (MQ=255)
                                                                 aTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                                 <  2:733266/106‑1 (MQ=255)
                                                                 aTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                                 >  1:733266/1‑106 (MQ=255)
                                                                       tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGc                                                         <  2:778233/139‑1 (MQ=255)
                                                                         tGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCa                                                       >  1:242647/1‑139 (MQ=255)
                                                                             aGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCca                                                   <  2:714301/139‑1 (MQ=255)
                                                                                        ccTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                                        <  2:720587/139‑1 (MQ=255)
                                                                                         ctctTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGcc                                       >  2:643880/1‑139 (MQ=255)
                                                                                                  tatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg                              >  1:517568/1‑139 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                                    >  1:65835/1‑36 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                                    <  2:65835/37‑2 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGcc                                                        >  1:567157/1‑112 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGcc                                                        <  2:567157/112‑1 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGc                             >  1:614887/1‑139 (MQ=255)
                                                                                                   atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGc                             >  1:421722/1‑139 (MQ=255)
                                                                                                    ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGcc                            >  1:276786/1‑139 (MQ=255)
                                                                                                               gATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATa                 <  1:227057/139‑1 (MQ=255)
                                                                                                                          ttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCAt      >  2:668880/2‑139 (MQ=255)
                                                                                                                              gTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTTCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAAGc  >  2:102164/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAAGC  >  NC_000913/1655225‑1655490

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