Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 738,636 A→G E182G (GAG→GGG)  ybfD → H repeat‑associated putative transposase YbfD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913738,6360AG97.4% 116.5 / ‑4.1 38E182G (GAG→GGG) ybfDH repeat‑associated putative transposase YbfD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (20/17);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATT  >  NC_000913/738500‑738766
                                                                                                                                        |                                                                                                                                  
aGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 >  1:175904/1‑140 (MQ=38)
aGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 <  1:438626/140‑1 (MQ=255)
  tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaaga                                                                                                                               >  1:274159/1‑140 (MQ=38)
  tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 <  1:522944/138‑1 (MQ=255)
  tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 <  1:80837/138‑1 (MQ=37)
  tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 >  2:522944/1‑138 (MQ=37)
  tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGggaa                                                                                                                                 >  2:80837/1‑138 (MQ=37)
    tGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             <  1:121299/140‑1 (MQ=255)
    tGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             <  1:149544/140‑1 (MQ=255)
          tCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             <  1:268568/134‑1 (MQ=255)
          tCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             >  2:268568/1‑134 (MQ=38)
            ggACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             >  1:364970/1‑132 (MQ=37)
            ggACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATa                                                                                                                             <  2:364970/132‑1 (MQ=255)
                aGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGa                                                                                                                 >  1:528624/1‑140 (MQ=21)
                 gATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAg                                                                                                                <  2:340930/140‑1 (MQ=25)
                  aTCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAgg                                                                                                               >  2:399674/1‑140 (MQ=21)
                    cAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTg                                                                                                             <  1:225582/140‑1 (MQ=25)
                            aTAAGAAATCCAATGAGATCATAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGAttattta                                                                                                     >  1:322664/1‑140 (MQ=19)
                            aTAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGAttattta                                                                                                     >  2:42903/1‑140 (MQ=21)
                            aTAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGAttattta                                                                                                     >  2:143376/1‑140 (MQ=21)
                              aaGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATtatttatt                                                                                                   <  1:450237/140‑1 (MQ=21)
                                         tGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAg                                                                                        <  2:374878/140‑1 (MQ=21)
                                         tGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAg                                                                                        <  2:175904/140‑1 (MQ=21)
                                          gagaTCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAgg                                                                                       >  2:239930/1‑140 (MQ=25)
                                                    cTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAgggg                                                                             <  1:42903/140‑1 (MQ=21)
                                                       tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGggcgg                                                                          >  2:195016/1‑140 (MQ=25)
                                                        cccTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGggcggc                                                                         <  1:237215/140‑1 (MQ=21)
                                                            gAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCtaaa                                                                     >  1:108502/1‑140 (MQ=255)
                                                            gAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCtaaa                                                                     >  1:130292/1‑140 (MQ=255)
                                                                  cttAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGc                                                               >  1:492646/1‑140 (MQ=255)
                                                                  cttAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGc                                                               >  2:344967/1‑140 (MQ=255)
                                                                          gCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAAGGGCGGCTAAATAAAGCCTTCGAgg                                                       >  2:269463/1‑140 (MQ=255)
                                                                            tGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGaa                                                     >  2:49086/1‑140 (MQ=255)
                                                                                        gAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAATAAAGCCTTCGAGGAAAAATTTCCGCTGa                                        >  2:320654/1‑140 (MQ=255)
                                                                                                   aaaaCCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTaaa                              <  2:322664/140‑1 (MQ=255)
                                                                                                     aaCCGATGCGATGGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAaata                            >  2:113321/1‑140 (MQ=38)
                                                                                                               atgGGTTGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCa                  >  2:281156/1‑140 (MQ=255)
                                                                                                                    ttGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGa               >  1:218523/1‑138 (MQ=255)
                                                                                                                    ttGCCAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGa               <  2:218523/138‑1 (MQ=255)
                                                                                                                        cAGAAAGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTa         <  1:195016/140‑1 (MQ=255)
                                                                                                                        cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTa         <  2:31283/140‑1 (MQ=255)
                                                                                                                             aGATATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCaa    <  1:113321/140‑1 (MQ=255)
                                                                                                                               atatTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAAtt  <  2:108502/140‑1 (MQ=255)
                                                                                                                                 atTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGa               <  2:316818/125‑1 (MQ=255)
                                                                                                                                 atTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGa               >  1:316818/1‑125 (MQ=255)
                                                                                                                                        |                                                                                                                                  
AGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTCGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATT  >  NC_000913/738500‑738766

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