Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,000,655 A→G E340G (GAA→GGA)  amyA → alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,000,6550AG100.0% 72.8 / NA 21E340G (GAA→GGA) amyAalpha‑amylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/12);  total (9/12)

CATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGG  >  NZ_CP009273/2000532‑2000782
                                                                                                                           |                                                                                                                               
cATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAAcc                                                                                                                  <  2:378439/139‑1 (MQ=255)
   cgcgCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTg                                                                                                               >  1:297312/1‑139 (MQ=255)
                      cGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTa                                                                                             <  2:235005/138‑1 (MQ=255)
                                             gTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGaaaa                                                                     >  2:364149/1‑139 (MQ=255)
                                                        gAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCtt                                                          >  1:141763/1‑139 (MQ=255)
                                                        gAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCtt                                                          <  1:288729/139‑1 (MQ=255)
                                                          aGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCg                                                        <  1:322333/139‑1 (MQ=255)
                                                          aGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCg                                                        <  2:297312/139‑1 (MQ=255)
                                                           gCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCgg                                                       >  1:74187/1‑139 (MQ=255)
                                                           gCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCAGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCgg                                                       >  2:61454/1‑139 (MQ=255)
                                                             cGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTa                                                     <  1:364149/139‑1 (MQ=255)
                                                               aTCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTAtt                                                   >  1:501632/1‑139 (MQ=255)
                                                                     tCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCg                                                                                                     <  2:476988/83‑1 (MQ=255)
                                                                     tCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCg                                                                                                     >  1:476988/1‑83 (MQ=255)
                                                                        aCGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGa                                          <  2:141763/139‑1 (MQ=255)
                                                                              tGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTa                                    <  1:26228/139‑1 (MQ=255)
                                                                                   cTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTg                               <  2:74187/139‑1 (MQ=255)
                                                                                   cTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTg                               >  1:47206/1‑139 (MQ=255)
                                                                                            aaTCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACg                      <  1:153209/139‑1 (MQ=255)
                                                                                                           ccGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGgtggtg       <  2:425249/139‑1 (MQ=255)
                                                                                                                 cAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACggg  >  2:549498/1‑138 (MQ=255)
                                                                                                                           |                                                                                                                               
CATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGG  >  NZ_CP009273/2000532‑2000782

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