Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,108,830	A→G	H143H (CAT→CAC)	rraA ←	ribonuclease E activity regulator RraA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,108,830	0	A	G	87.5%	53.9 / 2.7	24	H143H (CAT→CAC)	rraA	ribonuclease E activity regulator RraA
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (16/5);  total (18/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACC  >  NZ_CP009273/4108743‑4108911                                                                                        |                                                                                  gAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTg                                                                                 >  1:372477/1‑90 (MQ=255)             cTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaa                                                                     >  2:65461/1‑90 (MQ=255)             cTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaa                                                                     >  2:60073/1‑90 (MQ=255)              ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                    >  2:119692/1‑90 (MQ=255)                      ccTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGac                                                            >  1:214582/1‑90 (MQ=255)                        tttCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGacac                                                          >  2:20374/1‑90 (MQ=255)                                           cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                       >  2:208814/1‑90 (MQ=255)                                            aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                      >  2:547323/1‑90 (MQ=255)                                             gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                     <  1:15254/90‑1 (MQ=255)                                             gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                     <  2:224427/90‑1 (MQ=255)                                                gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                           >  2:274212/1‑81 (MQ=255)                                                gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                           <  1:274212/81‑1 (MQ=255)                                                  tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                >  1:573984/1‑90 (MQ=255)                                                  tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                >  1:11290/1‑90 (MQ=255)                                                  tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                <  2:144663/90‑1 (MQ=255)                                                         aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                         <  2:255799/90‑1 (MQ=255)                                                         aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                         <  1:296487/90‑1 (MQ=255)                                                          aaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAt                        >  2:527282/1‑90 (MQ=255)                                                             aataatCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc                     >  1:40228/1‑90 (MQ=255)                                                                           tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc       >  2:323168/1‑90 (MQ=255)                                                                           tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc       >  1:209214/1‑90 (MQ=255)                                                                            gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc      >  2:255916/1‑90 (MQ=255)                                                                            gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc      >  2:507556/1‑90 (MQ=255)                                                                                gCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAAcc  >  2:385481/1‑90 (MQ=255)                                                                                        |                                                                                  GAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACC  >  NZ_CP009273/4108743‑4108911

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