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_CP009273	2,000,655	0	A	G	97.1%	115.2 / ‑4.5	35	E340G (GAA→GGA)	amyA	alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/17);  total (17/18)
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.59e-01

CATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAC  >  NZ_CP009273/2000525‑2000787                                                                                                                                   |                                                                                                                                     cATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCgg                                                                                                                              <  2:150875/139‑1 (MQ=255)        cATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAAcc                                                                                                                       <  2:707093/139‑1 (MQ=255)         aTCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCg                                                                                                                      >  2:373340/1‑139 (MQ=255)           cgcgCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTg                                                                                                                    >  2:353811/1‑139 (MQ=255)                aTGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTa                                                                                                               <  1:468109/139‑1 (MQ=255)                        cgACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTg                                                                                                       >  1:678350/1‑139 (MQ=255)                           cTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCg                                                                                                    <  1:82413/139‑1 (MQ=255)                              cGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTa                                                                                                  <  1:353811/138‑1 (MQ=255)                               gACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATg                                                                                                >  2:560464/1‑139 (MQ=255)                                  aTGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCt                                                                                             <  2:468110/139‑1 (MQ=255)                                    gACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCtt                                                                                            >  2:621103/1‑138 (MQ=255)                                           aTTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTttgt                                                                                    >  1:570070/1‑139 (MQ=255)                                                                gAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCtt                                                               >  1:246625/1‑139 (MQ=255)                                                                  aGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCg                                                             <  1:238520/139‑1 (MQ=255)                                                                   gCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGc                                                                      <  1:273919/129‑1 (MQ=255)                                                                   gCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGc                                                                      >  2:273919/1‑129 (MQ=255)                                                                   gCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCgg                                                            >  2:147357/1‑139 (MQ=255)                                                                              ccACGCCGTGACGCTCGTTGCCAAACACGGCACCGAAGCGTAGCAAGCCCACGAAGCGCCGGGCGAACCGAGGTGTAAACCGCTGGCGTATGCCTTAAGGTTGGGGCGGGAAAATGGCGTACCTACGGGAATCTATCCg                                                 <  2:157500/139‑1 (MQ=255)                                                                               cACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                 <  1:373340/138‑1 (MQ=255)                                                                                      tGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                 >  1:519453/1‑131 (MQ=255)                                                                                      tGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                 <  2:519453/131‑1 (MQ=255)                                                                                      tGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTa                                         >  1:593933/1‑139 (MQ=255)                                                                                       gACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTAc                                        >  1:744442/1‑139 (MQ=255)                                                                                              gTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCgg                                                            >  1:605115/1‑112 (MQ=255)                                                                                              gTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCgg                                                            <  2:605115/112‑1 (MQ=255)                                                                                              gTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTgcg                                  >  2:408053/1‑138 (MQ=255)                                                                                              gTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTgcg                                  <  2:570070/138‑1 (MQ=255)                                                                                                                 aaCCGTTGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCggtgg              >  2:456586/1‑139 (MQ=255)                                                                                                                      ttGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACct                                            >  1:680221/1‑104 (MQ=255)                                                                                                                      ttGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACct                                            <  2:680221/104‑1 (MQ=255)                                                                                                                      ttGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGaa                         >  1:386325/1‑123 (MQ=255)                                                                                                                      ttGCAAGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGaa                         <  2:386325/123‑1 (MQ=255)                                                                                                                           aGCCCTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCaa    <  1:147357/139‑1 (MQ=255)                                                                                                                             cccTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTTGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAc  <  2:274790/139‑1 (MQ=255)                                                                                                                             cccTCGGAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAc  <  2:678350/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     CATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAC  >  NZ_CP009273/2000525‑2000787

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