Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,954,912	G→A	9.7%	E67E (GAG→GAA)	CJ019_05576 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,954,912	0	G	A	9.7%	119.4 / 2.7	41	E67E (GAG→GAA)	CJ019_05576	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTGTATCTGCCTGAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCACCCCCACCTACGTGTACTCGCGCGC  >  KT2440/5954787‑5955043                                                                                                                              |                                                                                                                                    ccGCTGTATCTGCCTGAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAg                                                                                                                           <  8:86214/136‑1 (MQ=255)   gCTGTATCTGCCTGAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAgcg                                                                                                                         >  8:199277/1‑136 (MQ=255)        aTCTGCCTGAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCa                                                                                                                       <  3:205180/133‑1 (MQ=255)                gAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAACAGTCGCCAGCGCAGTAAGGaacc                                                                                                             >  7:104846/1‑133 (MQ=255)                gAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAACAGTCGCCAGCGCAGTAAGGAAAc                                                                                                             >  7:108592/1‑135 (MQ=255)                 aaGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGaaa                                                                                                              >  1:236613/1‑133 (MQ=255)                    ctGGCGATAACGGCACAGCGAACACGTCGCACCAGCACCAGACAAAAGCCACGCCGGCCCCGGCCCAAGAACAGCAGACCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGGAAGGAACCCAGAt                                                                                                        >  8:65398/2‑136 (MQ=255)                    aTGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGa                                                                                                         <  3:142294/135‑1 (MQ=255)                     tGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGAt                                                                                                        <  4:223036/135‑1 (MQ=255)                         gATAACGGCAACGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAAc                                                                                                    >  6:251038/1‑135 (MQ=255)                                    aGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTa                                                                                         >  6:222489/1‑135 (MQ=255)                                      cccccacGGCCCACCACCACCACCCAAACCCCACGCCGGCCCCGGCACAACAACCGCCGCCCGAGCTGCAGTCCGAGCCCGCGCCCGAGCAGGCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACAAcc                                                                                       <  2:236613/131‑1 (MQ=255)                                           acGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACg                                                                                 >  8:36897/1‑136 (MQ=255)                                           acGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACg                                                                                 >  8:255124/1‑136 (MQ=255)                                                   ccagcaccagcCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCt                                                                          >  5:106261/1‑135 (MQ=255)                                                          cagcCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGCCGGTGAGCTGTTcgc                                                                    >  4:141739/1‑134 (MQ=255)                                                             ccAAACGCCAAGCCGGCCCCGGCACAACAACATCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAACAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAg                                                               <  8:108592/136‑1 (MQ=255)                                                             caaaaaCCCAACCCCGCCCCGCAACAACAACAGCAGCCCTAGCTGAAGTCCGAGTCCGAGCCCGAACAGTCCCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACAACCGCGACGGTGAGCTGTTCGCGGAAg                                                               <  8:104846/134‑1 (MQ=255)                                                               aaaaGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGt                                                                                                                     >  3:100587/1‑80 (MQ=255)                                                               aaaaGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGt                                                                                                                     <  4:100587/80‑1 (MQ=255)                                                                 acgccacgccGGCCCCGGAACAAGACCAGCAGCCCGAGCTGCAGTCCGAGCCCGCGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATCACCGCTTTCAACAACCGCGACGGTGAGCTGTTCGCGGAAGGCg                                                            <  6:153763/133‑1 (MQ=255)                                                                    ccaagccGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGa                                                                                                                >  4:108222/1‑80 (MQ=255)                                                                    ccaagccGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGa                                                                                                                <  3:108222/80‑1 (MQ=255)                                                                    ccaagccGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTgg                                                         >  2:178584/1‑135 (MQ=255)                                                                    ccaagccGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGc                                                        >  8:213941/1‑136 (MQ=255)                                                                      aagccGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAACCCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGcc                                                       <  4:42989/135‑1 (MQ=255)                                                                          cGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTcgc                                                                    >  8:142197/1‑118 (MQ=255)                                                                          cGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTcgc                                                                    <  7:142197/118‑1 (MQ=255)                                                                          cGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGt                                                   >  2:79385/1‑135 (MQ=255)                                                                                cGGCACAAGAACAGCAGCCCGCGCTGCAGTCCGAGGCCGAGCCCGAGCAGTCGCCAGCGCACTAAGTAAACCACATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCa                                             <  5:62528/135‑1 (MQ=255)                                                                                  gCACAAGAACAGCAGCCCGAGCGGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCAtt                                           <  3:91563/135‑1 (MQ=255)                                                                                   cacaAGACCCGCAGCCCGAGCTGCAGTCCGTCGCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGCAAACCAGATCAACGCATTCAACTCCCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTg                                          <  2:80252/135‑1 (MQ=255)                                                                                             gccgcCCGAGCGGCAGTCCGCGGCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCATATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTa                                <  1:168515/132‑1 (MQ=255)                                                                                              cagcCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTAc                               >  5:111045/1‑135 (MQ=255)                                                                                                gcCCGAGCTGCAGTCCGAGCCCGAGCCCGAGCAGTCGCCAGCTCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACgg                             <  3:124653/135‑1 (MQ=255)                                                                                                   cGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTc                                                                                              <  6:145448/67‑1 (MQ=255)                                                                                                   cGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTc                                                                                              >  5:145448/1‑67 (MQ=255)                                                                                                   cGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGGCGCCAGCGCAGTACGGAAACCAGATTACCGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCAc                          <  1:40531/135‑1 (MQ=255)                                                                                                      gcggCAGTCCGAGTCCGCGCCCTTGCTGCCGCCTGCGATTTAAGCAAACCTGTTTAACGCTTTCAACTCCCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCAcccc                       <  1:147160/132‑1 (MQ=255)                                                                                                      gCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCAccccc                      >  7:227456/1‑136 (MQ=255)                                                                                                        tGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCACCCCCa                     >  4:25046/1‑135 (MQ=255)                                                                                                                         cccgagcAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCTGCCATTGCTGAACGTTACGGCACCCCCACCTACGTGTACTcgcgc    <  5:222489/135‑1 (MQ=255)                                                                                                                          ccgagcAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACgg                             >  3:221402/1‑109 (MQ=255)                                                                                                                          ccgagcAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACgg                             <  4:221402/109‑1 (MQ=255)                                                                                                                          cccagcAGTCGCCAGCGCAGTAGGGAAACCAGATGAAAGCTTTCAACTACCGCGACCGTGACCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCCTTGCTGAAAGTTCCCGCCCCCCCCCCAACGAGTACTCgcgcgc  >  8:179886/4‑136 (MQ=255)                                                                                                                              |                                                                                                                                    CCGCTGTATCTGCCTGAAGATGGCGATAACGGCAAAGCGCACACGTCGCACCAGCACCAGCCAAAAGCCAAGCCGGCCCCGGCACAAGAACAGCAGCCCGAGCTGCAGTCCGAGTCCGAGCCCGAGCAGTCGCCAGCGCAGTAAGGAAACCAGATGAACGCTTTCAACTACCGCGACGGTGAGCTGTTCGCGGAAGGCGTGGCCCTGTCGGCCATTGCTGAACGTTACGGCACCCCCACCTACGTGTACTCGCGCGC  >  KT2440/5954787‑5955043

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