Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 46,651 G→A G290G (GGC→GGT)  SAUSA300_RS00200 ← recombinase RecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779346,6510GA98.3% 162.7 / ‑6.5 59G290G (GGC→GGT) SAUSA300_RS00200recombinase RecB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/48);  minor base T (0/1);  total (10/49)
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.95e-01

TAGTCCCGTCTTTAAGAGTATTCGTGGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAGTGAATAATGCCTGTAAGCAGATTCGTTCCTTTGCCATGGACTTGGGGTTTTTGACTGACTTGTTTCTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTC  >  NC_007793/46554‑46743
                                                                                                 |                                                                                            
tggTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTTGTTCCTTTAcc                                                                                            <  1:5381287/98‑1 (MQ=255)
tggTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTAcc                                                                                            <  1:3208722/98‑1 (MQ=255)
 ggTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCa                                                                                           <  1:3360678/99‑1 (MQ=255)
 ggTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCa                                                                                           <  1:5092173/99‑1 (MQ=255)
 ggTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCa                                                                                           >  1:579453/2‑100 (MQ=255)
  gTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCAt                                                                                          <  1:1280475/100‑1 (MQ=255)
  gTCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCCGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCAt                                                                                          <  1:3954231/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCAta                                                                                         <  1:2789511/100‑2 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:3896637/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:4018326/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:2804207/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:601284/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:786003/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATg                                                                                         <  1:751703/100‑1 (MQ=255)
   tCCCATCTTTAAGTGTATTCATTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCAtt                                                                                         <  1:255606/100‑2 (MQ=255)
    cccATCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATga                                                                                        <  1:136429/100‑2 (MQ=255)
     ccatCTTTAAGTGTATTCGTTGTGTTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCaagag                                                                                       <  1:2189396/97‑5 (MQ=255)
     ccatCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTAcccacac                                                                                       <  1:2937690/97‑6 (MQ=255)
     ccatCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATgag                                                                                       <  1:3953920/97‑3 (MQ=255)
      catCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTAccccttgg                                                                                      <  1:337996/98‑7 (MQ=255)
       atCTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCTACATTGGGGACAATGGATAATGACTGTAAGCAGATTCTTTCCTTTACCATGAACt                                                                                     <  1:222153/99‑1 (MQ=255)
         cTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTg                                                                                   <  1:539376/100‑1 (MQ=255)
         cTTTAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTg                                                                                   <  1:5708575/100‑1 (MQ=255)
          tttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTgg                                                                                  <  1:1992922/100‑1 (MQ=255)
           ttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTggg                                                                                 >  1:2815035/1‑100 (MQ=255)
           ttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTggg                                                                                 <  1:1353182/100‑1 (MQ=255)
           ttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTggg                                                                                 <  1:3387143/100‑1 (MQ=255)
           ttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTggc                                                                                 <  1:3389355/100‑2 (MQ=255)
           ttAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTgat                                                                                 <  1:777491/100‑3 (MQ=255)
            tAAGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTgggg                                                                                <  1:1500225/100‑1 (MQ=255)
             aaGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTggggg                                                                               <  1:2039236/100‑2 (MQ=255)
             aaGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGt                                                                               >  1:4959878/1‑100 (MQ=255)
             aaGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGt                                                                               <  1:885988/100‑1 (MQ=255)
             aaGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGt                                                                               <  1:971613/100‑1 (MQ=255)
             aaGTGTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGACTGTAAGCAGATTCGTTCCTTTACCATGAACTTggggg                                                                               <  1:5685764/100‑2 (MQ=255)
              agtgTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGtt                                                                              <  1:900392/97‑1 (MQ=255)
              agtgTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGtt                                                                              <  1:2708269/97‑1 (MQ=255)
               gtgTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGttt                                                                             <  1:3717183/98‑1 (MQ=255)
               gtgTATTCGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCATTTACCATGAACTTGGGGtaa                                                                             <  1:746529/98‑3 (MQ=255)
                      cGTTGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTtgactg                                                                      <  1:2384093/100‑1 (MQ=255)
                            aTTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTgg                                                                <  1:2002620/100‑2 (MQ=255)
                                    ctgcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACg                                                        >  1:3294330/1‑100 (MQ=255)
                                     tgcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGc                                                       <  1:156571/100‑1 (MQ=255)
                                      gcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCa                                                      >  1:446024/1‑100 (MQ=255)
                                      gcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCa                                                      >  1:5728413/1‑100 (MQ=255)
                                       ccgtAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCAt                                                     >  1:2791684/4‑100 (MQ=255)
                                               gcgcCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTTCCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACta                                             <  1:404656/100‑2 (MQ=255)
                                               gcgcCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACtt                                             >  1:2553641/1‑100 (MQ=255)
                                                    acaTTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATcc                                        <  1:3217511/100‑1 (MQ=255)
                                                       ttGGGGACAATTGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGAGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCAt                                     <  1:3718240/100‑1 (MQ=255)
                                                       ttGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCAt                                     <  1:3788460/100‑1 (MQ=255)
                                                              caatggataatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCtt                              <  1:3758372/94‑1 (MQ=255)
                                                               aatggataatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTg                             <  1:1912494/95‑1 (MQ=255)
                                                                     caatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATtaat                       >  1:4542211/2‑100 (MQ=255)
                                                                     agatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATtaat                       >  1:829127/3‑100 (MQ=255)
                                                                          ccTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATgg                  <  1:3902717/100‑1 (MQ=255)
                                                                                      ttCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTAc      <  1:1781260/100‑1 (MQ=255)
                                                                                         gTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCtt   <  1:3947526/100‑1 (MQ=255)
                                                                                          ttCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTAccctg  <  1:1625192/100‑4 (MQ=255)
                                                                                                 |                                                                                            
TAGTCCCGTCTTTAAGAGTATTCGTGGTATTGCTTGCTGCCATAGGTGCGCCACATTGGGGACAGTGAATAATGCCTGTAAGCAGATTCGTTCCTTTGCCATGGACTTGGGGTTTTTGACTGACTTGTTTCTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTC  >  NC_007793/46554‑46743

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