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.2% 138.2 / ‑3.2 36coding (318/327 nt)ynfAUPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (20/15);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAG  >  NC_000913/1655220‑1655470
                                                                                                                                        |                                                                                                                  
aaaaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                  <  1:1244710/139‑1 (MQ=255)
aaaaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                  >  2:143274/1‑139 (MQ=255)
  aaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                <  1:438361/139‑2 (MQ=255)
                        aaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATg                                                                                         >  1:1150780/1‑139 (MQ=255)
                          ggCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                <  2:317654/115‑2 (MQ=255)
                          ggCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                >  1:317654/1‑114 (MQ=255)
                          ggCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGcc                                                                                       >  1:635177/1‑139 (MQ=255)
                                     gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                >  2:244024/1‑119 (MQ=255)
                                     gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGAt                                                                                                <  1:244024/119‑1 (MQ=255)
                                                   tgtcAAAACATGAAAGAGAATACCTTTTGTCAGCTGACTATTCCTCTTATCTTATAGATAAAACGGATTTATGTTCTTACGTGCCCCCCAGCCCGCAACAATGACCAACCTGCGGCAAAGCGCAATCACCGCACCCGTc                                                              <  2:1150780/136‑1 (MQ=255)
                                                     aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                            >  2:98755/1‑139 (MQ=255)
                                                     aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                            >  1:860173/1‑139 (MQ=255)
                                                                 aTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGt                                                <  2:635177/139‑1 (MQ=255)
                                                                 aTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGt                                                <  1:143274/139‑1 (MQ=255)
                                                                       tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATc                                                                                               <  1:1064422/86‑1 (MQ=255)
                                                                       tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATc                                                                                               >  2:1064422/1‑86 (MQ=255)
                                                                       tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGttt                                          >  1:1253175/1‑139 (MQ=255)
                                                                       tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGttt                                          >  2:1215218/1‑139 (MQ=255)
                                                                          cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCAc                                       <  2:860173/139‑1 (MQ=255)
                                                                               gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                              <  2:13031/111‑1 (MQ=255)
                                                                               gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                              >  1:13031/1‑111 (MQ=255)
                                                                               gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                            <  1:86392/113‑1 (MQ=255)
                                                                               gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                            >  2:86392/1‑113 (MQ=255)
                                                                                        cacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgc                         <  2:1253175/139‑1 (MQ=255)
                                                                                             ccTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGCTTCACGCCATCCACAACGCGCAGc                    <  1:921506/139‑1 (MQ=255)
                                                                                              ctctTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGcc                   >  2:856466/1‑139 (MQ=255)
                                                                                                     cttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                     >  1:696870/1‑50 (MQ=255)
                                                                                                     cttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                     <  2:696870/50‑1 (MQ=255)
                                                                                                        atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAg                                                                                 <  1:983327/67‑1 (MQ=255)
                                                                                                        atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAg                                                                                 >  2:983327/1‑67 (MQ=255)
                                                                                                        atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGc         >  1:301854/1‑139 (MQ=255)
                                                                                                         ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCg                                                                                                          <  1:582797/41‑1 (MQ=255)
                                                                                                         ttGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCg                                                                                                          >  2:582797/1‑41 (MQ=255)
                                                                                                          tGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                                    <  1:510780/78‑1 (MQ=255)
                                                                                                          tGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                                    >  2:510780/1‑78 (MQ=255)
                                                                                                           gATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCa                                                                                    >  1:1492344/1‑61 (MQ=255)
                                                                                                           gATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCa                                                                                    <  2:1492344/61‑1 (MQ=255)
                                                                                                              aaaaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa   >  1:123713/1‑139 (MQ=255)
                                                                                                               aaaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAg  >  1:581985/1‑139 (MQ=255)
                                                                                                                              cttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                              <  2:1463003/62‑1 (MQ=39)
                                                                                                                              cttACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                              >  1:1463003/3‑64 (MQ=39)
                                                                                                                                        |                                                                                                                  
AAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAG  >  NC_000913/1655220‑1655470

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: