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.97.4% 154.8 / ‑2.1 38coding (318/327 nt)ynfAUPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (18/19);  total (19/19)
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.13e-01

GGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAA  >  NC_000913/1655228‑1655488
                                                                                                                                |                                                                                                                                    
ggAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGc                                                                                                                           <  2:207897/139‑1 (MQ=255)
  aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                         <  1:560705/139‑1 (MQ=255)
       tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                                    >  2:484755/1‑124 (MQ=255)
       tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                                    <  1:484755/124‑1 (MQ=255)
              aaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                         <  2:516100/127‑1 (MQ=255)
              aaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                         >  1:516100/1‑127 (MQ=255)
                  ggCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGcc                                                                                                         >  2:335508/1‑139 (MQ=255)
                             gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTGTCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCaa                                                                                              >  2:477227/1‑139 (MQ=255)
                                   tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                                                        >  2:9276/1‑139 (MQ=255)
                                   tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                                                        <  2:738464/139‑1 (MQ=255)
                                        ggTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAccc                                                                                   >  1:107538/1‑139 (MQ=255)
                                           tGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                                <  1:9276/139‑1 (MQ=255)
                                             aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                         <  2:655629/96‑1 (MQ=255)
                                             aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                         >  1:655629/1‑96 (MQ=255)
                                             aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                              >  1:363443/1‑139 (MQ=255)
                                             aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                              >  1:130122/1‑139 (MQ=255)
                                                    aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                                     >  1:487834/1‑125 (MQ=255)
                                                    aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                                     <  2:487834/125‑1 (MQ=255)
                                                       aaaTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                               <  1:445688/81‑5 (MQ=255)
                                                       aaaTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                               >  2:445688/1‑77 (MQ=255)
                                                                  cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGTCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCAc                                                         <  2:384742/139‑1 (MQ=255)
                                                                    tttGTCAGCTGCCACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTTCGCCCCCGCCCGCAACAATGATCAACATGCCGCAATGCGCAATCAACGCACCCGTCCAGTCATAAAGAGGCAGTTTCACGc                                                       >  1:308801/1‑139 (MQ=255)
                                                                    tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCGAAGCGCAATCAACGCACCCg                                                                                  <  2:363189/112‑1 (MQ=255)
                                                                    tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCGAAGCGCAATCAACGCACCCg                                                                                  >  1:363189/1‑112 (MQ=255)
                                                                                               tatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg                            >  2:260222/1‑139 (MQ=255)
                                                                                               tatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg                            <  1:477227/139‑1 (MQ=255)
                                                                                                 ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                                       <  1:652638/46‑1 (MQ=255)
                                                                                                 ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                                       >  2:652638/1‑46 (MQ=255)
                                                                                                 ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                                      <  2:254398/127‑1 (MQ=255)
                                                                                                 ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                                      >  1:254398/1‑127 (MQ=255)
                                                                                                 ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGcc                          >  2:599732/1‑139 (MQ=255)
                                                                                                    aTAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg                       <  2:130122/139‑1 (MQ=255)
                                                                                                    aTAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg                       <  1:335508/139‑1 (MQ=255)
                                                                                                      aaaaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa                     >  2:420909/1‑139 (MQ=255)
                                                                                                        aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACCATGCGCAGCCACATCACCACCGTg                       <  2:490916/135‑1 (MQ=255)
                                                                                                        aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAATGCGCAGCCACATCAACACCGTg                       >  1:490916/1‑135 (MQ=255)
                                                                                                              tttATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                                  >  2:642540/1‑38 (MQ=255)
                                                                                                              tttATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                                  <  1:642540/38‑1 (MQ=255)
                                                                                                                      cttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                                                              <  2:105831/80‑1 (MQ=255)
                                                                                                                      cttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                                                              >  1:105831/3‑82 (MQ=255)
                                                                                                                       ttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCAt    <  1:554878/138‑1 (MQ=255)
                                                                                                                         aCGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATaa  <  2:107538/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                    
GGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAA  >  NC_000913/1655228‑1655488

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