Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,022,294 G→A intergenic (+171/‑100) SAUSA300_RS05005 → / → SAUSA300_RS05010 CAAX amino protease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,022,2940GA100.0% 99.7 / NA 37intergenic (+171/‑100)SAUSA300_RS05005/SAUSA300_RS05010CAAX amino protease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/27);  total (10/27)

ATAGTGGTAATATTTTTTAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGGTATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTATGGCAT  >  NC_007793/1022208‑1022362
                                                                                      |                                                                    
ctaGTGGTAATATTTTTTAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTg                                                         <  1:293321/99‑1 (MQ=255)
    tGGTAATATTTTTTAACAGCAAGCGGTTATTATTGAGCTTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCa                                                     >  1:909467/1‑100 (MQ=255)
    tGGTAATATTTTTTAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCa                                                     <  1:2454060/100‑1 (MQ=255)
       tAATATTTTTTAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTGATAAAAATGATATAATATTTTTGAGCAAGt                                                  <  1:2138177/100‑1 (MQ=255)
                 tAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtg                                        <  1:6454200/100‑1 (MQ=255)
                  aaCAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtgg                                       <  1:2053838/100‑1 (MQ=255)
                  aaCAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTGTGAGCAAGTTGAATTGAtgg                                       <  1:2244469/100‑1 (MQ=255)
                   ccAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggt                                      >  1:6434676/2‑100 (MQ=255)
                   ccAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggt                                      >  1:3901193/2‑100 (MQ=255)
                   aCAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggg                                      <  1:4198359/100‑2 (MQ=255)
                    cAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTGCTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggag                                     <  1:1834102/100‑3 (MQ=255)
                    cAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtg                                     >  1:3545399/1‑100 (MQ=255)
                    cAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtg                                     <  1:407071/100‑1 (MQ=255)
                    cAGCAAGAGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtg                                     >  1:1640648/1‑100 (MQ=255)
                     tgCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtgg                                    >  1:5823188/2‑100 (MQ=255)
                      gCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtgga                                   <  1:3213248/100‑2 (MQ=255)
                      gCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGAtggtgag                                   >  1:300506/1‑98 (MQ=255)
                      gCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGc                                   <  1:1593476/100‑1 (MQ=255)
                      gCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATAGTGGc                                   >  1:3122214/1‑100 (MQ=255)
                       cAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCt                                  <  1:5056125/100‑1 (MQ=255)
                       cAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCt                                  >  1:4460622/1‑100 (MQ=255)
                       cAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATA‑TTTTGAGCAAGTTGAATTGATGGTGGCTa                                 <  1:2224281/100‑1 (MQ=255)
                        aaGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTa                                 >  1:3856777/1‑100 (MQ=255)
                            ggTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTg                             <  1:3054535/100‑1 (MQ=255)
                              ttattaTTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTggg                           <  1:2522778/100‑3 (MQ=255)
                                     tGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATA‑TTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAgggg                   <  1:3228180/100‑1 (MQ=255)
                                        gaTTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGGGGCTATCTGAGTAAAgggggg                 <  1:3863443/100‑1 (MQ=255)
                                         aTTAGGTTTTAGTCAG‑TTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGggtg               <  1:4930955/100‑1 (MQ=255)
                                          ttAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGggtg               <  1:5875491/100‑1 (MQ=255)
                                          ttAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGggtg               <  1:2102054/100‑1 (MQ=255)
                                              gTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGc           <  1:4356584/100‑1 (MQ=255)
                                                ctctgTCAGTTTTTATTACTAAAGATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGT‑AAGGGGGGTGGTGCCTa        <  1:6370726/96‑1 (MQ=38)
                                                atctgTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCt         <  1:4186800/96‑1 (MQ=255)
                                                 ttAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTa        <  1:5985623/100‑1 (MQ=255)
                                                  tcgTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTAt       <  1:3733960/98‑1 (MQ=255)
                                                  tAGTCAGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTag       <  1:3471887/100‑2 (MQ=255)
                                                       aGTTTTTATTACTAAATATTTTATAAAAATGATATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTATGGCag  <  1:2483055/100‑2 (MQ=255)
                                                                                      |                                                                    
ATAGTGGTAATATTTTTTAACAGCAAGCGGTTATTATTGAGATTAGGTTTTAGTCAGTTTTTATTACTAAATATTTTATAAAAATGGTATAATATTTTTGAGCAAGTTGAATTGATGGTGGCTATCTGAGTAAAGGGGGGTGGTGCCTATGGCAT  >  NC_007793/1022208‑1022362

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