Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,532,756 A→G D314G (GAC→GGC)  ydcC → H repeat‑associated putative transposase YdcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,532,7560AG97.8% 142.7 / ‑4.2 45D314G (GAC→GGC) ydcCH repeat‑associated putative transposase YdcC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/22);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCA  >  NC_000913/1532620‑1532866
                                                                                                                                        |                                                                                                              
ataataGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGt                                                                                                             >  2:71644/1‑140 (MQ=255)
      gCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                       <  1:71644/140‑1 (MQ=255)
       cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                       <  2:134351/139‑1 (MQ=255)
       cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                       >  1:134351/1‑139 (MQ=255)
         gAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaatgaa                                                                                                    >  1:374066/1‑140 (MQ=255)
                    agagCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACt                                                                                         <  2:418966/140‑1 (MQ=255)
                    agagCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACt                                                                                         >  2:659144/1‑140 (MQ=25)
                      agCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaatga                                                                                                     >  2:465349/1‑126 (MQ=255)
                      agCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaatga                                                                                                     <  1:465349/126‑1 (MQ=255)
                                    gTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                                                                         <  1:746689/140‑1 (MQ=21)
                                       aGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                                                                         >  1:522369/1‑137 (MQ=25)
                                       aGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                                                                         <  2:522369/137‑1 (MQ=21)
                                               tatCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGcagcag                                                              >  2:634979/1‑140 (MQ=255)
                                               tatCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGcagcag                                                              >  1:34264/1‑140 (MQ=255)
                                               tatCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGcagcag                                                              >  2:742744/1‑140 (MQ=255)
                                                   aGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatt                                                          >  2:245256/1‑140 (MQ=255)
                                                    gTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                                         >  1:688478/1‑140 (MQ=255)
                                                     ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcag                                                                 <  2:247867/131‑1 (MQ=255)
                                                     ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcag                                                                 >  1:247867/1‑131 (MQ=255)
                                                     ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAattat                                                        >  1:402970/1‑140 (MQ=255)
                                                       ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAttt                                                      >  1:831953/1‑140 (MQ=255)
                                                       ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAttt                                                      >  1:670076/1‑140 (MQ=255)
                                                        tgctgcTTTCACCGCTGAGAGGTTCTGCACCGCCCTCCCAATCCACTGGCGTGGGGGAAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAATACGCCTGCAAAATGAGAAGAGGGAATGCAGCAGAAATAtttt                                                     >  2:422500/1‑140 (MQ=2)
                                                         gctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                                    <  2:670076/140‑1 (MQ=255)
                                                               ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                  >  2:306703/1‑136 (MQ=255)
                                                               ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                  <  1:306703/136‑1 (MQ=255)
                                                                tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGAtg                                             <  2:322996/140‑2 (MQ=255)
                                                                  aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACg                                           >  1:558902/1‑140 (MQ=255)
                                                                     gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                        <  2:831953/140‑1 (MQ=255)
                                                                      cAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcac                                       <  1:235801/140‑1 (MQ=255)
                                                                        gagaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                     <  2:143430/140‑1 (MQ=255)
                                                                         agaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                    <  2:402970/140‑1 (MQ=255)
                                                                                cGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa                             >  1:630428/1‑140 (MQ=255)
                                                                                   cacaGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                     <  2:87992/129‑1 (MQ=38)
                                                                                   cacaGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                     >  1:87992/1‑129 (MQ=255)
                                                                                   cacaGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAata                          >  2:216692/1‑140 (MQ=255)
                                                                                     caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtt                        >  2:293907/1‑140 (MQ=255)
                                                                                     caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtt                        >  2:225218/1‑140 (MQ=255)
                                                                                        cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTg                     >  2:864563/1‑140 (MQ=38)
                                                                                          aTCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                   <  2:724975/140‑1 (MQ=21)
                                                                                           tCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                  >  2:451002/1‑140 (MQ=37)
                                                                                                  ccACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATa           <  1:864563/140‑1 (MQ=18)
                                                                                                    aCTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAg         <  2:558902/140‑1 (MQ=17)
                                                                                                    aCTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAg         <  1:742744/140‑1 (MQ=17)
                                                                                                         cATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAtt    <  2:34264/140‑1 (MQ=17)
                                                                                                         cATGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAtt    <  2:630428/140‑1 (MQ=17)
                                                                                                           tgtgGAGAATAAGCTGCACTGGCGTCTGGGCGTGTTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCa  <  1:366066/140‑1 (MQ=11)
                                                                                                                                        |                                                                                                              
ATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCA  >  NC_000913/1532620‑1532866

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