Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	5,386,506	1 bp→GT	intergenic (‑515/+1136)	CJ019_05070 ← / ← nptII	hypothetical protein/neomycin phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,386,505	1	.	G	90.9%	66.6 / 0.9	22	intergenic (‑514/+1137)	CJ019_05070/nptII	hypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (11/9);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	ACB143	5,386,506	0	C	T	90.5%	44.2 / ‑0.5	22	intergenic (‑515/+1136)	CJ019_05070/nptII	hypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (10/9);  minor base A (0/2);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

TTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CA‑GGATC‑‑A‑‑TC‑CGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  ACB143/5386414‑5386607                                                                                                  ||                                                                                                           ttttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGcccc                                                                 >  8:102260/1‑135 (MQ=255)         aTAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGa                                                          <  8:262590/134‑1 (MQ=255)          tAACAAAAATTGATCGAATTCTTTCATTAAATACCCTAATATCACAGCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATc                                                        <  5:32403/135‑1 (MQ=255)              aaaaaTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCt                                                    >  8:110250/1‑135 (MQ=255)                      aTCGAATTCTTTCATTTACGACCCTAATATGACAGCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGAACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                            <  4:185565/135‑1 (MQ=255)                                  cATTAAAGACCCTAATATGCCATCTCGGCCGGCTGGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGAAC‑‑AACCC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                  <  1:371333/133‑1 (MQ=255)                                   aTTTAAGACCCTAATATCACATCTCGGCCGGCTGGGCGACGACGGCGGTCTCAGG‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                              <  4:366827/136‑1 (MQ=255)                                      tAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGGGT‑‑‑AC‑GGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTg                            <  3:296930/135‑1 (MQ=255)                                               aaTATGACATCTCGGCCGGCTTGTCGACGAAGGCGTTCTCAGT‑‑GGTGTACG‑GTACA‑‑A‑‑ACACCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCATCCTGGGGATTTTTTATTTCTCTGCGGGTGCCAg                  >  4:125691/1‑136 (MQ=255)                                                atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                 <  4:112251/136‑1 (MQ=255)                                                 tatGACATCTCGGCCGGCTTGTCGACGAGGGCGGGCTCTGTGGGGT‑‑‑AC‑GGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                 <  7:82993/135‑1 (MQ=255)                                                        tctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                  <  5:369642/111‑1 (MQ=255)                                                        tctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                  >  6:369642/1‑111 (MQ=255)                                                         ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTAC‑‑AACCC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGTGTGCCAGGGCGTGccc         >  8:154/1‑135 (MQ=255)                                                            ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt       <  4:311429/134‑1 (MQ=255)                                                            ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg     >  4:362330/1‑136 (MQ=255)                                                            ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg     >  4:377128/1‑136 (MQ=255)                                                            ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTg      >  1:199257/1‑135 (MQ=255)                                                             gccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCttcg     >  7:258659/1‑133 (MQ=255)                                                             gccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt       <  5:66238/133‑1 (MQ=255)                                                              ccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTggg    >  2:284182/1‑135 (MQ=255)                                                                ggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAA‑‑CC‑CCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCt  >  2:375146/1‑135 (MQ=255)                                                                                                  ||                                                                                                           TTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CA‑GGATC‑‑A‑‑TC‑CGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  ACB143/5386414‑5386607

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