Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	737,899	A→G	N38N (AAT→AAC)	dtpD ←	dipeptide permease DtpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	737,899	0	A	G	93.8%	44.9 / ‑3.2	16	N38N (AAT→AAC)	dtpD	dipeptide permease DtpD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (7/8);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  NZ_CP009273/737781‑738033                                                                                                                       |                                                                                                                                       ccATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATGGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcag                                                                                                                    <  1:188933/139‑1 (MQ=255)    tGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcagcag                                                                                                                 >  2:110923/1‑139 (MQ=255)                                                     gACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                                <  2:83819/139‑1 (MQ=255)                                                      aCATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                                <  2:425883/138‑1 (MQ=255)                                                       cATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGa                                                               >  2:286831/1‑138 (MQ=255)                                                                   aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgc                                                  >  1:415526/1‑139 (MQ=255)                                                                      aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                                                                        <  2:104673/82‑1 (MQ=255)                                                                      aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                                                                        >  1:104673/1‑82 (MQ=255)                                                                      aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                            >  2:182252/1‑126 (MQ=255)                                                                      aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                            <  1:182252/126‑1 (MQ=255)                                                                      aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                               <  1:61994/139‑1 (MQ=255)                                                                          ggCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                           <  1:60977/139‑1 (MQ=255)                                                                                aaaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGAt                                     >  2:150443/1‑139 (MQ=255)                                                                                      cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt                               <  1:110923/139‑1 (MQ=255)                                                                                           aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt                          >  2:215631/1‑139 (MQ=255)                                                                                                  aGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                    >  2:224860/1‑138 (MQ=255)                                                                                                                    ttggttggTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc  >  1:162203/1‑138 (MQ=255)                                                                                                                       |                                                                                                                                       CCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  NZ_CP009273/737781‑738033

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