Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,184,142	A→T	9.2%	I52N (ATC→AAC)	CJ019_03932 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,184,142	0	A	T	9.2%	152.9 / 5.6	52	I52N (ATC→AAC)	CJ019_03932	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base T (2/3);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

AGTGGTGCAAGACTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGTGCTCTGGATGATGCGATCGAGCT  >  ACB143/4184010‑4184274                                                                                                                                     |                                                                                                                                     aGTGGTGCAAGACTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGAtct                                                                                                                                   <  7:81843/136‑1 (MQ=255)        cAAGACTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAActa                                                                                                                            >  1:341409/1‑135 (MQ=255)         aaGACTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAActct                                                                                                                           >  2:127866/1‑136 (MQ=255)             cTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCg                                                                                                                       >  1:161889/1‑136 (MQ=255)                cATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTAACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCt                                                                                                                    <  4:302434/136‑1 (MQ=255)                 ataataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtt                                                                                                                   >  1:222728/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCATGATCTCTAACTCTTTCGCCtttt                                                                                                                 >  4:395963/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 >  5:30709/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 >  3:235168/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 >  7:27273/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 <  2:50743/136‑1 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 >  2:31108/1‑136 (MQ=255)                   aataTTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCcttt                                                                                                                 >  5:360759/1‑136 (MQ=255)                    atatTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCttttt                                                                                                                <  4:374148/136‑1 (MQ=255)                    atatTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 <  1:348360/135‑1 (MQ=255)                    atatTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 <  5:53704/135‑1 (MQ=255)                    atatTTTTCTCCGACGAAGCTCACGAAGCCGGAGACACTTGTCTCTCACATTAAAATAACGGGAACTGAATGTTTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCtttt                                                                                                                 <  6:21135/135‑1 (MQ=255)                     tatTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTATCGCCTTTTTc                                                                                                               >  7:351278/1‑136 (MQ=255)                             tcCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTa                                                                                                       >  5:140359/1‑136 (MQ=255)                             tcCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTa                                                                                                       >  7:84451/1‑136 (MQ=255)                             tcCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTa                                                                                                       >  2:198774/1‑136 (MQ=255)                             tcCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGAACTCTAACTCTTTCGCCTTTTTCGTGGTGTa                                                                                                       >  3:5206/1‑136 (MQ=255)                             tcCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGAGTTAATTTTGCGGATCATCAACTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTGTTTAGTGGTGTa                                                                                                       >  3:10474/1‑136 (MQ=255)                                 ccgAACCTCACGAAGCCGGAGACACTGGTCTCGGACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAATCTCGAGCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTAcgag                                                                                                   <  1:98164/136‑1 (MQ=255)                                  cgAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTAcgagt                                                                                                  >  2:402381/1‑136 (MQ=255)                                   gAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTAcgagtc                                                                                                 <  2:161889/136‑1 (MQ=255)                                   gAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTAcgagtc                                                                                                 <  2:272313/136‑1 (MQ=255)                                                       cacTGGTCTCGCACTTTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCt                                                                              >  5:126458/1‑135 (MQ=255)                                                           ggTCTCGCACATTAACATAATGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAAc                                                                         >  6:169238/1‑136 (MQ=255)                                                            gTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAAc                                                                         <  4:81054/135‑1 (MQ=255)                                                            gTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACt                                                                        <  2:211522/136‑1 (MQ=255)                                                             tctcGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTc                                                                       >  7:197649/1‑136 (MQ=255)                                                                                  tCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCtttt                                                  >  5:43619/1‑136 (MQ=255)                                                                                    aCTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCttttt                                                 <  8:193611/135‑1 (MQ=255)                                                                                    aCTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTc                                                <  1:138364/136‑1 (MQ=255)                                                                                        aaTGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGt                                            >  2:216106/1‑136 (MQ=255)                                                                                            aTTTAATTTTTCGGAGTATCACCTTGTAAAAATCGATCAAGATCTCTTACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCtt                                         <  6:141005/135‑1 (MQ=255)                                                                                             tttAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTg                                        <  7:252440/135‑1 (MQ=255)                                                                                               tAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGTTCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTg                                        <  1:385494/133‑1 (MQ=255)                                                                                               tAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGTTCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTg                                        >  2:390831/1‑133 (MQ=255)                                                                                               tAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGTTCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTg                                        <  1:390831/133‑1 (MQ=255)                                                                                               tAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGTTCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTg                                        >  2:385494/1‑133 (MQ=255)                                                                                                aaTTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTaa                                     <  4:235168/135‑1 (MQ=255)                                                                                                  ttttGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGt                                  <  4:250669/136‑1 (MQ=255)                                                                                                       cGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTATCTCTTTCGCCTTTTTCGTGGTGTACTAGTCGATTTTAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTagaga                             <  7:346500/136‑1 (MQ=255)                                                                                                       cGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTagaga                             <  6:43619/136‑1 (MQ=255)                                                                                                       cGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTagag                              <  2:372966/135‑1 (MQ=255)                                                                                                          agaATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCg                          >  2:243818/1‑136 (MQ=255)                                                                                                           gaATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGt                         <  2:222728/136‑1 (MQ=255)                                                                                                                cccTTGTAAAAATCGATCAAGTTCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGTGctct                    <  5:162253/136‑1 (MQ=255)                                                                                                                               aTCAAGATCTCTAACTCTTTCGCCTTTTTCGTCGGGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGTGCTCTGGATGATGCGATCGa     <  3:100269/136‑1 (MQ=255)                                                                                                                                  aaGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGTGCTCTGGATGATGCGATCGAGCt  <  6:159419/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     AGTGGTGCAAGACTCCATAATATTTTTCTCCGCCGAAGCTCACGAAGCCGGAGACACTGGTCTCGCACATTAACATAACGGTCACTGAATGATTTAATTTTGCGGAGAATCCCCTTGTAAAAATCGATCAAGATCTCTAACTCTTTCGCCTTTTTCGTGGTGTACGAGTCGAGTTGAGTATCTGCGCCTAAAACTCGCATCTCTGAGACTTGCTTTTTCAGGTCTTGTAAAGTAGAGAGCGTGCTCTGGATGATGCGATCGAGCT  >  ACB143/4184010‑4184274

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