Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	336,483	A→T	21.0%	P4P (CCT→CCA)	CJ019_00309 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	336,483	0	A	T	21.0%	84.6 / 7.5	52	P4P (CCT→CCA)	CJ019_00309	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/24);  new base T (8/3);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCCTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTC  >  ACB122/336352‑336615                                                                                                                                      |                                                                                                                                      aTGGCCTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAg                                                                                                                                      <  2:247942/133‑1 (MQ=255)     ccgtgttgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTgg                                                                                                                                <  2:364177/132‑1 (MQ=255)     ccTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTgg                                                                                                                                <  4:216567/135‑1 (MQ=255)     ccTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTgg                                                                                                                                <  8:127235/135‑1 (MQ=255)     ccTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTg                                                                                                                                 <  4:46238/134‑1 (MQ=255)     acTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGAGGTTCTGCTTGGCCTTGCTGGCGGCGCGCGTGTCGCGGTTCTGTTT‑GGTTGTTAgggggg                                                                                                                                 >  5:73784/2‑133 (MQ=255)          ttgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATg                                                                                                                           <  5:372599/135‑1 (MQ=255)          ttgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCGTGGCGCGGGTGTG‑TT‑GGTTGGGTGggggggggctt                                                                                                                           >  5:299360/1‑130 (MQ=255)           tgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATgg                                                                                                                          <  5:325135/135‑1 (MQ=255)           tgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATgg                                                                                                                          >  4:81225/1‑135 (MQ=255)           tgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATg                                                                                                                           <  1:161394/134‑1 (MQ=255)           tgtACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGGTTGTG‑TT‑TGGTGGGAGGGGTGGGCCTgg                                                                                                                          >  3:375693/1‑135 (MQ=255)                  cTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCAGGCGGTTCTGCTTGGCCTTGCTGGCGGCGCCCGTTGGGGGTTTCTG‑TTGGGTTGGT‑GGGGGGGCCATgggggtctc                                                                                                                   >  5:70426/1‑130 (MQ=255)                     cTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGGGTGGGGCGGTGCGG‑GT‑GCGTGGGTGGGGTGGGCCTGGGGGAAtgcgt                                                                                                                >  4:217216/1‑131 (MQ=255)                      tCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTg                                                                                                               <  5:195153/135‑1 (MQ=255)                      tCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTg                                                                                                               <  5:17464/135‑1 (MQ=255)                      tCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTg                                                                                                               <  8:55539/135‑1 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGGGTGGCGCGGTTGTG‑GT‑GGTTGGTAGGGGTGGGCCTGGGGACATGCGTgg                                                                                                              >  3:245620/1‑135 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGGGCGGTTGTG‑GT‑GGTTGGTAGGGGTGGGCATGGGGAAATGCGTgg                                                                                                              >  7:181048/1‑135 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAg                                                                                                                                      <  1:14774/111‑1 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAg                                                                                                                                      >  2:14774/1‑111 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTgg                                                                                                              >  2:113422/1‑135 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCGTGGGGCGGTTGTG‑TT‑GGTTGGTTGGGGGGGCGCTGGGGGAATCCGTgg                                                                                                              >  6:372599/1‑135 (MQ=255)                       cGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCGTGGCGCGTTTCTG‑TT‑GGATGGGTTGGGTGGCCCTGGGGTACTCCTTgg                                                                                                              >  8:378166/1‑135 (MQ=255)                         aCGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCTGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTggg                                                                                                             <  3:176763/134‑1 (MQ=255)                               gcgcTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCt                                                                                                      <  5:154211/135‑1 (MQ=255)                                        ccGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTGTG‑TT‑GCGTGGGAGGGGTGGCCCTTGGGCAAGACGTGGGTGCGTCTGTAAGTggg                                                                                             >  1:6546/1‑133 (MQ=255)                                                              ggCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GGTTGGGAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTTGAGCG‑CAATTGCCCg                                                                      >  7:37133/1‑135 (MQ=255)                                                                                ggTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑GT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGgcgcg                                                    >  7:301324/1‑135 (MQ=255)                                                                                ggTTCTGCTTGGCCTTGCTGGCGGCGCGCGTGGCGCGGTTGTT‑GT‑GGGTGGTTGGGGTGGGCATGGGGCAATCCGTGGGGGAGTCGGAAAGTGCGCGGATTGAAGGGCCAATGGCCGCATAAACCAAAccgcgcg                                                     >  6:69016/1‑131 (MQ=255)                                                                                   tCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTGTG‑TT‑TGTTGGTTGGGGTGGGCATGGGGGAATACATTGGTGAGTATTTAAATTCGCGTTTTTAAGCG‑CAAATTCCCCATAAAACCAATCGGCGCGCCt                                                 >  6:85147/1‑135 (MQ=255)                                                                                          ggCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCGTGGTAGGGGTGGCCATGGGGAAATCCGTGGGTGGGTCTGAAAGTGCGCGTATTGCAGCC‑CAATTTCCCGATAAACCCAAGCGGCGCGCCTGGgcgcg                                          >  2:81449/1‑134 (MQ=255)                                                                                           gCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGGGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGgcgc                                           <  1:113422/133‑1 (MQ=255)                                                                                               tGCTGGCGGCGCGCTTGGCGCGGTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGGGAGTCTGAAAGTGCGCGTATTGCAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGc                                     >  7:322049/1‑135 (MQ=255)                                                                                               tGCTGGCGGCGCGCTTGGCGCGGTTCTG‑GT‑GGGTGGTTGGGGTCGCCATGGGGAAATGCTTGGGTGACGCTGAAATTGCGCGTATTGCAGCG‑CAAATGGCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGc                                     >  8:222781/1‑135 (MQ=255)                                                                                                 cTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCAGGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCt                                    <  8:120308/134‑1 (MQ=255)                                                                                                     cggagCGCTTGGCGAGTTTCTG‑GT‑GATTGGTCGGTGTGGCCCTGGGGAAATCCTTTGGTGAGTCTTAACTTGCGCGTATTGAAGCG‑CAATTGACCGATAAACCAAATCGGCTCGCCTGGGCGCTGCAGCTTgaga                               <  5:311717/135‑1 (MQ=255)                                                                                                         gcgcTTGGCGCGTTTATG‑TT‑TGTTTGTTTGGGTGGCCATTGGGAAATCCTTTGGTGAGTCTGAAAATTCGCGTATTTTAGCG‑CAATTGCCCGATAAAACCAACCGGCGCGCCAGGGCGCTGCAGCTTGAGATGtg                           >  4:164664/1‑134 (MQ=255)                                                                                                          cgcTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCGTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCt                          <  6:102424/135‑1 (MQ=255)                                                                                                             ttGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCtgtg                       <  3:81225/135‑1 (MQ=255)                                                                                                                gcggGTTGCTT‑TT‑GCTTGTTTGGGGTGGCGATGGGGAATTCTTTGGGTGAGTGTTAAAGTGGGCTTTTTGAAGCG‑AATTTGCCGTAAAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTggg                     <  3:253608/134‑1 (MQ=255)                                                                                                                gcgcGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGTGAAATCCTTGGGTGAGTCTGAAAGGGCGCTTATTGAAGCG‑AAATTGCCCGATAACCCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTggg                     <  5:85147/134‑1 (MQ=255)                                                                                                                gcgcGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTggg                     <  3:164664/134‑1 (MQ=255)                                                                                                                gcgcGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑AAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTggg                     <  7:281852/134‑1 (MQ=255)                                                                                                                gcgcGTTTCTG‑TT‑GCTTGGTAGGGGGGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTggg                     <  5:146465/134‑1 (MQ=255)                                                                                                                    gTTTCTG‑TT‑GCGTGGTAGGGGGGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑AAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACa                <  7:153300/135‑1 (MQ=255)                                                                                                                    gTTTCTG‑TG‑GCGTGGTTGGGGTGGCCTTGGGTAAATCGTTGGGTGAGTCTGAAAGTGCGCGTATTGCAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACa                <  7:378166/135‑1 (MQ=255)                                                                                                                         tG‑TT‑GCTTGGTCGGGGGGGCCATGGGGAAAGCCTTGGGTGAGTCTGAAAGGGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATcc           <  2:293805/135‑1 (MQ=255)                                                                                                                              ggCTTGGTTGGGGCGGCCATGGGTAAATCGTTGGTTGTGTCTGAAAGTGCGCGTATTGAAGCG‑AAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATcc           <  3:140907/131‑1 (MQ=255)                                                                                                                                   ggTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTga                                 >  3:108939/1‑105 (MQ=255)                                                                                                                                   ggTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTga                                 <  4:108939/105‑1 (MQ=255)                                                                                                                                     tAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTc  >  1:313266/1‑134 (MQ=255)                                                                                                                                      |                                                                                                                                      ATGGCCTTGTTGTACACCTGCTCGACGCTGGCGCTGTCGCCGGCACTGGCCTTGACCGCTTGGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTG‑TT‑GCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCG‑CAATTGCCCGATAAACCCAAGCGGCGCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTC  >  ACB122/336352‑336615

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