Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	5,124,825	A→G	85.1%	intergenic (‑103/‑141)	recQ ← / → CJ019_04820	ATP‑dependent DNA helicase RecQ/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,124,825	0	A	G	85.1%	126.4 / 12.7	47	intergenic (‑103/‑141)	recQ/CJ019_04820	ATP‑dependent DNA helicase RecQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (17/23);  total (20/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 = 8.12e-01

ACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGG  >  ACB143/5124700‑5124955                                                                                                                              |                                                                                                                                   aCGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTg                                                                                                                          >  4:2124/1‑136 (MQ=255) aCGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCAGTACCATCCAGGATAAAGGCTAGACGAAAAACACTCCACGAACACGCCAACTGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCATg                                                                                                                          >  6:195918/1‑136 (MQ=255)          gCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGaaaaa                                                                                                                 <  5:124211/136‑1 (MQ=255)          gCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGaaaaa                                                                                                                 <  2:134940/136‑1 (MQ=255)            cTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATc                                                                                                               <  5:39583/136‑1 (MQ=255)            cTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATc                                                                                                               <  5:139005/136‑1 (MQ=255)                  gAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGa                                                                                                          <  2:218257/135‑1 (MQ=255)                  gAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGa                                                                                                          <  3:213546/135‑1 (MQ=255)                      aTGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGtata                                                                                                     <  6:13819/136‑1 (MQ=255)                      aTGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGtata                                                                                                     >  8:143753/1‑136 (MQ=255)                           cAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCgca                                                                                                <  6:115626/136‑1 (MQ=255)                           cAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCgca                                                                                                <  3:217410/136‑1 (MQ=255)                                 cAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCg                                                                                          <  5:127552/136‑1 (MQ=255)                                 cAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCg                                                                                          >  4:53880/1‑136 (MQ=255)                                      cGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                     >  7:190403/1‑136 (MQ=255)                                      cGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                     >  5:136766/1‑136 (MQ=255)                                      cGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGtgctg                                                                                     >  5:115921/1‑136 (MQ=255)                                               aTCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCAc                                                                            <  6:204006/136‑1 (MQ=255)                                                      aTAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCgcggc                                                                      >  7:12726/1‑135 (MQ=255)                                                       tAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGa                                                                    <  2:290507/136‑1 (MQ=255)                                                        aaaGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGa                                                                    >  1:240097/1‑135 (MQ=255)                                                        aaaGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGa                                                                    >  1:286598/1‑135 (MQ=255)                                                          aGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAg                                                                 <  7:143753/136‑1 (MQ=255)                                                             cTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATc                                                                                                               <  6:273616/87‑1 (MQ=255)                                                             cTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATc                                                                                                               >  5:273616/1‑87 (MQ=255)                                                                  cGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTccc                                                         >  4:245481/1‑136 (MQ=255)                                                                  cGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTcc                                                          >  6:16048/1‑135 (MQ=255)                                                                              cccgcACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctg                                             <  4:219295/136‑1 (MQ=255)                                                                                   acacGCCAACTGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctgccgct                                        <  8:61645/136‑1 (MQ=255)                                                                                   acacGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctgccgct                                        <  8:190403/136‑1 (MQ=255)                                                                                     acGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCccgctgccgctg                                       >  1:143548/1‑135 (MQ=255)                                                                                               cGTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGc                               >  2:96044/1‑133 (MQ=255)                                                                                                gTAAGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGa                           <  6:274180/136‑1 (MQ=255)                                                                                                  aaGCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAAt                         >  1:243636/1‑136 (MQ=255)                                                                                                    gCGGTTTTTGCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTc                       <  7:234479/136‑1 (MQ=255)                                                                                                    gCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATc                                                                                                               <  4:260733/48‑1 (MQ=255)                                                                                                    gCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATc                                                                                                               >  3:260733/1‑48 (MQ=255)                                                                                                    gCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCATGGATGAAAAATCGCGGAGTATACCGCAGCAACGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTc                       >  6:109618/1‑136 (MQ=255)                                                                                                            tGCCAGGCAGCCATTGAAGCCTGCCCCTGGATG‑AAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGt               <  5:219127/135‑1 (MQ=255)                                                                                                             gCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGt               <  1:11539/135‑1 (MQ=255)                                                                                                             gCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGt               <  1:9871/135‑1 (MQ=255)                                                                                                             gCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGt               <  6:272034/135‑1 (MQ=255)                                                                                                             gCCAGGCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGt               <  5:16048/135‑1 (MQ=255)                                                                                                                 ggCAGCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTAtt          <  5:36692/136‑1 (MQ=255)                                                                                                                     gCCATTGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGTTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTccc       >  6:283611/1‑135 (MQ=255)                                                                                                                         ttGAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAgg  <  6:127330/136‑1 (MQ=255)                                                                                                                           gAAGCCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAgg  >  2:237090/1‑134 (MQ=255)                                                                                                                              |                                                                                                                                   ACGCGCTGAGCCTGTTCGAGCATGTGCAACTCCAAATCGCCGTACCATCCTGGATAAAGGCTGGACGAAAAACCCGCCCCGCACACGCCAAATGCGTAAGCGGTTTTTGCCAGGCAGCCATTGAAACCTGCCCCTGGATGAAAAATCGCGGAGTATACCGCAGCACCGTGCTGCACATGCACCGCGGCGATAGACGTTCCCTTTGCCCCGCTGCCGCTGCAAGGTGCTAGAATTCGTTATCGTTTATTCCCCAAGG  >  ACB143/5124700‑5124955

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