Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	455,525	C→T	H394H (CAC→CAT)	lon →	endopeptidase La

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	455,525	0	C	T	95.0%	62.6 / ‑3.5	20	H394H (CAC→CAT)	lon	endopeptidase La
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (12/7);  total (13/7)
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.19e-01

GTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGC  >  NZ_CP009273/455440‑455598                                                                                      |                                                                          gTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCAtcgccg                                                                      <  1:401405/91‑1 (MQ=255)      tCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACtt                                                                 <  1:309253/91‑1 (MQ=255)      tCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACtt                                                                 <  2:266260/91‑1 (MQ=255)           tgacAAAGCCGCAGGGCATAAAGATGTCCGTATGGCGCTGGGCGGCCTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATc                                                            <  2:312739/88‑1 (MQ=255)           tGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATc                                                            <  1:670495/91‑1 (MQ=255)           tGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATc                                                            <  2:360426/91‑1 (MQ=255)               aaaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGtt                                                        >  2:694043/1‑91 (MQ=255)               aaaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGtt                                                        >  2:420112/1‑91 (MQ=255)                     aCCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGc                                                  >  2:531632/1‑91 (MQ=255)                      ccGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGcc                                                 >  2:68998/1‑91 (MQ=255)                       cGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCg                                                <  1:653515/91‑1 (MQ=255)                            gTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTaa                                           >  2:14525/1‑91 (MQ=255)                                   tGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATc                                    >  1:448420/1‑91 (MQ=255)                                   tGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATc                                    >  2:576339/1‑91 (MQ=255)                                   tGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATc                                    >  2:335808/1‑91 (MQ=255)                                    gTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                   >  2:367076/1‑91 (MQ=255)                                        gTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGaa                               >  2:255334/1‑91 (MQ=255)                                            ggCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                           >  1:174982/1‑91 (MQ=255)                                             gcgcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATgg                          >  1:480843/1‑91 (MQ=255)                                                 tGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGaa                      >  2:769330/1‑91 (MQ=255)                                                                     gCGGAAATCCGTGGTCATCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGc  <  1:423319/91‑1 (MQ=255)                                                                                      |                                                                          GTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGC  >  NZ_CP009273/455440‑455598

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