Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,074,815 G→A D71N (GAT→AAT)  SAUSA300_RS05275 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,074,8150GA100.0% 104.2 / NA 37D71N (GAT→AAT) SAUSA300_RS05275hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/12);  total (25/12)

GTGGTAAAAAAGTTGATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTGATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACTA  >  NC_007793/1074720‑1074909
                                                                                               |                                                                                              
gtggtAAAAAAGTTGATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAAttt                                                                                            >  1:3570499/1‑100 (MQ=255)
           gTTGATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTAtagtta                                                                                 <  1:832843/100‑1 (MQ=255)
              gATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTg                                                                              <  1:864531/100‑1 (MQ=255)
              gATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTg                                                                              >  1:22435/1‑100 (MQ=255)
               aTACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGc                                                                             <  1:1127973/100‑1 (MQ=255)
               aTACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGc                                                                             <  1:196865/100‑1 (MQ=255)
               aTACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGc                                                                             <  1:789064/100‑1 (MQ=255)
               aTACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGc                                                                             <  1:3501985/100‑1 (MQ=255)
                tACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCa                                                                            <  1:1174165/100‑1 (MQ=255)
                  cacaTGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCagt                                                                          <  1:994633/100‑3 (MQ=255)
                  cacaTGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCaag                                                                          <  1:1049327/100‑2 (MQ=255)
                   acaTGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAAtt                                                                         >  1:2570881/1‑100 (MQ=255)
                    caTGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATtg                                                                        >  1:3684913/1‑100 (MQ=255)
                         ttgtAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATa                                                                   >  1:5590177/4‑100 (MQ=255)
                                ttgtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTatat                                                            >  1:4641726/2‑100 (MQ=255)
                                gtgtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTatat                                                            >  1:2322419/1‑100 (MQ=255)
                                gtgtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTatat                                                            >  1:514460/1‑100 (MQ=255)
                                 tgtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATGGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATAtt                                                           >  1:4223330/1‑100 (MQ=255)
                                 tgtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATAtt                                                           >  1:1020697/1‑100 (MQ=255)
                                  gtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTc                                                          >  1:3329949/1‑100 (MQ=255)
                                  gtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTc                                                          >  1:5685040/1‑100 (MQ=255)
                                  gtgGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTc                                                          >  1:4117363/1‑100 (MQ=255)
                                    gggAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCgg                                                        >  1:4728015/1‑100 (MQ=255)
                                    gggAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCgg                                                        >  1:4729205/1‑100 (MQ=255)
                                    gggAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATACTCgg                                                        <  1:1506515/100‑1 (MQ=255)
                                      gAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTa                                                      <  1:4766233/100‑1 (MQ=255)
                                              gAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGtt                                              >  1:3191117/1‑100 (MQ=255)
                                              gAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGtt                                              >  1:3764617/1‑100 (MQ=255)
                                              aaaTGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGtt                                              >  1:3428169/2‑100 (MQ=255)
                                                  gCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGt                                          <  1:5489653/100‑1 (MQ=255)
                                                        ttCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTg                                    >  1:5380217/1‑100 (MQ=255)
                                                                             ctGAGTCGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTa               >  1:3250729/2‑100 (MQ=255)
                                                                                   cGATTTATAGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAga         >  1:742578/1‑99 (MQ=255)
                                                                                         ataGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACt   >  1:3604557/1‑100 (MQ=255)
                                                                                         ataGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACt   >  1:2279425/1‑100 (MQ=255)
                                                                                         ataGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACt   >  1:1812268/1‑100 (MQ=255)
                                                                                          taGGTAATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACTa  >  1:3464418/1‑100 (MQ=255)
                                                                                               |                                                                                              
GTGGTAAAAAAGTTGATACACATGTTCATAATGTGTGGGAATTTCTGAATGCATTTTTCAGTGAAATTATTAAAGTAATGAGTCGATTTATAGGTGATTTCCCTATAGTTAGTGCAATTGTGATAATTATATTCGGTATTTTAGTTATGTTGATTGGTCATACATTACTTAGAACTATTAAGTATGACTA  >  NC_007793/1074720‑1074909

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