Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 796,114 C→A 22.6% S345Y (TCT→TAT)  recQ1 → ATP‑dependent helicase RecQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730796,1140CA22.6% 36.1 / 4.3 17S345Y (TCT→TAT) recQ1ATP‑dependent helicase RecQ
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base A (2/2);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTTGTTATTCATTATAATATGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTAGAAGAATGTGAACAATGTAGT  >  CP000730/795966‑796259
                                                                                                                                                    |                                                                                                                                                 
cGCTTTGTTATTCATTATAATATGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGtctc                                                                                                                                                   >  1:100818/1‑149 (MQ=255)
 tcTTTGTGATTCAGTATAATATGCCTGGAGATTTAGAATCTTATTATCAAGAATCGGGTCGTGCATGTCGTGACGGGTTGAAAAGTGAATGTATTTTATTATTTAAAGGACGTGATATCAAGTTACACGAGTATTTTATAACAGtctct                                                                                                                                                  <  2:241073/148‑1 (MQ=255)
          ttCATTANAATATGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGTAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCtgat                                                                                                                                         <  2:226296/149‑1 (MQ=255)
                 taataTGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACt                                                                                                                                  <  2:292448/149‑1 (MQ=255)
                  aataTGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTATCAAGCTGATGATGACta                                                                                                                                 >  2:17604/1‑149 (MQ=255)
                  aataTGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTATCAAGCTGATGATGACta                                                                                                                                 >  2:17690/1‑149 (MQ=255)
                    tatGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACtata                                                                                                                               >  2:237885/1‑149 (MQ=255)
                     atGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATaa                                                                                                                              <  1:258525/149‑1 (MQ=255)
                                     aaTCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGaa                                                                                                              >  1:278084/1‑149 (MQ=255)
                                        cTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTGACACTAGTATTTGATAACAGTCTCTCAAGCTTATGATGACTATAAAGATAAAATGGGCGAAAAg                                                                                                           >  1:301902/1‑149 (MQ=255)
                                              atCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACt                                                                                                     >  2:114447/1‑149 (MQ=255)
                                                             gTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATACAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAatat                                                                                      <  2:107354/149‑1 (MQ=255)
                                                                cAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATaca                                                                                   >  1:14535/1‑149 (MQ=255)
                                                                           gggTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACaaaaaa                                                                        >  2:169124/1‑149 (MQ=255)
                                                                                              ttttGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAAtt                                                     >  2:181593/1‑149 (MQ=255)
                                                                                                           cGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTg                                        >  1:262483/1‑149 (MQ=255)
                                                                                                                                cGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTagaag                   >  1:215651/1‑149 (MQ=255)
                                                                                                                                   gTATTTTATAACAGTCTATCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTAGAAGAAt                <  1:17690/149‑1 (MQ=255)
                                                                                                                                   gTATTTTATAACAGTCTATCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTAGAAGAAt                <  1:17604/149‑1 (MQ=255)
                                                                                                                                                 tctctcAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTAGAAGAATGTGAACAATgtagt  >  1:25541/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCTTTGTTATTCATTATAATATGCCTGGAGATTTAGAATCTTATTATCAAGAAGCGGGTCGTGCAGGTCGTGACGGGTTGAAAAGTGAATGTATTTTGTTATTTAGCGAACGCGATATCAATTTACACGAGTATTTTATAACAGTCTCTCAAGCTGATGATGACTATAAAGATAAAATGGGCGAAAAGTTAACTAAAATGATTCAATATACAAAAACAAAAAAATGTCTAGAAGCAACAATTGTCCATTATTTTGAACCGAATGAAAAATTAGAAGAATGTGAACAATGTAGT  >  CP000730/795966‑796259

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