Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,401,096	T→C	G19G (GGT→GGC)	alaA →	alanine transaminase AlaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,401,096	0	T	C	89.4%	122.7 / 8.8	47	G19G (GGT→GGC)	alaA	alanine transaminase AlaA
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (24/18);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

ACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGTCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAACCTGCCTAC  >  NZ_CP009273/2400959‑2401230                                                                                                                                          |                                                                                                                                       ccacCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGTc                                                                                                                                       <  1:763941/138‑1 (MQ=255)                 gTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAaagaag                                                                                                                        >  2:965426/1‑138 (MQ=255)                 gTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAaagaag                                                                                                                        >  1:1264590/1‑138 (MQ=255)                  tACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGTCCGGTGCTGAAAGAAGc                                                                                                                       >  2:754134/1‑138 (MQ=255)                  tACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGc                                                                                                                       >  1:8250/1‑138 (MQ=255)                       tGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCaaaa                                                                                                                   >  2:631341/1‑137 (MQ=255)                       tGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCaaaa                                                                                                                   <  1:631341/137‑1 (MQ=255)                         ttCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGcc                                                                                                               >  1:74171/1‑139 (MQ=255)                          tCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAaaga                                                                                                                          <  1:109503/127‑1 (MQ=255)                          tCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAaaga                                                                                                                          >  2:109503/1‑127 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  2:398541/1‑138 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  1:629411/1‑138 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  2:406804/1‑138 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  1:942243/1‑138 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  2:305601/1‑138 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      <  2:1264590/138‑1 (MQ=255)                                   aCTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGgaagaag                                                                                                      >  2:154524/1‑138 (MQ=255)                                       ccACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGTCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAgg                                                                                                  <  2:176665/138‑1 (MQ=255)                                        cACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGt                                                                                                 >  2:829231/1‑138 (MQ=255)                                        cACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGt                                                                                                 >  2:751840/1‑138 (MQ=255)                                         aCGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTaa                                                                                               <  1:1213634/139‑1 (MQ=255)                                          cGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAAc                                                                                              <  1:1166237/139‑1 (MQ=255)                                               gggAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAgg                                                                                          >  2:348684/1‑138 (MQ=255)                                                 gAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTac                                                                                       <  2:890244/139‑1 (MQ=255)                                                     ggCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTactga                                                                                    >  2:94850/1‑138 (MQ=255)                                                      gCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTactgaaa                                                                                  <  2:1172601/139‑1 (MQ=255)                                                      gCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTactgaaa                                                                                  <  2:230318/139‑1 (MQ=255)                                                      gCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTactgaaa                                                                                  >  1:259495/1‑139 (MQ=255)                                                            aCAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAAc                                                                            <  1:154524/139‑1 (MQ=255)                                                                cAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCg                                                                        >  2:547548/1‑139 (MQ=255)                                                                 aaaaaaGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCgg                                                                       <  1:1060621/139‑1 (MQ=255)                                                                          cAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGc                                                              >  1:147178/1‑139 (MQ=255)                                                                           aaGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGcc                                                             <  1:598872/139‑1 (MQ=255)                                                                                 cATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGtt                                                        <  2:336776/138‑1 (MQ=255)                                                                                  aTGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCg                                                      >  2:542885/1‑139 (MQ=255)                                                                                  aTGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCg                                                      >  2:1209727/1‑139 (MQ=255)                                                                                       ccccATTGAAAAATCCAGCAAAATAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGtttt                                                 >  2:1137937/1‑139 (MQ=255)                                                                                           aTTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGAcg                                             <  2:259495/139‑1 (MQ=255)                                                                                                 aaaTCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAg                                       <  2:942243/139‑1 (MQ=255)                                                                                                    tCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGAt                                     >  2:450510/1‑138 (MQ=255)                                                                                                                           ttATGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGc             <  2:297458/139‑1 (MQ=255)                                                                                                                              tGACATCCGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAAc          <  2:745620/139‑1 (MQ=255)                                                                                                                                   tCCGTGGTCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAAt                                <  2:422594/112‑1 (MQ=255)                                                                                                                                   tCCGTGGTCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAAt                                >  1:422594/1‑112 (MQ=255)                                                                                                                                     cGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAACCTGCCTa   <  1:829231/139‑1 (MQ=255)                                                                                                                                     cGTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAACCTGCCTa   >  1:216409/1‑139 (MQ=255)                                                                                                                                      gTGGCCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAACCTGCCTAc  <  1:348684/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       ACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAAAAAGGTCAAGATTCATGTCCCCCATTGAAAAATCCAGCAAATTAGAGAATGTCTGTTATGACATCCGTGGTCCGGTGCTGAAAGAAGCAAAACGCCTGGAAGAAGAAGGTAACAAGGTACTGAAACTGAACATCGGCAACCCAGCCCCGTTCGGTTTTGACGCGCCAGATGAAATCCTCGTTGACGTGATACGCAACCTGCCTAC  >  NZ_CP009273/2400959‑2401230

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