Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 966,112 A→T 12.2% intergenic (+35/‑17) ycaQ → / → ycaR DUF1006 family protein with C‑terminal wHTH domain/peroxide and acid resistance protein, UPF0434 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273966,1120AT12.2% 77.7 / 4.1 35intergenic (+35/‑17)ycaQ/ycaRDUF1006 family protein with C‑terminal wHTH domain/peroxide and acid resistance protein, UPF0434 family
Reads supporting (aligned to +/- strand):  ref base A (18/12);  new base T (2/2);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

GCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTGGACAACCTCGCTTTCCCGCTGCGTGAT  >  CP009273/965964‑966248
                                                                                                                                                    |                                                                                                                                        
gCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCCTAACCCGa                                                                                                                                          >  1:66144/1‑149 (MQ=255)
          aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGCTGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGGTTTATGCATTCATAACCCGATCCa                                                                                                                                      >  2:266350/1‑143 (MQ=255)
          aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCa                                                                                                                                      >  2:443952/1‑143 (MQ=255)
           tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGTTCCATCTggag                                                                                                                               >  1:263588/1‑149 (MQ=255)
           tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGTTCCATCTggag                                                                                                                               >  1:263699/1‑149 (MQ=255)
             tCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGCTGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGCTAAGATTTATGCATTCATAACCCGATCCATCt                                                                                                                                   >  2:334251/1‑143 (MQ=255)
               gCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTgg                                                                                                                                 <  2:179246/143‑1 (MQ=255)
                 tAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGTACCGGCGGGGAAATAGACCCCGTCGCATATATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTggag                                                                                                                               <  2:370989/143‑1 (MQ=255)
                     tGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGa                                                                                                                     <  1:328463/149‑1 (MQ=255)
                          ggagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATAATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATgg                                                                                                                      <  2:21113/141‑1 (MQ=255)
                            agGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCGTTCATAACCCGATCCATCTGGAGGAACTATGGag                                                                                                                    >  2:108747/1‑142 (MQ=255)
                                aaCGCGGGTGACATTAGGATGCTGCACGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCACGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACACGATCCATCTGGAGGAAGGATGGATCATCGACTGc                                                                                                          >  1:266164/1‑149 (MQ=255)
                                  cgcgGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTCTGATAAGATTTATGCATTAATAACCCGATCGATCTGGCGGAGCTATGGACCATCGt                                                                                                              >  2:100087/1‑143 (MQ=255)
                                      ggTGACATTAGGATGCCGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGaaa                                                                                                    <  1:476617/149‑1 (MQ=255)
                                                     ctgcCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCgg                                                                                     <  1:142041/149‑1 (MQ=255)
                                                      tgcCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGAGCATCGTCTGCTTGAAATCATtgcctg                                                                                          >  2:468584/1‑143 (MQ=255)
                                                             aaGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATGTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAAc                                                                             <  1:369885/149‑1 (MQ=255)
                                                              aGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACg                                                                            <  1:195384/149‑1 (MQ=255)
                                                                            aTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACtt                                                                    >  2:171119/1‑143 (MQ=255)
                                                                                  gCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGtat                                                              >  2:40888/1‑143 (MQ=255)
                                                                                       ggCTGGGTAATCGACCCCGTCGCATAAATAATTATGACAACACTTATGCATTCATACCCCGATCCATCTGCAGGAACTCTGTATCCTCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCa                                                         <  2:465934/143‑1 (MQ=255)
                                                                                           gggAAATAGACCCCGTCGCATAAATGATTATGCTAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATACCCAGGAAAAACaa                                               <  1:360233/149‑1 (MQ=255)
                                                                                                  aGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTa                                        >  1:320460/1‑149 (MQ=255)
                                                                                                    aCCCCGTCGCCTAAATGAATATGATTAGATTTATGCATTCACACCCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAagag                                            <  2:340734/143‑1 (MQ=255)
                                                                                                      cccGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCTTTGCCTGCCTGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTg                                    >  1:283358/1‑149 (MQ=255)
                                                                                                      cccGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTg                                    >  1:271285/1‑149 (MQ=255)
                                                                                                          tCGGATAGATGATTAAGATAAGAATTAAAAAATAATAACCCGTTCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAAACAATGCCTGCCCGGTTTGCAACGGAAAACTTTGAAAAAAAAAAAAAAAACAAGAGCTTAtt                                      <  2:231750/143‑1 (MQ=255)
                                                                                                          tCGCATAAATGATGATGATAATATTTATGCATTCATAACCCGTTCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTAtt                                      <  2:231845/143‑1 (MQ=255)
                                                                                                             cATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTg                             >  1:250174/1‑149 (MQ=255)
                                                                                                               tAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCaa                                 >  2:356004/1‑143 (MQ=255)
                                                                                                               tAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCATATTTGCAAACTGGa                           >  1:74537/1‑149 (MQ=255)
                                                                                                                    gATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTgg                            >  2:396881/1‑143 (MQ=255)
                                                                                                                    gATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAAATTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTgg                            >  2:49335/1‑143 (MQ=255)
                                                                                                                          gATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTGGACAAcc                      <  2:132372/143‑1 (MQ=255)
                                                                                                                                        aTTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTGTAACCAGGAAATACAAGAGCTTATTTGCAAACTGAGCAACCTAGCTTTCTCGCTGAGTAAt  >  1:67075/1‑149 (MQ=255)
                                                                                                                                           cgtAACCCGGTCCATCTGGAGGAACTATGGATAATCGTCTGCTTGAAAGCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGAAAAACAAAAAAAAACAAGAGCTTATTTGCAAACTGGACAACCTCGCTTTCCCGCTGCGt     <  2:138748/141‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
GCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTGCCTGCCCGGTTTGCAACGGAAAACTTTGGTATAACCAGGAAAAACAAGAGCTTATTTGCAAACTGGACAACCTCGCTTTCCCGCTGCGTGAT  >  CP009273/965964‑966248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: