Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,249,251 A→G 19.1% N140D (AAT→GAT)  yeiL → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,249,2510AG19.1% 27.6 / 4.6 21N140D (AAT→GAT) yeiLputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (7/10);  new base G (2/2);  total (9/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 = 1.80e-01

CCTTGCGCTCCCTATGAAACATTACCGTCCGCTGTTATTAAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAA  >  CP009273/2249112‑2249398
                                                                                                                                           |                                                                                                                                                   
ccTTGCGCTCCCTATGAAACATTACCGTCCGCTGTTATTAAACGACACACTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGGATCGCCTgg                                                                                                                                            >  1:206748/1‑149 (MQ=255)
ccTTGCGCTCCCTATGAAACATTACCGTCCGCTGTTATTAAACAACACGCTATATTTAAGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGGATCGCCTgg                                                                                                                                            >  1:206813/1‑149 (MQ=255)
     ctcTCCCTCTGCAACATTACCGTCCGCTGTTATTAATCAACACGGTTTGTGAACGGAGTATGGGCATCAACTGGATTCATAAGCATAACCATAATAGTGGTAATTAAACTCTGAATCAATTATTTCCGTTAGAGAATCGCCTGgcggca                                                                                                                                       <  2:158789/147‑1 (MQ=255)
        tCCCTATGAAACATTACCGTCCGCTGTTATTAAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCAttt                                                                                                                                    <  2:61083/149‑1 (MQ=255)
          ccTATGAAACATTACCGTCCGCTGTTATTAAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATAATTTCCGTTAGTGAATCGCCTGGCAGCATTtat                                                                                                                                  >  1:204499/1‑149 (MQ=255)
                      tACCGTCCGCTGTTATTAAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCa                                                                                                                      <  2:311647/149‑1 (MQ=255)
                                    attaAACGGCACGCTATTTTGTCGAACACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCttt                                                                                                        <  2:220763/149‑1 (MQ=255)
                                     ttaAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTa                                                                                                       >  1:206345/1‑149 (MQ=255)
                                        aaCGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCa                                                                                                    >  1:78102/1‑149 (MQ=255)
                                                     ttttACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTACTATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAATGTGATATTTATCACGgaaacaaaacg                                                                                       >  2:261676/1‑142 (MQ=255)
                                                        tACGAAAACTCTGCGTCACCTTAAGTCATAAATATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCaa                                                                                    >  1:163601/1‑149 (MQ=255)
                                                               cctctGCGTCACCTTAAATCATAAAAATTATCTTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCTCTTGGCATCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCag                                                                             <  1:152058/148‑1 (MQ=255)
                                                                ctctGCGTCACCTTACGTCATAAAAATTATCGTAATATTGTTTCTTTTACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGTATTTATATTACTCTCGCAGGAAGGTGAGCTTTATCACGAAAAGCATACGCAAGCGGCaga                                                                            <  1:28981/149‑1 (MQ=255)
                                                                ctctGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCaga                                                                            <  1:141659/149‑1 (MQ=255)
                                                                    gCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTACGTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGGGATCTTTATCACGAGAAGCATCCGCAAGCGCCAGATTAt                                                                        >  2:29028/1‑149 (MQ=255)
                                                                                 tCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATAATTTCCGGCAGAAAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCtt                                                           >  1:317349/1‑149 (MQ=255)
                                                                                       aaaTTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTTATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGAc                                                     <  1:108421/149‑1 (MQ=255)
                                                                                            aTCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGGATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCtt                                                <  2:206813/149‑1 (MQ=255)
                                                                                            aTCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGGATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCtt                                                <  2:206748/149‑1 (MQ=255)
                                                                                             tcggtttatTGTTTCTTCAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTTCTCTCGCAGTAAGGTGATCTTTATCACGAAAAGTATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCttt                                               <  1:136015/143‑1 (MQ=255)
                                                                                                                                          gAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTaaaa  <  2:274650/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CCTTGCGCTCCCTATGAAACATTACCGTCCGCTGTTATTAAACGACACGCTATTTTTACGAAAACTCTGCGTCACCTTAAGTCATAAAAATTATCGTAATATTGTTTCTTTAACTCAGAATCAATCATTTCCGTTAGTGAATCGCCTGGCAGCATTTATATTACTCTCGCAGGAAGGTGATCTTTATCACGAAAAGCATACGCAAGCGGCAGAGTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAA  >  CP009273/2249112‑2249398

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