Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	36,268	T→C	6.7%	D72G (GAC→GGC)	CJ019_00044 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	36,268	0	T	C	6.7%	186.3 / 4.1	60	D72G (GAC→GGC)	CJ019_00044	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/22);  new base C (2/2);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

TGTTCTGTTCAAATCGTTTACTACTTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGA  >  ACB122/36139‑36399                                                                                                                                  |                                                                                                                                    tgttctgttcAAATCGTTTACTACTTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGc                                                                                                                                <  5:187695/135‑1 (MQ=255)       gttcAAATCGTTTACTACTTGTTGTTGATGCATTGCGAATACATCGAGCAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAATGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACtca                                                                                                                          >  4:179369/1‑133 (MQ=255)        ttcAAATCGTTTACTACTTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCGTGCGACTGAc                                                                                                                         >  8:142639/1‑135 (MQ=255)                  ttACTACTTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAAttt                                                                                                               <  4:95216/135‑1 (MQ=255)                    actactTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTc                                                                                                              <  3:166280/134‑1 (MQ=255)                              gttgATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTAt                                                                                                   >  3:107825/1‑135 (MQ=255)                                tgATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCGTGCGACTGACGTACGAATTTCGCAACGTCTATCt                                                                                                 <  7:142639/135‑1 (MQ=255)                                  aTGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCtt                                                                                                <  8:76032/134‑1 (MQ=255)                                         gCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAATGTCATCAACCGTGGATAGACCTGGAAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCat                                                                                        >  1:245307/1‑135 (MQ=255)                                            aaTACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATAtt                                                                                     <  7:54873/135‑1 (MQ=255)                                              tACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTg                                                                                   <  2:192270/135‑1 (MQ=255)                                                       tAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGa                                                                          >  7:52249/1‑135 (MQ=255)                                                            aaCAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCAATATTTGTAATCCTGAAtctc                                                                     >  7:53186/1‑135 (MQ=255)                                                            aaCAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGCCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAAtctc                                                                     >  3:141835/1‑135 (MQ=255)                                                                  cTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCAt                                                               >  6:267706/1‑135 (MQ=255)                                                                    tGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTc                                                                                                              >  6:330940/1‑86 (MQ=255)                                                                    tGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTc                                                                                                              <  5:330940/86‑1 (MQ=255)                                                                         gTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTg                                                        >  1:374408/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCACCCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCAGGAATCTCTGGCATTGGCTTGCGAtt                                                   >  2:226473/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCg                                                                                                             >  2:333303/1‑77 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCg                                                                                                             <  1:333303/77‑1 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  1:2801/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  6:83458/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  6:47238/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  6:33015/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  2:231181/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  2:44159/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  3:266314/1‑135 (MQ=255)                                                                              cTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCACGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCAtt                                                   >  7:274182/1‑135 (MQ=255)                                                                                   aGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTtct                                              <  4:142671/135‑1 (MQ=255)                                                                                        tCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCcgcg                                         <  7:381764/135‑1 (MQ=255)                                                                                          gAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTc                                                                                               <  2:78853/79‑1 (MQ=255)                                                                                          gAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTc                                                                                               >  1:78853/1‑79 (MQ=255)                                                                                          gAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTa                                       >  6:73368/1‑135 (MQ=255)                                                                                          gAAGGTCATCAACCGTGGATAGACATGGCAGGTACGGGGGTCCTGCGAATGACCTCCGAATTTCGCAACGTCTATCTTCATGCCCTATTTGTAATCATGAATATCTGGCATTGGCTTGCCATTTTTCTCCGCGt                                        >  4:350876/1‑134 (MQ=255)                                                                                             ggTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACtt                                    <  2:282665/135‑1 (MQ=255)                                                                                                    aaCCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCcgcg                                         <  4:340201/123‑1 (MQ=255)                                                                                                    aaCCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCcgcg                                         >  3:340201/1‑123 (MQ=255)                                                                                                    aaCCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACt                             >  5:248072/1‑135 (MQ=255)                                                                                                    aaCCGTGGATAGACCTGGCAGGTACGGGGGCCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCa                               >  5:51150/1‑133 (MQ=255)                                                                                                    aaCCGTGGATAGACCTGGCAGGTACGGGGGCCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCa                               <  6:51150/133‑1 (MQ=255)                                                                                                      ccGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGc                           <  5:47238/135‑1 (MQ=255)                                                                                                            aTAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCt                      <  2:2801/134‑1 (MQ=255)                                                                                                            aTAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTg                     >  8:315474/1‑135 (MQ=255)                                                                                                            aTAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTg                     >  5:186096/1‑135 (MQ=255)                                                                                                            aTAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTg                     >  4:144086/1‑135 (MQ=255)                                                                                                            aTAGACCTGGCAGGTACGGGGGTCATGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCACTGGAATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTg                     >  6:22761/1‑135 (MQ=255)                                                                                                                aCCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCTCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGgtgt                 <  8:362037/135‑1 (MQ=255)                                                                                                                 ccTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGtt                <  4:107825/135‑1 (MQ=255)                                                                                                                 ccTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGtt                <  3:17712/135‑1 (MQ=255)                                                                                                                 ccTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGtt                <  5:83458/135‑1 (MQ=255)                                                                                                                   tGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTcc              <  8:52249/135‑1 (MQ=255)                                                                                                                   tGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTCGTGTTcc              >  2:340494/1‑135 (MQ=255)                                                                                                                    ggCAGGTACGGGGGTCCTGCGACTGCCGTACGAATTTCGCAACGTCTATCGTCATGCCATATTTGTAATCATGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGGTCCt             >  6:229670/1‑135 (MQ=255)                                                                                                                        ggTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTcc              >  5:219242/1‑130 (MQ=255)                                                                                                                        ggTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTcc              <  6:219242/130‑1 (MQ=255)                                                                                                                        ggTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGgt         >  8:80670/1‑135 (MQ=255)                                                                                                                        ggTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGgt         >  7:353158/1‑135 (MQ=255)                                                                                                                             gggggTCCTGCGACTGACGTACGACGTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCAGTGGCGTGCCATTTGTCTCCGCGTCCTTGGCCACTGCGGGCGGGTGTGCCGGGgtggctt    >  3:216977/1‑131 (MQ=255)                                                                                                                               gggCCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGa  <  4:141835/135‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    TGTTCTGTTCAAATCGTTTACTACTTGTTGTTGATGCATTGCGAATACATCGAGTAACGAACAGTCTTGAGCGTGCCCTGGGAGTCCTCGAAGGTCATCAACCGTGGATAGACCTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGA  >  ACB122/36139‑36399

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