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

AGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGT  >  CP009273/877608‑877895
                                                                                                                                             |                                                                                                                                                   
aGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGa                                                                                                                                             >  1:87266/1‑149 (MQ=255)
      ttACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGATGGCGTGAAATTca                                                                                                                                       >  1:156952/1‑149 (MQ=255)
      ttACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGAACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGGTTGAAAATGG‑CGTTGATGGCGTGAAATTca                                                                                                                                       >  1:196449/1‑149 (MQ=255)
           ggCAACTCGTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTc                                                                                                                                        <  2:217769/143‑1 (MQ=255)
                 tctTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAg                                                                                                                            <  1:366536/149‑1 (MQ=255)
                 tcctCAAGATTAAGACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAg                                                                                                                            <  1:14144/146‑1 (MQ=255)
                   ttCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGaa                                                                                                                          >  1:170879/1‑149 (MQ=255)
                        gATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGACATAATCGGGTAGAAGTTATTCAGCGTTTGCTGACAATCGGTGGGTTTGAAAATGG‑CGGTGCTGGCGAGAGAATCACCGCCGCGTAGACAcgcg                                                                                                                     >  1:124153/1‑149 (MQ=255)
                                         ccTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTtctc                                                                                                    <  1:145539/149‑1 (MQ=255)
                                            tcaccgACATGCCCGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATCATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCttt                                                                                                 <  1:147933/144‑1 (MQ=255)
                                             tACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGACAGTGGACGAT‑CAGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTc                                                                                                >  1:86027/1‑149 (MQ=255)
                                                tGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAgc                                                                                             >  1:289301/1‑149 (MQ=255)
                                                     tGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGCTTCTCTTTCAg                                                                                              >  2:245594/1‑143 (MQ=255)
                                                                cgcAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGGTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAgg                                                                                   >  2:16797/1‑143 (MQ=255)
                                                                 gcAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGt                                                                                  >  2:160076/1‑143 (MQ=255)
                                                                  cAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGc                                                                           >  1:58559/1‑149 (MQ=255)
                                                                  cAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGc                                                                           >  1:357947/1‑149 (MQ=255)
                                                                                       ggCTAAAATTATTCAGCGGTTGCTGAAAATCGATGGGTTTGAAAAGGG‑CGTTGATGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAAcggcg                                                            <  2:173295/143‑1 (MQ=255)
                                                                                       ggCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGATGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAAcggcg                                                            <  2:138532/143‑1 (MQ=255)
                                                                                        gngaaaaTTATTCATCGTTTGCCGAAACTGGATTGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTGACAGCCGAACAGCGCCGGCTCTCAATTGGGTTTCTCGGTCAGCGCCAGAAAGGTACTTATCGAATAACAAAcggcgg                                                           <  2:318602/140‑1 (MQ=255)
                                                                                         cAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAgg                                                    <  1:181855/149‑1 (MQ=255)
                                                                                                  ttCCGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCAtt                                           <  1:172023/149‑1 (MQ=255)
                                                                                                    cAGCGTTTGCGGAAAATCGATGGGTTGTAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGAttt                                               <  2:10188/143‑1 (MQ=255)
                                                                                                                      gATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCaa                             <  2:77179/143‑1 (MQ=255)
                                                                                                                                  aaTGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCg                 >  2:219365/1‑143 (MQ=255)
                                                                                                                                    ggg‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCgtggtggtg         <  1:92266/148‑1 (MQ=255)
                                                                                                                                            gCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGt  <  1:274044/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGTGTTTTACCGGCAACTCTTCAAGATTAAAACAGGTAGCGCCTTTACTGACATGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGG‑CGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGT  >  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: ‑