Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,170,443 G→T 18.2% P60Q (CCG→CAG)  ycfT ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,170,4430GT18.2% 50.4 / 4.5 22P60Q (CCG→CAG) ycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (8/10);  new base T (2/2);  total (10/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 = 2.74e-01

TCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACATAACCCTTTGATCTGGTTAATCCATAGCTCTTTTTGTTTCATCAGGACAACACGC  >  CP009273/1170213‑1170634
                                                                                                                                                                                                                                      |                                                                                                                                                                                               
tcCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAACGCGCCAGCACCCAGCAGATGTTCCAGCTCCGTTTATCGAGCCAATTTCCCCACGGTACGCTTTCTATATCGCGGCGCATCAAATAGCCGCAAATAAAGAAAAAAACCGGCATACAAACGGGGGAAAGGAAAAGACTGAAATacaacaataa                                                                                                                                                      >  2:110637/1‑265 (MQ=255)
                    gTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGTCGAATCAAATAGCCGGATATAAAGAAAAAAACCGCCATAAGAAAGGGGGAAAGGTAAAGATTGAAATAGATCCAGCATTTGCTAAGGANTGCCGACaa                                                                                                                                 >  2:152748/1‑275 (MQ=255)
                                                          aTCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTa                                                                                                                                                                           >  2:11342/1‑195 (MQ=255)
                                                              cgcTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGa                                                                                                                                                                       <  1:11342/195‑1 (MQ=255)
                                                                                            cTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGc                                                                                                                                                                                               <  2:337680/141‑1 (MQ=255)
                                                                                            cTTAAAGCCAGCCACTGCACCACGCCCCAAAGCGACAGCACCCAGAAGATGTTCCAGATCCGTTGATCGAGACAATTTACCCACGGAACGCTGTCGATATAGAGGCGAAACAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACa                                                         <  1:110641/275‑1 (MQ=255)
                                                                                                        caCTGCACCACGACCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACATAACCCTTTGAt                                             <  1:152753/275‑1 (MQ=255)
                                                                                                           tGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCTGCATACGAAAgg                                                                                                                                                                                     >  1:6519/1‑136 (MQ=255)
                                                                                                               ccacGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCTGCATACGAAAGGGGGc                                                                                                                                                                                 <  2:6519/136‑1 (MQ=255)
                                                                                                                cacGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTg                                                                                                                        >  1:364876/1‑192 (MQ=255)
                                                                                                                    ccccAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGaaaa                                                                                                                    <  2:364870/192‑1 (MQ=255)
                                                                                                                                                                  gaCAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAgg                                                                                                                                                                                     >  2:231802/1‑81 (MQ=255)
                                                                                                                                                                      aTTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGc                                                                                                                                                                                 <  1:231807/81‑1 (MQ=255)
                                                                                                                                                                                      cTGTCGATATAGCGGCGAATCAAATAGCCNGATATAAANAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCATGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACACATACATACCCCTTTGTTCTGGTTAATCCATAGCTCTTTTTGt                    <  1:274079/222‑1 (MQ=255)
                                                                                                                                                                                              ataGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTa                                                                                                                                                                           >  2:128306/1‑63 (MQ=255)
                                                                                                                                                                                                  cggcgAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGa                                                                                                                                                                       <  1:128310/63‑1 (MQ=255)
                                                                                                                                                                                                     cgAATCAAATAGCCGGATATAAAGAAAAAAACCTGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAg                                                                                                                   >  1:74616/1‑112 (MQ=255)
                                                                                                                                                                                                         tCAAATAGCCGGATATAAAGAAAAAAACCTGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGt                                                                                                               <  2:74613/112‑1 (MQ=255)
                                                                                                                                                                                                         tCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTa                                                                                                                                                                           >  1:279560/1‑52 (MQ=255)
                                                                                                                                                                                                             aTAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGa                                                                                                                                                                       <  2:279555/52‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaaCCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACATAACCCTTTGATCTGGTTAATCCATAGCTCTTTTTGTTTCATCAGGacaac      >  1:112018/1‑195 (MQ=255)
                                                                                                                                                                                                                                   aCCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACATAACCCTTTGATCTGGTTAATCCATAGCTCTTTTTGTTTCATCAGGACAACACGc  <  2:112014/195‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                               
TCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCCAAAGCGCCAGCACCCAGAAGATGTTCCAGATCCGTTTATCGAGACAATTTCCCCACGGTACGCTGTCGATATAGCGGCGAATCAAATAGCCGGATATAAAGAAAAAAACCGGCATACGAAAGGGGGCAAGGTAAAGATTGAAATAGATCCAGCATTTGCTCAGGACTTCCGATAACGGATGCTGAAAAGTGGTCAGATGCGGATAAAAGGTAATGACCGAGTGATAAATCACCACCAGACAAATACATAACCCTTTGATCTGGTTAATCCATAGCTCTTTTTGTTTCATCAGGACAACACGC  >  CP009273/1170213‑1170634

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