Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,618,403	T→G	6.2%	M166L (ATG→CTG)	PP_1418 ←	tricarboxylate transport protein TctC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,618,403	0	T	G	6.2%	197.0 / 3.5	65	M166L (ATG→CTG)	PP_1418	tricarboxylate transport protein TctC
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base G (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

AGTCGGAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGC  >  NC_002947/1618265‑1618551                                                                                                                                            |                                                                                                                                                     aGTCGGAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTgg                                                                                                                                            <  5:364995/149‑1 (MQ=255)     ggAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGc                                                                                                                                        >  5:231717/1‑149 (MQ=255)     ggAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCAAGTAGCGCAGTTCGCGCGGGATG‑ATGCCAGCGGACTTGGCGATCAACGCAGTCTGCTTCCAGTCATGGCTGc                                                                                                                                        >  8:273325/1‑149 (MQ=255)       aaaTGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGcc                                                                                                                                       >  4:312103/1‑148 (MQ=255)          tGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCAGCCAGGACTGGATGCAGAct                                                                                                                                   <  2:138453/149‑2 (MQ=255)             cGGTTGAGCCGACCTGGATGTAGCCGCCCAGCAATGGCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCATGTCGCGCGGGTTGCATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACgg                                                                                                                                  <  1:286555/148‑1 (MQ=255)                     ccGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGa                                                                                                                        <  1:208988/149‑1 (MQ=255)                      cGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGTCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCgg                                                                                                                         <  4:10531/147‑1 (MQ=255)                        aCCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAAcc                                                                                                                     >  1:338282/1‑149 (MQ=255)                          cTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCg                                                                                                                    >  2:170971/1‑148 (MQ=255)                             gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGAt                                                                                                                  <  1:62348/147‑1 (MQ=255)                             gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATc                                                                                                                 <  2:303549/148‑1 (MQ=255)                             gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCa                                                                                                                >  7:54892/1‑149 (MQ=255)                                gtgGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGa                                                                                                             >  4:363720/1‑149 (MQ=255)                                  ggCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGAcc                                                                                                           >  2:88950/1‑149 (MQ=255)                                     cgccCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTg                                                                                                        <  1:42126/149‑1 (MQ=255)                                        ccAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTc                                                                                                     <  4:70135/149‑1 (MQ=255)                                         cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGTCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACTGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                    <  7:44467/149‑1 (MQ=255)                                         cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                    >  4:20321/1‑149 (MQ=255)                                         cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                    >  4:252347/1‑149 (MQ=255)                                             aaTGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTAc                                                                                                                               >  3:55902/1‑118 (MQ=255)                                             aaTGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTAc                                                                                                                               <  4:55902/118‑1 (MQ=255)                                                gCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGtcttt                                                                                             >  4:79142/1‑149 (MQ=255)                                                    tGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTttcttc                                                                                         <  3:312103/149‑1 (MQ=255)                                                    tGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTttcttc                                                                                         <  3:20321/149‑1 (MQ=255)                                                      gCGATTTCGCCGCCGCCCTCCAGGGCAACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAg                                                                                       >  1:5688/1‑149 (MQ=255)                                                           ttCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTg                                                                                  >  1:201934/1‑149 (MQ=255)                                                           ttCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTg                                                                                  >  1:364279/1‑149 (MQ=255)                                                              gccgccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGa                                                                               <  5:41971/149‑1 (MQ=255)                                                               ccgccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGAc                                                                              >  1:14443/1‑149 (MQ=255)                                                                  ccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                           <  3:363720/149‑1 (MQ=255)                                                                     cccTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTGTCTTCAGCGCTGCGACCAGGt                                                                         >  2:355789/1‑148 (MQ=255)                                                                      ccTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGTCCTTGGAGCTAAACGCAGTCTGCATCCAGTCCTGGCTGCAGACGGTAACGCCGGAACCGATAAAGACCTTGCTACGGTCTTTCTTCAGAGCTGCTACCACGtc                                                                        >  1:329144/1‑148 (MQ=255)                                                                       cTCCAGGGCCACGTAGGGAAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGGCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcat                                                                      <  7:273325/149‑1 (MQ=255)                                                                       cTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcat                                                                      <  2:201934/149‑1 (MQ=255)                                                                        tCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                           <  6:68161/143‑1 (MQ=255)                                                                        tCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                           >  5:68161/1‑143 (MQ=255)                                                                          cAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcatc                                                                     <  3:252347/147‑1 (MQ=255)                                                                          cAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCaa                                                                   <  1:170971/149‑1 (MQ=255)                                                                            gggCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                 <  8:36037/149‑1 (MQ=255)                                                                            gggCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                 <  3:79142/149‑1 (MQ=255)                                                                            gggCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                 <  1:88950/149‑1 (MQ=255)                                                                                    tAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTa                                                         <  1:355789/149‑1 (MQ=255)                                                                                      gcgcAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                          <  4:100567/146‑1 (MQ=255)                                                                                      gcgcAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                          >  3:100567/1‑146 (MQ=255)                                                                                       cgcAGGTCGCGGGGGTTG‑ATGCCTGCGGCCTTGGGGATCAAGGATTTGTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGTACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCATGGGTCTTGTAggg                                                      <  8:119400/149‑1 (MQ=255)                                                                                       cgcAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAggg                                                      <  2:14443/149‑1 (MQ=255)                                                                                                       tttaTGCCAGCGGCCTTGGCTATAAACGCAGTCTTAATCCAGTCCTGGCTGCTGACGTTACCGCTGGAACCGATCATGACCTTGCTCGGGTCTTTCTTCAGCGCTTCGACCAGGTCATCAATGGTCTTGTAGGGCGAGTCGCTCTTCAc                                       <  2:329144/146‑1 (MQ=255)                                                                                                        g‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACt                                      <  8:54892/148‑1 (MQ=255)                                                                                                              cAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATg                                >  3:96805/1‑149 (MQ=255)                                                                                                                  ggCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAc                            <  6:334908/149‑1 (MQ=255)                                                                                                                           gATCAACTCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTgg                   >  1:104089/1‑149 (MQ=255)                                                                                                                           gATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTgg                   >  1:142935/1‑149 (MQ=255)                                                                                                                                aCGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCa                             >  6:321681/1‑134 (MQ=255)                                                                                                                                aCGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCa                             <  5:321681/134‑1 (MQ=255)                                                                                                                                 cGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGACGGAACCGATCACGACCTTGCTAGGGTCTTTCTTCAGCGCTGAGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACAGTAGCTGGTGCCGa             >  1:243793/1‑149 (MQ=255)                                                                                                                                  gCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAcc                           >  2:52649/1‑134 (MQ=255)                                                                                                                                  gCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAcc                           <  1:52649/134‑1 (MQ=255)                                                                                                                                    aGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACgg          <  2:104089/149‑1 (MQ=255)                                                                                                                                     gTCGGCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACgg          >  1:63187/1‑148 (MQ=255)                                                                                                                                      tcggCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                 >  7:32136/4‑92 (MQ=38)                                                                                                                                      tcggCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                 <  8:32136/89‑1 (MQ=38)                                                                                                                                      tCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGca        <  2:364279/149‑1 (MQ=255)                                                                                                                                      tCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGAGCCCGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCGTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGca        <  6:231717/149‑1 (MQ=255)                                                                                                                                        tGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGAACCGTAGCTGGTGCCGACGGcagc      >  1:109130/1‑149 (MQ=255)                                                                                                                                            tCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGcagccagc  >  7:93815/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     AGTCGGAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGC  >  NC_002947/1618265‑1618551

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