Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	990,702	T→A	18.0%	intergenic (‑54/+175)	CJ019_00935 ← / ← bfrD	hypothetical protein/putative TonB‑dependent receptor BfrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	990,702	0	T	A	18.0%	73.4 / 3.7	43	intergenic (‑54/+175)	CJ019_00935/bfrD	hypothetical protein/putative TonB‑dependent receptor BfrD
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base A (2/6);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAATGCCCATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAG  >  KT2440/990568‑990830                                                                                                                                        |                                                                                                                                  ggTAATGCCCATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATa                                                                                                                                  >  4:345914/1‑136 (MQ=255) ggTAATGCCCATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCAGTATGCCTTGAt                                                                                                                                   >  8:361600/1‑135 (MQ=255)   tAATGCCCATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATac                                                                                                                                 >  5:172357/1‑135 (MQ=255)        cccATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGCCCCCCGCGGCGT‑GCATTTTTCCTTTAAACACaa                                                                                                                             >  2:342808/1‑134 (MQ=255)           aTCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGcc                                                                                                                        >  3:43373/1‑136 (MQ=255)             cAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCTGCGCCCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGccc                                                                                                                       >  8:16820/1‑135 (MQ=255)               aGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCCCGCGGCGTGGCATT‑TGCCGTGATAGACCAAACccccc                                                                                                                     >  2:137031/1‑134 (MQ=255)                ggCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGCCCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGc                                                                                                                    >  8:38004/1‑135 (MQ=255)                ggCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGc                                                                                                                    <  1:368445/135‑1 (MQ=255)                ggCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGc                                                                                                                    <  5:285480/135‑1 (MQ=255)                ggCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCGTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGGCATT‑TGGCTTTATGCACAAAAGCCCCGc                                                                                                                    >  1:250113/1‑135 (MQ=255)                     gCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACg                                                                                                               <  2:105767/135‑1 (MQ=255)                      cACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATAAACAAAAGCCCCGCTCACg                                                                                                               >  6:299305/1‑134 (MQ=255)                         cAGGCCAGCGTTGATGCCAAAGAACCAGTGCAATTGTAACAGCGTCTTCTTCACCACGTTTCTCACATTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GAATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAAt                                                                                                          <  4:184796/136‑1 (MQ=255)                           ggCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAACGCCCCGCTCACGGGAATGa                                                                                                        <  3:345914/136‑1 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGCCCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGa                                                                                                        >  3:298637/1‑133 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGa                                                                                                        <  4:298637/133‑1 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGCAc                                                                                                      >  6:344249/1‑135 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAAc                                                                                                      >  8:229786/1‑135 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAAc                                                                                                      >  2:288608/1‑135 (MQ=255)                              cAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACg                                                                                                     >  3:128291/1‑136 (MQ=255)                                         ccAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGc                                                                                                                    >  6:69521/1‑110 (MQ=255)                                         ccAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGc                                                                                                                    <  5:69521/110‑1 (MQ=255)                                                 ccAGTGCAATTGGAACAGCGTCTTCTTCACCAGGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GAATTATGCCTTGATACACAAACGCCCCGCTCACGGGAATGACCGGGGCCGGGAGGGCTGCtt                                                                                   <  7:257243/135‑1 (MQ=255)                                                  cAGTGCAATTGCAACAGCGTCTTCTCCCCCACGTTCGACACCTTGCGTGGCTGCCGCGCCCCGCGCGGCGT‑GCATTATGCCTTGATCCACAAACCCCCCGCTCCCGGGAATCACCGGGGCCGGGAGGGCTGCtttt                                                                                 <  4:210433/136‑1 (MQ=255)                                                        aaTTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTAAGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTAGtgtt                                                                           >  3:93577/1‑136 (MQ=255)                                                        aaTTGCAACAGCGGCTTCTTCACCACGTTGCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GAATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCgtgt                                                                            <  2:274078/135‑1 (MQ=255)                                                             gAACAGCGGCTTCTTCACCACGTTTCTCCCCTTGAGTGGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCt                                                                       <  4:153817/135‑1 (MQ=255)                                                                  gCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATg                                                                                                         >  2:59861/1‑96 (MQ=255)                                                                  gCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATg                                                                                                         <  1:59861/96‑1 (MQ=255)                                                                  gCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCGCACGGGAATTAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCgg                                                                 >  3:316701/1‑136 (MQ=255)                                                                     tcttctTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGc                                                               >  7:200719/1‑135 (MQ=255)                                                                               cacGTTACTCACCTTTAGTTGCTGCCCCGCCCCGCGCGGCGT‑GCTTTTTGCCTTGAACCACAAAAGCCCCGCTCACGGCAATGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGc                                                   <  4:100180/136‑1 (MQ=255)                                                                                acctttATCACCGTGAGTTGCTGCCGCGACCCGCGCGGCCT‑GCAATATACCATGAATAACCAAACCCCCCCACACCGGAAAGAACCGGGCCCCGCGGGCAGGCTTTAGTTTTGTAATTTCGGGCCTCATTCGCgg                                                    >  8:322438/4‑135 (MQ=255)                                                                                  gttgcccACCTTTAGTGGCGGCGGCGACCCGGGCGGCGT‑GCATTATGCCTGTATACACAAAAGCCCCGCTCCCGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTaa                                                <  3:37987/130‑1 (MQ=255)                                                                                   tttCTCCCCTTGAGGTGCGGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTCAAACACAAAAGCCCCGCTCACGGGAATGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTaa                                                <  5:265816/135‑1 (MQ=255)                                                                                               tgTTGCGGCCGCGCCCCGCGCGGCGT‑GCATTATGCCTTGATACACAAACGCCCCGCGCACGGGAATGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGc                                    <  5:282853/134‑1 (MQ=255)                                                                                                          cgACCCGCGCGGGGT‑GCATTATGCCTTGAAACACAAAAGCCCCGCTCACGGGAATGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGGCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                        <  4:316701/136‑1 (MQ=255)                                                                                                          cgACCCGCGCGGCGT‑GCATTATGCCTGGAAACACAAACGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                        <  4:182763/136‑1 (MQ=255)                                                                                                           gACCCGCGCGGCGT‑GCATTATGCCTTGAAACACAAAACCCCCGCCCACGGGAATGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                        <  6:259257/135‑1 (MQ=255)                                                                                                           gACCCGCGCGGCGT‑GCATTAGGCCTTGAACCACAACCCCCCCGCCCCCGGAAAGGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                        <  1:42290/135‑1 (MQ=255)                                                                                                                     gcgT‑GAATTATCCCTTGATACACAAATGCCCCCCTCACGGCATTGACCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTc             <  3:171343/136‑1 (MQ=255)                                                                                                                                  cctgcaTACACAAACGCCCCGCTCACGGGAAGTACCGGGGCCGGGAGGGCTGCTTTTCGGGTGGTCGGTGCGGCCCTC‑TTCGCGGCTAAACCCGCCCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAg  <  2:11208/130‑1 (MQ=255)                                                                                                                                  ccTTCATCCCCAAACCCCCCCCCCCCGGCAATCACCGGGGCCGGGAGGGCGGCTTTTCGTGTTGGCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAg  <  6:172357/135‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GGTAATGCCCATCAAGGCCAGCACCAGGCCAGCGGTGATGCCAAAGAACCAGTGCAATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGT‑GCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTC‑TTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAG  >  KT2440/990568‑990830

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