Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,048,228	G→A	5.6%	P104P (CCG→CCA)	aroC →	Chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,048,228	0	G	A	5.6%	213.9 / 3.1	71	P104P (CCG→CCA)	aroC	Chorismate synthase
Reads supporting (aligned to +/- strand):  ref base G (30/37);  new base A (2/2);  total (32/39)
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-01

CAGGAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCACCCAGGG  >  ACB122/2048094‑2048362                                                                                                                                       |                                                                                                                                       cAGGATGCCGAGGAGGGGGAAATCCTTTCGGGCGTGTGCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCg                                                                                                                                        <  8:170390/135‑1 (MQ=255)    gAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGcc                                                                                                                                     >  2:368946/1‑135 (MQ=255)    gAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGcc                                                                                                                                     <  8:84802/135‑1 (MQ=255)    gAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGcc                                                                                                                                     >  8:406977/1‑135 (MQ=255)    gAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGcc                                                                                                                                     >  2:225320/1‑135 (MQ=255)     aaGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGccc                                                                                                                                    <  7:218708/135‑1 (MQ=255)        ccGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACg                                                                                                                                 <  2:337736/135‑1 (MQ=255)           acgaGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCg                                                                                                                              <  8:96606/135‑1 (MQ=255)           acgaGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCg                                                                                                                              >  2:379534/1‑135 (MQ=255)           acgaGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCg                                                                                                                              >  8:303104/1‑135 (MQ=255)                   gAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACAcc                                                                                                                      <  5:106383/135‑1 (MQ=255)                      aTCCTTTCCGGCGTGTTCGAAGGGCGCACCACCGGCTGGTCGTTCGGCCGGCTGATCCGCAATTCCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGTCCTGTTCCGCCCGGCCCTGGCCGACTACACCTac                                                                                                                   <  1:52531/135‑1 (MQ=255)                      aTCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTac                                                                                                                   <  6:378575/135‑1 (MQ=255)                      aTCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGATGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTac                                                                                                                   <  8:354321/135‑1 (MQ=255)                          tttCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTaccac                                                                                                                <  2:105279/134‑1 (MQ=255)                              cGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACaa                                                                                                           >  3:305331/1‑135 (MQ=255)                               ggCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACaa                                                                                                           <  4:390214/134‑1 (MQ=255)                               ggCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCCACTACACCTACCACCACaa                                                                                                           >  5:329935/1‑134 (MQ=255)                                    gtTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTaca                                                                                                                        >  3:427901/1‑116 (MQ=255)                                    gtTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTaca                                                                                                                        <  4:427901/116‑1 (MQ=255)                                      tCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTa                                                                                                   >  4:236071/1‑135 (MQ=255)                                        gAAGCCCGCACCACCGGCGGCTCGATCGGCCTGCTGATCCGCAATACCGATCATAAGGCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATc                                                                                                 <  3:253964/135‑1 (MQ=255)                                           ggCCGCACCCCCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCcgc                                                                                              <  6:359288/135‑1 (MQ=255)                                           gcccGCACCACCGGCTGCTCGATCGGCCTGCTGAGCCGCAAAACCGTTCATAACTCCAAGGACAACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCcgc                                                                                              <  5:265012/133‑1 (MQ=255)                                             ccGCAGCACCGGCTGCTCGATCGGCCTGGTGATCCGAATTACCGATCTGATGTCCATGCACTACTCGGCCATCATGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGa                                                                                            <  6:270640/135‑1 (MQ=255)                                                  cccccGGCTGCGCGAGCGGCCGGCGGATCCGCAATCCCCATCAGACGGCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTAcc                                                                                       <  5:281896/132‑1 (MQ=255)                                                    accGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACcgc                                                                                     <  6:20854/135‑1 (MQ=255)                                                    accGGCTGCTCGATCGGCCTGCTCATCCGCAATACCCATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCAGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACcgc                                                                                     <  2:359368/135‑1 (MQ=255)                                                    accGGCTGCTCGAACGGCCTGCTGATCCGCAAAACCCAACAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGCTCCGCCCAGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACcgc                                                                                     <  2:353260/135‑1 (MQ=255)                                                          tgctCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCggtggt                                                                               <  4:135744/135‑1 (MQ=255)                                                          tgctCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCggtggt                                                                               <  7:406977/135‑1 (MQ=255)                                                                aTCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACGACCCCTACCACCACACGTACGGTATCCGCGACTACCGCGGTGGTGGCCGc                                                                         >  5:445155/1‑135 (MQ=255)                                                                   ggCCTGCTGAGCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAACTACGGTATCCCCGACTACCGCCGAGGTGGCcgccgc                                                                      >  1:86521/1‑135 (MQ=255)                                                                          tgATCCGCAATCCCCTTCAGTAGTCCAAGGACTACTCGGCCATCAAGGAGCTGTTCCGCCCGGCCCACGCCGACTACACCTCCCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGc                                                                 <  1:314233/133‑1 (MQ=255)                                                                          tgATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGccc                                                               <  1:245707/135‑1 (MQ=255)                                                                            aTCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCGCGGCCCGt                                                             >  1:431495/1‑135 (MQ=255)                                                                                 cAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGCCCCACGCCGCCTACACCTACCACCACAACTACGGTATCCGCGACTACCGCGGTGGTGGCCGCCGCACGGCCCGGGAAAc                                                        >  4:41995/1‑135 (MQ=255)                                                                                   aTACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCg                                                      <  4:173394/135‑1 (MQ=255)                                                                                      ccGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGCTGGTGGCCGCACCTCGGCCCGTGAACCCGCCa                                                   >  4:281989/1‑135 (MQ=255)                                                                                         aTCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATgc                                                >  3:107064/1‑135 (MQ=255)                                                                                            aggaGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGcgcg                                             >  6:230058/4‑135 (MQ=255)                                                                                            aGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACCCCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGcgcg                                             >  2:49703/1‑135 (MQ=255)                                                                                            aGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCCTGcgcg                                             >  7:222424/1‑135 (MQ=255)                                                                                            aGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGcgcg                                             >  8:242285/1‑135 (MQ=255)                                                                                            aGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGACCGCAGCTCGGc                                                                 >  6:374309/1‑115 (MQ=255)                                                                                            aGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACAACAAGTACGGTATCCGCGACTACCGCGGTGGTGACCGCAGCTCGGc                                                                 <  5:374309/115‑1 (MQ=255)                                                                                                        actactCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGgcg                                 <  5:192403/135‑1 (MQ=255)                                                                                                        actactCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGgcg                                 <  7:19159/135‑1 (MQ=255)                                                                                                              cGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGCCTACACCTACCCCCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCg                           >  8:229736/1‑135 (MQ=255)                                                                                                              cGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCg                           >  8:57596/1‑135 (MQ=255)                                                                                                              cGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCg                           >  1:110694/1‑135 (MQ=255)                                                                                                              cGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGACATGCGCGTGGCCGCTGGCGCTATCg                           >  4:345316/1‑135 (MQ=255)                                                                                                              cGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGGGAAACCGCCATGCGCGTGGCCGCGGGCGCTCTCg                           >  2:247722/1‑135 (MQ=255)                                                                                                               ggCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCg                           >  6:226476/1‑134 (MQ=255)                                                                                                                    ttaaGGACGTTTTCCGTCCAGGCCATGCCGATTACACCTACGAACAAAAATACGGTCTGCGCGATTATCGCGGCGGTGGACGTTCTTCCGCCCGCGAAACCGCCATGCGCGTGGCGGCAGGAGCTATTGCCaaa                      >  6:166557/3‑133 (MQ=255)                                                                                                                    ttaaGGACGTTTTCCGTCCAGGCCATGCCGATTACACCTACGAACAAAAATACGGTCTGCGCGATTATCGCGGCGGTGGACGTTCTTCCGCCCGCGAAACCGCCATGCGCGTGGCGGCAGGAGCTATTGCCaaa                      >  6:158772/3‑133 (MQ=255)                                                                                                                    tCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCaa                       >  6:157258/1‑133 (MQ=255)                                                                                                                         gACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTc                <  5:344374/135‑1 (MQ=255)                                                                                                                          aCCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCTCAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAATCCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTcc               <  7:257734/135‑1 (MQ=255)                                                                                                                          aCCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTcc               <  4:21114/135‑1 (MQ=255)                                                                                                                               ttCCGCCCGGCCCACGCCGACTACACCTACCACCACCAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTTGCGCTATCGCCAAGAAGTTCCTGGcc          <  6:186293/135‑1 (MQ=255)                                                                                                                               ttCCGCCCGGCCCACGCCGACTACACCTACCACCACAACTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  7:303104/135‑1 (MQ=255)                                                                                                                                tCCGCCCGGCCCCCGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGct          <  1:197050/134‑2 (MQ=255)                                                                                                                                tCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  7:34346/134‑1 (MQ=255)                                                                                                                                tCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  7:5078/134‑1 (MQ=255)                                                                                                                                tCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  7:242285/134‑1 (MQ=255)                                                                                                                                tCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  7:12041/134‑1 (MQ=255)                                                                                                                                gccGCCCGGCCCACGCCGACTACACCTCCCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGcc          <  6:56537/133‑1 (MQ=255)                                                                                                                                 ccGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCAc        >  1:362102/1‑135 (MQ=255)                                                                                                                                  cGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGCCTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCAGGCGCTATCGCCAAGAAGTTCCTGGCCAcc       >  4:61524/1‑135 (MQ=255)                                                                                                                                  cGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCAcc       >  2:252226/1‑135 (MQ=255)                                                                                                                                  cGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCAc        <  7:30084/134‑1 (MQ=255)                                                                                                                                       ggCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCCGCTCGGCCCGGGAAACCCCCATGCGCGTGGCCGCAGGCGCTTTAGCCAAGAAGTTACAGGccccccccg   >  5:398710/1‑131 (MQ=255)                                                                                                                                       ggCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCACCCAggg  >  7:370326/1‑135 (MQ=255)                                                                                                                                       ggCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCACCCAgg   >  6:163530/1‑134 (MQ=255)                                                                                                                                       ggCCCACGCCGACTACACCTACCACCACAACTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCACCCAggg  >  6:243450/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                       CAGGAAGCCGACGAGGTGGAAATCCTTTCCGGCGTGTTCGAAGGCCGCACCACCGGCTGCTCGATCGGCCTGCTGATCCGCAATACCGATCAGAAGTCCAAGGACTACTCGGCCATCAAGGACCTGTTCCGCCCGGCCCACGCCGACTACACCTACCACCACAAGTACGGTATCCGCGACTACCGCGGTGGTGGCCGCAGCTCGGCCCGTGAAACCGCCATGCGCGTGGCCGCTGGCGCTATCGCCAAGAAGTTCCTGGCCACCCAGGG  >  ACB122/2048094‑2048362

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