Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,865,753	C→T	6.8%	E342K (GAA→AAA)	glllA ←	gallate dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,865,753	0	C	T	6.8%	173.2 / 2.4	59	E342K (GAA→AAA)	glllA	gallate dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/2);  total (31/28)
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.54e-01

TAAATGTCCAGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGATGCTGCATGTG  >  NC_002947/2865605‑2865900                                                                                                                                                      |                                                                                                                                                     tAAATGTCCAGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTc                                                                                                                                                      >  4:240070/1‑149 (MQ=255)    aTGTCCAGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCggg                                                                                                                                                   <  8:143742/149‑1 (MQ=255)       tCCAGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTc                                                                                                                                                <  3:261460/149‑1 (MQ=255)          aGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCt                                                                                                                                             >  6:111291/1‑149 (MQ=255)                gACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGa                                                                                                                                        >  4:154008/1‑148 (MQ=255)                 acacacCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGaaa                                                                                                                                      >  8:240140/1‑149 (MQ=255)                 acacacCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGaaa                                                                                                                                      >  3:207464/1‑149 (MQ=255)                     acCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCt                                                                                                                                  >  2:66352/1‑149 (MQ=255)                       ccAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTgg                                                                                                                                >  3:39549/1‑149 (MQ=255)                        cAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGgc                                                                                                                               <  8:3963/149‑1 (MQ=255)                        cAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGgc                                                                                                                               <  6:273926/149‑1 (MQ=255)                        cAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGgc                                                                                                                               <  3:240070/149‑1 (MQ=255)                                 gcCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCg                                                                                                                      >  5:98803/1‑149 (MQ=255)                                  ccTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCgg                                                                                                                     <  4:207464/149‑1 (MQ=255)                                           tcCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCAt                                                                                                            >  5:70927/1‑149 (MQ=255)                                             cAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATg                                                                                                           <  6:112853/148‑1 (MQ=255)                                             cAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATg                                                                                                           <  1:66352/148‑1 (MQ=255)                                             cAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGc                                                                                                          <  7:116856/149‑1 (MQ=255)                                              aGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGc                                                                                                                    >  3:269443/1‑138 (MQ=255)                                                   aaaGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGc                                                                                                     >  3:334397/1‑148 (MQ=255)                                                     agaagaTCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGc                                                                                                     >  8:14079/1‑146 (MQ=255)                                                        agaTCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAAc                                                                                               <  4:39549/149‑1 (MQ=255)                                                         gaTCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACt                                                                                              >  8:190651/1‑149 (MQ=255)                                                             aCCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGt                                                                                           >  2:258852/1‑148 (MQ=255)                                                             aCCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGt                                                                                           <  7:190651/148‑1 (MQ=255)                                                                       ggATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAg                                                                                <  7:240140/149‑1 (MQ=255)                                                                             gCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTg                                                                         >  8:72607/1‑149 (MQ=255)                                                                                       gTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCt                                                               >  5:255835/1‑149 (MQ=255)                                                                                            ggcggcgCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAgg                                                          <  6:98803/149‑1 (MQ=255)                                                                                            ggcggcgCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAgg                                                          <  5:111291/149‑1 (MQ=255)                                                                                            ggcggcgCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAgg                                                          <  4:334397/149‑1 (MQ=255)                                                                                                  gcagcaGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAg                                                    >  4:25057/1‑149 (MQ=255)                                                                                                    agcagGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAggg                                                  <  3:154008/149‑1 (MQ=255)                                                                                                             cTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGcc                                         >  1:251925/1‑149 (MQ=255)                                                                                                                tCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGCTTTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGt                                                >  1:61826/1‑140 (MQ=255)                                                                                                                tCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGt                                                <  2:61826/139‑1 (MQ=255)                                                                                                                   tctTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCTGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAg                                   >  7:272359/1‑149 (MQ=255)                                                                                                                    ctTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgc                                  >  2:194733/1‑149 (MQ=255)                                                                                                                       cGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTTGGGGTAGGTGCCTTCAAgctg                                <  3:8936/148‑1 (MQ=255)                                                                                                                       cGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgctg                                <  6:255835/148‑1 (MQ=255)                                                                                                                        gCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgctgc                               <  2:251925/148‑1 (MQ=255)                                                                                                                          tCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTc                             >  4:76603/1‑148 (MQ=255)                                                                                                                            aGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAGCCCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCa                            <  1:258852/148‑1 (MQ=255)                                                                                                                               cccGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGcccc                       >  2:24394/1‑149 (MQ=255)                                                                                                                               cccGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGcccc                       <  3:25057/149‑1 (MQ=255)                                                                                                                                ccGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCt                                        >  1:332388/1‑131 (MQ=255)                                                                                                                                ccGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCt                                        <  2:332388/131‑1 (MQ=255)                                                                                                                                ccGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCt                      >  3:70596/1‑149 (MQ=255)                                                                                                                                   ccTCGCGGTACAGCGC‑TTTGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgc                                  <  4:184953/134‑1 (MQ=255)                                                                                                                                   ccTCGCGGTACAGCGC‑TTTGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgc                                  >  3:184953/1‑134 (MQ=255)                                                                                                                                   ccTCGCGGTACAGCGC‑TTTGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgc                                  >  3:189290/1‑134 (MQ=255)                                                                                                                                   ccTCGCGGTACAGCGC‑TTTGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAgc                                  <  4:189290/134‑1 (MQ=255)                                                                                                                                    cTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGc                    >  7:296223/1‑147 (MQ=255)                                                                                                                                        cgGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGc                <  3:16493/147‑1 (MQ=255)                                                                                                                                        cgGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGc                >  4:16493/1‑147 (MQ=255)                                                                                                                                         ggTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTCTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGa             >  8:113032/1‑149 (MQ=255)                                                                                                                                              aGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGAtgctg        <  3:272153/149‑1 (MQ=255)                                                                                                                                               gcgc‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGAtgctgc       <  7:14079/149‑1 (MQ=255)                                                                                                                                                  c‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGATGCTGCAtg    >  3:290541/1‑149 (MQ=255)                                                                                                                                                    ttCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTTGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGATGCTGCAtgt   >  3:54559/1‑149 (MQ=255)                                                                                                                                                      cGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCTGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGATGCTGCAtgtg  <  8:272359/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TAAATGTCCAGGTTGGACACACCCAGTACGGCGCCTAGCTTCTCCAGTACAAAGAAGATCACCCCATAGTGGATCAGCCCACGCCAGTCGCGGCGGCGCAGCAGGTCGCTTTCCTCTTCGCTCAGGCCCGCCTCGCGGTACAGCGC‑TTCGGGTTCGCTCAGGAAACGCTGGCGCCACTGCGGCTCAATCATGCGATGCAAGAACTTGTTCAGGCGGTAG‑CCCTTGGCACTGCGCTCCAGGGTGTAGGGGTAGGTGCCTTCAAGCTGCTCAGCCCCTGCCAGCTGATGCTGCATGTG  >  NC_002947/2865605‑2865900

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