Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,287,236	C→T	6.4%	A24T (GCC→ACC)	ggt ←	gamma‑glutamyltranspeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,287,236	0	C	T	6.4%	196.9 / 2.0	67	A24T (GCC→ACC)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TTACCGGCCTCTGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAACCGGCGGAAGTT  >  NC_002947/5287091‑5287384                                                                                                                                                  |                                                                                                                                                     ttACCGGCCTCTGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTggcggc                                                                                                                                                   <  2:112804/149‑1 (MQ=255)      ggCCTCTGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTaaa                                                                                                                                              >  7:153884/1‑149 (MQ=255)       gCCTCTGGGTACGTCCCCGCCAGGGTGAAGCCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTATCAGAATATCGGCGGCCACCTGTGCAGCAAATAGATCGGGGGCGGCCAACGCGCCGCCCTCCAGGGTGGCGGCGTAAAc                                                                                                                                             <  1:266908/149‑1 (MQ=255)         ctctGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAAcc                                                                                                                                            >  6:106051/1‑148 (MQ=255)            tGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAg                                                                                                                                         >  8:27465/1‑148 (MQ=255)                tACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGaa                                                                                                                                      >  7:288782/1‑147 (MQ=255)                    tCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAAc                                                                                                                                             <  3:6740/136‑1 (MQ=255)                    tCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAAc                                                                                                                                             >  4:6740/1‑136 (MQ=255)                      aCCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAAc                                                                                                                                             <  3:92790/134‑1 (MQ=255)                      aCCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAAc                                                                                                                                             >  4:92790/1‑134 (MQ=255)                      aCCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCa                                                                                                                               <  2:114549/148‑1 (MQ=255)                        cGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAg                                                                                                                                         >  8:106993/1‑136 (MQ=255)                        cGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAg                                                                                                                                         <  7:106993/136‑1 (MQ=255)                        cGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGa                                                                                                                            >  2:75727/1‑149 (MQ=255)                            aGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCTCCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATcgc                                                                                                                        <  6:362528/149‑1 (MQ=255)                                 gAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTg                                                                                                                   <  8:288782/149‑1 (MQ=255)                                    ggCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCga                                                                                                                <  7:41081/149‑1 (MQ=255)                                    ggCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCga                                                                                                                <  5:106051/149‑1 (MQ=255)                                            cAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAgg                                                                                                        >  5:258646/1‑149 (MQ=255)                                             aaCTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAgg                                                                                                        >  6:110690/1‑148 (MQ=255)                                               cTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTAc                                                                                                     >  7:17455/1‑149 (MQ=255)                                               cTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTAc                                                                                                     <  7:27465/149‑1 (MQ=255)                                                  ccGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTg                                                                                                   <  6:258646/148‑1 (MQ=255)                                                    gcgTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGAGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCga                                                                                                                >  2:305549/1‑133 (MQ=255)                                                    gcgTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGAGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCga                                                                                                                <  1:305549/133‑1 (MQ=255)                                                       tCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACgg                                                                                             >  5:86880/1‑149 (MQ=255)                                                         cACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGa                                                                                           <  1:75727/149‑1 (MQ=255)                                                          aCCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGa                                                                                                                            <  2:115497/115‑1 (MQ=255)                                                          aCCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGa                                                                                           >  6:305489/1‑148 (MQ=255)                                                                 tGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACg                                                                                    >  3:147967/1‑148 (MQ=255)                                                                   ccgccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCa                                                                                  <  5:110690/148‑1 (MQ=255)                                                                    cgccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGTGGCGTAAACCGAg                                                                                                                                         <  4:92876/92‑1 (MQ=255)                                                                    cgccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATa                                                                                      >  6:334637/1‑143 (MQ=255)                                                                    cgccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATa                                                                                      <  5:334637/143‑1 (MQ=255)                                                                    cgccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACAGCGCCGCCCTCCAGGGTGGTGGCGTAAACCGAg                                                                                                                                         >  3:92876/1‑92 (MQ=255)                                                                     gccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTc                                                                                                                                <  2:125388/100‑1 (MQ=255)                                                                     gccCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTc                                                                                                                                >  1:125388/1‑100 (MQ=255)                                                                        cTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGacac                                                                            >  1:105069/1‑149 (MQ=255)                                                                              cAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCt                                                                       <  8:17455/148‑1 (MQ=255)                                                                               aGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTg                                                                     <  2:105069/149‑1 (MQ=255)                                                                                            ccACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTttgt                                                        <  2:83657/149‑1 (MQ=255)                                                                                             cACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTttgtt                                                       <  4:147967/149‑1 (MQ=255)                                                                                               ccTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGa                                                     <  8:153884/149‑1 (MQ=255)                                                                                                 tgtgCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTg                                                                     >  2:74779/1‑131 (MQ=255)                                                                                                 tgtgCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTg                                                                     <  1:74779/131‑1 (MQ=255)                                                                                                   tgCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCCGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGca                                                  <  5:305489/148‑1 (MQ=255)                                                                                                     cACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGcagt                                                <  2:33069/148‑1 (MQ=255)                                                                                                      aCCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGcagtca                                              <  6:86880/149‑1 (MQ=255)                                                                                                        cATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGa                                                     >  5:291766/1‑140 (MQ=255)                                                                                                        cATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGa                                                     <  6:291766/140‑1 (MQ=255)                                                                                                           atTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGgtg                                                                  >  4:34009/1‑124 (MQ=255)                                                                                                           atTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGgtg                                                                  <  3:34009/124‑1 (MQ=255)                                                                                                               aTCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACCCTTCCTTGGTGTTGTTTTTGTTGCGCAGTCAGTAAGGcc                                      >  5:266193/1‑148 (MQ=255)                                                                                                                    gggCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAgcg                                >  2:109734/1‑149 (MQ=255)                                                                                                                        ggcCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTc                            >  8:223875/1‑149 (MQ=255)                                                                                                                          ccACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCAGAAGCGCTTCGa                          >  4:47939/1‑149 (MQ=255)                                                                                                                           cACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGAc                         >  8:55994/1‑149 (MQ=255)                                                                                                                                cgccgcCACCCTGGAGGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGc                                                   >  7:275339/1‑118 (MQ=255)                                                                                                                                cgccgcCACCCTGGAGGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGc                                                   <  8:275339/118‑1 (MQ=255)                                                                                                                                  ccgccCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCaa                  >  1:103511/1‑149 (MQ=255)                                                                                                                                      ccTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCaa              >  4:104485/1‑149 (MQ=255)                                                                                                                                      ccTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCaa              >  1:163645/1‑149 (MQ=255)                                                                                                                                         ccAGGGTGGTGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTc                            <  2:181393/132‑1 (MQ=255)                                                                                                                                         ccAGGGTGGTGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTc                            >  1:181393/1‑132 (MQ=255)                                                                                                                                            gggTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACggg                                                                                            >  2:178836/1‑65 (MQ=255)                                                                                                                                            gggTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACggg                                                                                            <  1:178836/65‑1 (MQ=255)                                                                                                                                             ggTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAACcggcg        >  2:104833/1‑148 (MQ=255)                                                                                                                                              gTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAACcggcgg       >  6:124501/1‑148 (MQ=255)                                                                                                                                                  cggcgTAAACCGAGGAACAGCTCAGGCTCGCGGCTGCGAGAGCCAGGGACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAACCGGCGGAAGtt  <  8:15324/149‑1 (MQ=255)                                                                                                                                                  cggcgTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAAGCGGCGGAAGtt  <  8:10092/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TTACCGGCCTCTGGGTACGTCACCGCCAGGGTGAAGGCGGTCGCAACTGCCGCGTCCACCGCGTTGCCGCCCTTTTTCAGAATATCGGCGGCCACCTGTGCACCATATTGATCGGGGGCGGCCACCGCGCCGCCCTCCAGGGTGGCGGCGTAAACCGAGGAACAGCTCAGGATCGCGGCTGCGAGAGCCAGGTACTGGAACGGGACAATACGCATGACACTTCCTTGGTGTTGTTTTTGTTGAGCAGTCAGTAAGGCCGAAGCGCTTCGACTGTGCAATCAACCGGCGGAAGTT  >  NC_002947/5287091‑5287384

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