Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,589,954 A→C K28N (AAA→AAC)  SAUSA300_RS07845 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,589,9540AC100.0% 124.1 / NA 47K28N (AAA→AAC) SAUSA300_RS07845hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (31/16);  total (31/16)

GTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACA  >  NC_007793/1589857‑1590047
                                                                                                 |                                                                                             
gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACtt                                                                                             >  1:2350372/1‑100 (MQ=255)
gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGATACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACtt                                                                                             >  1:4786197/1‑100 (MQ=255)
  ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGAAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTaa                                                                                           >  1:5107228/1‑100 (MQ=255)
  ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAACCTCCTCGTTGTAACCCCATTAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACTTaa                                                                                           >  1:3111504/1‑100 (MQ=255)
            ttATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTAAAAAGCACttt                                                                                 >  1:3428327/1‑100 (MQ=255)
            ctATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACttt                                                                                 >  1:3021768/2‑100 (MQ=255)
             tgtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                                                                >  1:1009260/8‑100 (MQ=255)
             tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                                                                >  1:4147905/1‑100 (MQ=255)
             tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                                                                >  1:4847659/1‑100 (MQ=255)
              gtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:1486363/7‑100 (MQ=255)
              ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:4542906/7‑100 (MQ=255)
              ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:3590381/7‑100 (MQ=255)
              atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:5030204/7‑100 (MQ=255)
              atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:6167256/7‑100 (MQ=255)
              atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                                               >  1:847219/7‑100 (MQ=255)
               tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                                                              >  1:6307574/6‑100 (MQ=255)
               tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                                                              <  1:2546694/95‑1 (MQ=255)
               tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGATTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                                                              >  1:5653315/6‑100 (MQ=255)
               tGGCAGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                                                              <  1:1395034/100‑1 (MQ=255)
                ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                                                                             >  1:227284/5‑100 (MQ=255)
                ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                                                                             >  1:141162/5‑100 (MQ=255)
                          aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTa                                                                   <  1:6159713/100‑1 (MQ=255)
                          aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTCTATCAGAAAAAGCTa                                                                   <  1:4160601/100‑1 (MQ=255)
                               aaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAGAGCTAAGatc                                                              >  1:5303032/1‑98 (MQ=255)
                                    aagaagCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGa                                                         >  1:365721/1‑100 (MQ=255)
                                       tccgatctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTa                                                      <  1:1631534/94‑1 (MQ=255)
                                         gctactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat                                                    >  1:2743385/2‑98 (MQ=255)
                                         agtactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat                                                    >  1:5917729/3‑98 (MQ=255)
                                            actCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATatt                                                 >  1:2700871/1‑100 (MQ=255)
                                             ctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATattc                                                <  1:1716201/100‑1 (MQ=255)
                                                             aaaCAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATaaa                                >  1:3159360/1‑100 (MQ=255)
                                                                   aaGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCaag                          <  1:1511918/100‑2 (MQ=255)
                                                                            gtgtgcgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAATAGCTAAGCAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATttgt                 >  1:4814450/11‑100 (MQ=255)
                                                                                ccgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTgttgtt             <  1:2796794/94‑1 (MQ=255)
                                                                                ccgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTgttgtt             >  1:4454824/7‑100 (MQ=255)
                                                                                 cgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGttt            <  1:2144643/95‑1 (MQ=255)
                                                                                  gggtgtcAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTa           <  1:1604523/94‑1 (MQ=255)
                                                                                  gggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTa           <  1:5193890/96‑1 (MQ=255)
                                                                                      gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGTTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt       >  1:576476/1‑100 (MQ=255)
                                                                                      gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt       <  1:4842213/100‑1 (MQ=255)
                                                                                      gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt       >  1:4050606/1‑100 (MQ=255)
                                                                                      gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt       >  1:4010836/1‑100 (MQ=255)
                                                                                       ccAGCGAAAACTTAAAAAGCACTTTAGCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTc      >  1:6050457/1‑100 (MQ=255)
                                                                                        cAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCt     <  1:566346/100‑1 (MQ=255)
                                                                                         aGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTa    <  1:3789678/100‑1 (MQ=255)
                                                                                          gCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTAc   <  1:2217506/100‑1 (MQ=255)
                                                                                           cGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCtgaa  <  1:1576743/100‑4 (MQ=255)
                                                                                                 |                                                                                             
GTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACA  >  NC_007793/1589857‑1590047

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