Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 877,748 C→A 15.1% S148I (AGC→ATC)  deoR ← deoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273877,7480CA15.1% 51.8 / 3.5 27S148I (AGC→ATC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base A (2/2);  total (13/14)
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.32e-01

AGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGT  >  CP009273/877608‑877895
                                                                                                                                             |                                                                                                                                                   
aGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGa                                                                                                                                             >  1:87266/1‑149 (MQ=255)
      ttACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGATGGCGTGAAATTca                                                                                                                                       >  1:156952/1‑149 (MQ=255)
      ttACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGAACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGGTTGAAAATGGCGTTGATGGCGTGAAATTca                                                                                                                                       >  1:196449/1‑149 (MQ=255)
           ggCAACTCGTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTc                                                                                                                                        <  2:217769/143‑1 (MQ=255)
                 tctTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAg                                                                                                                            <  1:366536/149‑1 (MQ=255)
                 tcctCAAGATTAAGACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAg                                                                                                                            <  1:14144/146‑1 (MQ=255)
                   ttCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGaa                                                                                                                          >  1:170879/1‑149 (MQ=255)
                        gATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGACATAATCGGGTAGAAGTTATTCAGCGTTTGCTGACAATCGGTGGGTTTGAAAATGGCGGTGCTGGCGAGAGAATCACCGCCGCGTAGACAcgcg                                                                                                                     >  1:124153/1‑149 (MQ=255)
                                         ccTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTtctc                                                                                                    <  1:145539/149‑1 (MQ=255)
                                            tcaccgACATGCCCGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATCATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCttt                                                                                                 <  1:147933/144‑1 (MQ=255)
                                             tACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGACAGTGGACGATCAGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTc                                                                                                >  1:86027/1‑149 (MQ=255)
                                                tGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAgc                                                                                             >  1:289301/1‑149 (MQ=255)
                                                     tGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGCTTCTCTTTCAg                                                                                              >  2:245594/1‑143 (MQ=255)
                                                                cgcAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGGTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAgg                                                                                   >  2:16797/1‑143 (MQ=255)
                                                                 gcAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGt                                                                                  >  2:160076/1‑143 (MQ=255)
                                                                  cAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGc                                                                           >  1:58559/1‑149 (MQ=255)
                                                                  cAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGc                                                                           >  1:357947/1‑149 (MQ=255)
                                                                                       ggCTAAAATTATTCAGCGGTTGCTGAAAATCGATGGGTTTGAAAAGGGCGTTGATGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAAcggcg                                                            <  2:173295/143‑1 (MQ=255)
                                                                                       ggCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGATGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAAcggcg                                                            <  2:138532/143‑1 (MQ=255)
                                                                                        gngaaaaTTATTCATCGTTTGCCGAAACTGGATTGGTTTGAAAATGGCGTTGCTGGCGTGAAATTGACAGCCGAACAGCGCCGGCTCTCAATTGGGTTTCTCGGTCAGCGCCAGAAAGGTACTTATCGAATAACAAAcggcgg                                                           <  2:318602/140‑1 (MQ=255)
                                                                                         cAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAgg                                                    <  1:181855/149‑1 (MQ=255)
                                                                                                  ttCCGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCAtt                                           <  1:172023/149‑1 (MQ=255)
                                                                                                    cAGCGTTTGCGGAAAATCGATGGGTTGTAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGAttt                                               <  2:10188/143‑1 (MQ=255)
                                                                                                                      gATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCaa                             <  2:77179/143‑1 (MQ=255)
                                                                                                                                  aaTGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCg                 >  2:219365/1‑143 (MQ=255)
                                                                                                                                    gggCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCgtggtggtg         <  1:92266/148‑1 (MQ=255)
                                                                                                                                            gCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGt  <  1:274044/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGT  >  CP009273/877608‑877895

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