Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,937,715	C→T	6.2%	R109H (CGC→CAC)	CJ019_03722 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,937,715	0	C	T	6.2%	186.4 / 3.6	66	R109H (CGC→CAC)	CJ019_03722	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (2/2);  total (35/31)
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.20e-01

AGACCAGGTTCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTGGGAACAGGTCG  >  ACB122/3937584‑3937846                                                                                                                                    |                                                                                                                                    aGACCAGGTTCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGc                                                                                                                                  >  5:205671/1‑135 (MQ=255) aGACCAGGTTCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTGTTTCGCCGCTCTGCAGGAAGCCCTGCAGGTCGTCGCGGTCACCCCTGCggg                                                                                                                                  >  5:148406/1‑134 (MQ=255)      aGGTTCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCCCCCATGCGGCCGAt                                                                                                                              >  3:358241/1‑134 (MQ=255)        gTTCAGCTCGTTGGCCGATACCCGATCGAGGCGAAACGTGTCTTCGATCACGTCGCCTAGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGATGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGaa                                                                                                                           >  7:56329/1‑135 (MQ=255)          tCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGTGGCCGATGAACt                                                                                                                         >  5:403622/1‑135 (MQ=255)          tCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGTGGCCGATGAACt                                                                                                                         >  5:152797/1‑135 (MQ=255)           cAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCCCCCATGCGGCCGATGAACt                                                                                                                         >  1:92181/1‑134 (MQ=255)           cAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACt                                                                                                                         >  3:54880/1‑134 (MQ=255)              cTCGTGGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTGGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGaa                                                                                                                     <  1:314246/135‑1 (MQ=255)               tccgtGGCCGTTACGCGATCGAGGCGATAGGTGTCTTCGATCACGGCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGaaa                                                                                                                    <  2:188113/131‑1 (MQ=255)               tCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGaaa                                                                                                                    <  3:185776/135‑1 (MQ=255)                         aTACCCGATCGAGGCGATAGGTGTCTTCGATCACGGCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCg                                                                                                          <  6:205249/135‑1 (MQ=255)                           aCCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCggg                                                                                                        >  6:365111/1‑135 (MQ=255)                           aCCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCAACCATGCGGCCGCTGAACTGGAAAGGCAGTTCCggg                                                                                                        >  7:423798/1‑135 (MQ=255)                              cGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCCGGGcc                                                                                                     >  4:463414/1‑135 (MQ=255)                                   gAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCgg                                                                                                >  6:429378/1‑135 (MQ=255)                                   gAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCgg                                                                                                >  6:463695/1‑135 (MQ=255)                                       cGATAGGGGTCTTCGATCACGTCGCCTGGGCGCACGACGTAGAGTTGAGCGCCGCTCTCCAGCAAGCCCAGCAGGTCGTCGCGGCCACCCAGGCGGCCGATCAACTGGAAAGGCAGTTCCGGGGCCGGCGGGGcc                                                                                            <  8:90553/135‑1 (MQ=255)                                                  ttCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGcgacg                                                                                   >  7:29543/1‑133 (MQ=255)                                                   tCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGaaa                                                                                                                    <  4:149278/99‑1 (MQ=255)                                                   tCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGaaa                                                                                                                    >  3:149278/1‑99 (MQ=255)                                                   gcGATCACGTCGCCTGGGCGCACGACGTAGAGTTTTTCGCCGCGCTGCAGGACGACCTGCAGGGCGGCGCGGGCCCCCCTGCGGCCGATGAACTGGAACGGCCGGGCCGGGGCCGGCGGGGCCGCGGCGACGACa                                                                                <  4:54880/134‑1 (MQ=255)                                                     gATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAgg                                                                              >  7:139889/1‑135 (MQ=255)                                                         aCGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCg                                                                          >  3:119192/1‑135 (MQ=255)                                                             cGCCTTGGCGCAGGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGGCGCGGCGACGACAGGGGCGGtca                                                                      >  7:118542/1‑135 (MQ=255)                                                             cGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGGCGCCAGGTGCGGtca                                                                      >  1:233090/1‑135 (MQ=255)                                                             cGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGAAGAACTGGAAAGGCAGTTAACGGGCCGTCGGGGCCGCCGGGGCGACAGGTGGGGtga                                                                      >  5:66581/1‑133 (MQ=255)                                                              ccgttGGCGAAGTACGTAGTTTTTTTCGCGGCTCTGCTGGAAGACGGGCTGGTCGTCGCGGTCACCCATGCGGCCTATTAACTTGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGtcat                                                                     <  4:270543/132‑1 (MQ=255)                                                               ccTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGGCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGtcat                                                                     <  6:178456/134‑1 (MQ=255)                                                                 tgggCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGGCACCCATGCGTCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCAc                                                                  <  2:304866/133‑1 (MQ=255)                                                                    gcgcACGACGTATATTTTTGCGCCGCTCTGCAGGAAGACCTGCAGGGCGGCGCGGTCACCCATGCGGCCGATGAACTGTAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACg                                                                 <  3:346245/133‑1 (MQ=255)                                                                             gTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCgg                                                       >  3:100914/1‑134 (MQ=255)                                                                              tAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCAGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTg                                                     <  3:461151/135‑1 (MQ=255)                                                                              tAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCAGGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTg                                                     <  7:380719/135‑1 (MQ=255)                                                                                      ttCCCCGCTCTGCAGTAAGACCTGCAGGTCGGCGCGGACACCCATGCGGCCTATTAACTGGAACGGCAGTTCGGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTgg                                              <  5:365111/134‑1 (MQ=255)                                                                                      ttCCCCGCCCGGCAGGAACACCTGCTGGTCGTCGCGGTCACCCATGGGGCCCATGAACTGGAAACGCAGTTCCGGGGAAGTAGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTGCGGTGACAGTgg                                              <  6:435565/134‑1 (MQ=255)                                                                                         gccgcTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCg                                                                                                    <  6:259404/77‑1 (MQ=255)                                                                                         gccgcTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCg                                                                                                    >  5:259404/1‑77 (MQ=255)                                                                                         gccgcTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTTCGGTCATCACGGGGTGTTAGGTGACCGTGGCCAg                                          >  2:107454/1‑135 (MQ=255)                                                                                         gccgcTCTGCAGGAAGACCTGCAGCTCGTCGCGGTCACCCATGCGGCCGATGCACAGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGGTCGGTGACCGTGGCCAg                                          >  7:274023/1‑135 (MQ=255)                                                                                          ccgcTCTGCAGGAAGACCTGCAGGGCGTCGCGGTCACCCATGCGGCCGATGAACTGGAACGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAgg                                         <  6:205671/135‑1 (MQ=255)                                                                                                 gctggAAGACCTGATGGTGGTCGCGGTCACCCATGCGGCCTTTGAAGTGGAAGGGCAGTTCGGGGGCGGTTGGGGCGGGGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTttgg                                  <  4:200957/132‑1 (MQ=255)                                                                                                     gAAGACCTGCAGGTCGTCGCGGTCACCCATGTGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTttgggttg                              <  6:403622/135‑1 (MQ=255)                                                                                                      acgACCGGCTGGGCGTCGCGGTCACCCAGGCGGCCGATCAACTGGAAAGGCTGTGCGGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTttgggttgg                             <  5:371187/133‑1 (MQ=255)                                                                                                      aaGCCCTGCTGGGCGGCGCGGCCACCCTTGCGGCCGATGAACTGGAAAGGCCGTGCCGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTttgggttgg                             <  5:69427/135‑1 (MQ=255)                                                                                                      aaGACCTGCAGGTCGTCGCGGTCACCCATGTGGCCGATGAACTGGAAAGGCAGTGCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTttgggttg                              <  6:152797/134‑1 (MQ=255)                                                                                                        gACCTGCTGGTCGGCGCGTTCAGCCATGCGGCCGATGAACTGGAACGGGCGTTCGGGGGCGGTCGGGGCGGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTtgggttgggt                           <  8:423798/135‑1 (MQ=255)                                                                                                          ccTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGCCAGGTGCGGTCATCACGGGTTGTTCGGTGCCCGTGGCCAGGGTTTTGGGTTGGGTcc                         >  1:13498/1‑135 (MQ=255)                                                                                                           cTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGGCATCGCGGGTTGTTCGGTGCCCGTGGCCAGGGTTTTGGGTTGGGTcct                        >  3:105767/1‑134 (MQ=255)                                                                                                             gctggtggtCGCGGGCACCCATGCGGCGTTTCAAGTGTAATGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCAtt                      <  4:235011/132‑1 (MQ=255)                                                                                                             gCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCAtt                      >  7:447518/1‑135 (MQ=255)                                                                                                             gCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCAtt                      <  8:29543/135‑1 (MQ=255)                                                                                                               tgggggtcgCGGTCCCGCAGCCGCCGGATGAAGTGGAACTGCAGGTGCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTGCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCAtt                      <  4:324868/128‑1 (MQ=255)                                                                                                               aGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGc                    >  7:53879/1‑135 (MQ=255)                                                                                                               aGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGGCGGGGCCGCGGCGCCGCCAGGTGCGGTCCTCACGGGGTGTTGGGGGACCGGGGCCAGGGTTTTGGGGTGGGTGCCTTGc                    >  7:189168/1‑135 (MQ=255)                                                                                                               aGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGGCGGGGCCGCGGCGCCGACAGGTGCGGTCATCACGGGGTGTTGGGTGACCGGGGGCAGGGGTTTGGGTTGGGTCCGTTGc                    >  5:37507/1‑135 (MQ=255)                                                                                                                 gtcgtcGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTg                  >  5:162539/1‑135 (MQ=255)                                                                                                                      cgcgGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTg             >  1:19890/1‑135 (MQ=255)                                                                                                                       gcgGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGCCGCCAGGTGCGGTCATCACGGGTTGTTGCGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTgg            >  5:167583/1‑135 (MQ=255)                                                                                                                        cgGTCCCCCCGGCGGCCGAGGAACTGGACAGGCCGGGCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGGGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTggg           <  3:306659/135‑1 (MQ=255)                                                                                                                        cgGTCACCCATGCGGCCGATTAACTGGAAAGTCAGTGCCGGGGCCGTCGGGGGCGGGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTgg            <  8:447518/134‑1 (MQ=255)                                                                                                                        cgGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGTCGACGCCAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGGTGGGTCCATTGCTGGGTTgg            >  3:288606/1‑134 (MQ=255)                                                                                                                        cgGTCACCCATGCGGCCGATGAACTGGAAAGGCAGATCCGGGGCCGTCGGGGCCGCGGCGACGACAGCTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTACATTGCTGGGTTggg           >  6:335935/1‑135 (MQ=255)                                                                                                                            cACCCATGCGGCCTAGGCACTGGAACGCCCTTGCGGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACTGGTTGTTCGGTGACCGTGGCCAGGGTTGTGGGTTGGGTCCATTGCTTGGTTGGGAACa       <  1:24335/135‑1 (MQ=255)                                                                                                                            cACCCATGCGGCCGATGACCTGGAACGGCAGTTCCGGGGCCGGCGGGGCCGCGGCGACGACAGGTGCGGTCATCACAGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTGGGAACa       <  7:113394/135‑1 (MQ=255)                                                                                                                                 ctGCGGCAGTTTAACTGGAAAGGCAGTTCGGGGGCCGTGGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTGGGAACAGGTcg  <  2:168415/134‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    AGACCAGGTTCAGCTCGTTGGCCGATACCCGATCGAGGCGATAGGTGTCTTCGATCACGTCGCCTTGGCGCACGACGTAGAGTTTTTCGCCGCTCTGCAGGAAGACCTGCAGGTCGTCGCGGTCACCCATGCGGCCGATGAACTGGAAAGGCAGTTCCGGGGCCGTCGGGGCCGCGGCGACGACAGGTGCGGTCATCACGGGTTGTTCGGTGACCGTGGCCAGGGTTTTGGGTTGGGTCCATTGCTGGGTTGGGAACAGGTCG  >  ACB122/3937584‑3937846

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