Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 76,319 C→T A17A (GCC→GCT)  SAUSA300_RS00345 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779376,3190CT93.1% 107.1 / 0.7 58A17A (GCC→GCT) SAUSA300_RS00345transposase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (7/47);  total (9/49)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

AAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCG  >  NC_007793/76225‑76414
                                                                                              |                                                                                               
aaGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTgg                                                                                            <  1:2967658/100‑1 (MQ=32)
                                                         gATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACg                                   <  1:208621/100‑1 (MQ=17)
                                                            ataAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgg                                <  1:1305875/100‑1 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               <  1:4360896/100‑1 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               <  1:3382833/100‑1 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               >  1:3272674/1‑100 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               <  1:4851102/100‑1 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               <  1:571103/100‑1 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggg                               <  1:5153223/100‑2 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggc                               <  1:5434473/100‑2 (MQ=17)
                                                             taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgcg                               <  1:5752216/100‑3 (MQ=17)
                                                             aaaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                               >  1:270589/2‑100 (MQ=17)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTTCTATCTAAGATATGCATTGAGTTATCGTGTTATATCTGAAATATTAAGGGAACgtggtg                              <  1:2743422/100‑1 (MQ=255)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATTGTGATATATCTGAAATATTAAGGGAACgtggtg                              >  1:5126565/1‑100 (MQ=255)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                              <  1:4089008/100‑1 (MQ=17)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                              <  1:5743813/100‑1 (MQ=17)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                              <  1:1521004/100‑1 (MQ=17)
                                                              aaaCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                              >  1:4992168/1‑100 (MQ=32)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGTTATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:3363738/100‑1 (MQ=12)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtgg                             <  1:303990/100‑2 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgggta                             <  1:644375/100‑4 (MQ=12)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:5216920/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:5668262/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:4833222/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:600857/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:4297084/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:4150240/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:1204234/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:4011093/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:1134671/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:3090941/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:3309819/100‑1 (MQ=17)
                                                               aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                             <  1:3415491/100‑1 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgtg                            <  1:922632/100‑2 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                            <  1:4738863/100‑1 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                            <  1:1421987/100‑1 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                            <  1:5692981/100‑1 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                            <  1:199597/100‑1 (MQ=17)
                                                                aCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                            <  1:4334572/100‑1 (MQ=32)
                                                                aCAATTTAACAAGAATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgtg                            <  1:1388618/100‑2 (MQ=255)
                                                                 cAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTag                           <  1:407015/100‑2 (MQ=17)
                                                                 cAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaa                           <  1:1734429/100‑1 (MQ=17)
                                                                  aaTTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaag                          <  1:3903814/100‑2 (MQ=17)
                                                                        aCAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtca                    >  1:1578454/1‑100 (MQ=18)
                                                                          aaGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtcatc                  <  1:4985052/100‑1 (MQ=18)
                                                                                gTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCaa            >  1:5228943/1‑100 (MQ=18)
                                                                                  tATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACg          >  1:4576059/1‑100 (MQ=32)
                                                                                    tctcTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGt        >  1:5005924/4‑100 (MQ=12)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:5586135/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:2432066/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:2868129/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:6161324/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:350354/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:75056/100‑1 (MQ=17)
                                                                                      aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt      <  1:4304184/100‑1 (MQ=17)
                                                                                       cTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTa     <  1:2716940/100‑1 (MQ=17)
                                                                                        tGTAGCTGTTGGTTTCTATCTAAGATATGCATTTAGTTAACGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTag    <  1:2004645/100‑2 (MQ=11)
                                                                                        tGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTag    <  1:5469457/100‑2 (MQ=17)
                                                                                          tAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCg  >  1:2879779/1‑100 (MQ=17)
                                                                                              |                                                                                               
AAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCG  >  NC_007793/76225‑76414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑