Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	1,063,986	T→G	27.2%	D151A (GAC→GCC)	CJ019_01001 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	1,063,986	0	T	G	27.2%	41.0 / 2.3	40	D151A (GAC→GCC)	CJ019_01001	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/24);  new base G (9/3);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 3.93e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCAGGTCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAGCGCATACAC  >  KT2440/1063855‑1064116                                                                                                                                      |                                                                                                                                      ggCCCAGGTCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTAGCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCg                                                                                                                                    <  6:140678/135‑1 (MQ=255) ggCCCAGGGCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGGGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCg                                                                                                                                    <  6:172039/135‑1 (MQ=255) ggCCCAGGGCGGTGGCCACCTCGATGGCAGCTTCCAGAAGCAAGGGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCg                                                                                                                                    <  2:23026/135‑1 (MQ=255) ggCCCAGGGCGGTGCCCACCTCGATGGCAGCTGCCAGTCGCACGGGGCCCGCCAGGCGCAGCCGGCCAAGCCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCg                                                                                                                                    <  6:328226/135‑1 (MQ=255)   ccccGGTCGGTGGCCACCTCGTGGGCAGCTCCAATTAGCAAGTGGCCCGCCAGGCGCCGCAGCCCACGGCCGACCAGTAGCTGCAGCAGCGTCCCTGACCAGCCGCGCCCGG‑CCCCGGACCTGGC‑GGCGGCCg                                                                                                                                    <  8:301647/133‑1 (MQ=255)   cccTGGTCGGTGGCCACCTGGATGGCAGCTTCCAGTCGCTCGGGGCCCGCCAGGCTCAGCCGGCCATCGCCGACCAGTCTGCCCCGCCGCGTCCGTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCggg                                                                                                                                  <  4:260429/135‑1 (MQ=255)           ggTGGCCACCTCGATGTCAGCTTCCAGTCGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCACAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTa                                                                                                                          <  8:12143/135‑1 (MQ=255)             tGGCCACCGCGATGGCAGCTTCCAGTAGCACGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAgg                                                                                                                        <  2:276678/135‑1 (MQ=255)               gCCCCCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCACGGCCGCCCAGTATCTGCAGCAGCGTACCTGCCCAGAAGCGCCCGG‑CCCCCCCCCCGGC‑GGCGGCCCCGGGAGGGGCGgtg                                                                                                                      >  7:9870/1‑135 (MQ=255)               gCCACCTCGATGGCAGCTTCCCGTAGCAAGTGGCCCGCCAGGCGCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCACAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGGCCGGGGCAGGGTAGgtg                                                                                                                      <  1:215560/135‑1 (MQ=255)               gCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGCAGCGCCCGGCCCCCGGCCCGGGC‑GGGGGCcggggcgggggggtg                                                                                                                       >  2:319592/1‑130 (MQ=255)                          ggCAGCTTCCAGTAGCACGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGACGCGCCCGG‑CCCCGGACCTGGC‑GGCGGCCGGGGCAGGGTAGGTGCGTGCGTAATg                                                                                                           <  5:39533/135‑1 (MQ=255)                          ggCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATg                                                                                                           >  2:12683/1‑135 (MQ=255)                               cTTCCAGTAGCAAGGGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGACGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGcgcg                                                                                                      <  4:300205/135‑1 (MQ=255)                                 tCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCAGCG‑GCCGGCCGGGG‑AGGGTTGGTGCGTGGGTAATGGcgcgcc                                                                                                    >  4:67651/1‑132 (MQ=255)                                                       ggCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGGCCGGGC‑GGCGGGCGGGGGCGGGGTGGTGGGTGCGTTATTGCGCGGCGGATGCCCG‑CG‑CATTTCC‑acg                                                                            >  8:91140/1‑132 (MQ=255)                                                                  ccAAGGCCGACCAGTATCTGCAGCCGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGc                                                                <  8:40206/135‑1 (MQ=255)                                                                  ccAAGGCCGACCAGTAGCTGCCCCAGCGTACGTCCCCCCACGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGc                                                                <  7:70273/135‑1 (MQ=255)                                                                       gCCGACCAGTAGCTCCAGCAGCGTACCTGCCCAGACGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAg                                                           <  3:57347/135‑1 (MQ=255)                                                                         cGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCAGGCGGGGGGCCGGGGCGGGTTTGGTGGGGGGGTATTGG‑GCGACTGATTGCCG‑AG‑CAATTGC‑AGGTAACCGGCCACccccccc                                                          >  4:60455/1‑131 (MQ=255)                                                                              aGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGTCCGGGG‑GGCGGGCGGGGGCGGGGTGGGGGGGGGGGAATGGGGCGGCGGATGGCAG‑AGCCTTTGCC‑GGGCAACCGCccccccccgcaaaa                                                       >  7:27559/1‑127 (MQ=255)                                                                              aGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCGGGG‑GGCGGCCCGGGCGCGGGAGGGGCGGGCGGTTTGGCGCGGCAGCAGGCAGCGG‑CATTCCA‑CGGTACCCGGCC‑GCCCCCGCAAAAGc                                                    >  4:54164/1‑135 (MQ=255)                                                                              aGTAGCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGc                                                    <  7:57158/135‑1 (MQ=255)                                                                               gtagcggcagcaGCGCACCTGCCCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCa                                                   <  3:147203/129‑1 (MQ=255)                                                                               gTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCa                                                   <  6:250017/135‑1 (MQ=255)                                                                                tctCTGCAGCAGCGTACCTTCGCAGACGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGAGCGTGCGTAATCGCGCGACTGATGCGAG‑AT‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAg                                                  <  8:132472/133‑1 (MQ=255)                                                                                   cTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCAGGC‑GGCGGCCGGGGCAAGGGAGGTGGGTGCGTAAAGGGGGGCCTGATGCCAC‑AG‑CATTGCCAAGGAAACCGGCC‑GCCCCCCCAATAGCACTgg                                               >  2:59662/1‑134 (MQ=255)                                                                                     gcagcaGCGTACCTGACCAGAAGCGCCCGG‑CCCCCGCCCTGGC‑GGCGGCCGGGGGAGGGTAGGTGCGTGCGTAATTGCGCGCCTGATGCCAG‑AG‑CATTTCC‑AGGTAACCCGCCCGCCCCCGCAAAAGCAGTGccc                                             >  1:178088/1‑135 (MQ=255)                                                                                       agcagcGTACCTGCCCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGccc                                             <  6:71517/133‑1 (MQ=255)                                                                                          ccggaaCCTGACCCGACGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCa                                        <  6:184982/130‑1 (MQ=255)                                                                                            cGTACCTGACCAGAAGCGCCCGG‑CCCCGGCCCAGGC‑GGCGTCCGGGGCAGGGTATGTGCCTGCGTAATTGCGCGACTGATGCCAG‑AG‑CAATGCC‑AGGTAACCGGCCAGCCCCCGCAATAGCAGTGccccatacac                                      >  6:163164/1‑129 (MQ=255)                                                                                               aCCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCt                                   <  6:14501/135‑1 (MQ=255)                                                                                                 cTGACCAGAAGCGCCCGG‑CCCCGGCCCAGGC‑GGGGTGCCGGGCGGGGGTGGGGGGGGGGTATTGGCGCGCCTGATGCCAGCAG‑AATT‑CC‑AGGGAACCCGCCAGCCCCCCCAAAAGCAGTGGCCCATCAACCCGGc                                 >  6:147403/1‑135 (MQ=255)                                                                                                                 gg‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGc                 <  3:315523/135‑1 (MQ=255)                                                                                                                  g‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCg                <  7:20137/135‑1 (MQ=255)                                                                                                                  g‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCg                <  2:25547/135‑1 (MQ=255)                                                                                                                  g‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGGAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCg                <  3:54164/135‑1 (MQ=255)                                                                                                                    ccccGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCgg               <  7:91140/135‑1 (MQ=255)                                                                                                                        ggACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAG‑CCCAGCAATAGCAGAGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAgc          <  1:48373/135‑1 (MQ=255)                                                                                                                                  ggcgTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGCCTGATGCCAG‑AG‑CATTGAC‑AGGTAACCGGCCAGCCCCCGCAATAGCAGTGCCCCGTCAACCCTGCCCCGGTGGGCCTTGGCGGCCAgcgccgaccc  >  2:153525/1‑129 (MQ=255)                                                                                                                                      |                                                                                                                                      GGCCCAGGTCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGG‑CCCCGGACCTGGC‑GGCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAG‑AG‑CATTGCC‑AGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAGCGCATACAC  >  KT2440/1063855‑1064116

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: ‑