Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,618,418	C→A	7.8%	G161C (GGC→TGC)	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,418	0	C	A	7.8%	199.8 / 4.2	66	G161C (GGC→TGC)	PP_1418	tricarboxylate transport protein TctC
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base A (3/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

AAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACCGCGT  >  NC_002947/1618271‑1618564                                                                                                                                                     |                                                                                                                                                   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:364279/1‑149 (MQ=255)                                                     ttCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTg                                                                                               >  1:201934/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                                                                              <  3:79142/149‑1 (MQ=255)                                                                      gggCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                              <  1:88950/149‑1 (MQ=255)                                                                      gggCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                              <  8:36037/149‑1 (MQ=255)                                                                              tAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTa                                                                      <  1:355789/149‑1 (MQ=255)                                                                                gcgcAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                                       >  3:100567/1‑146 (MQ=255)                                                                                gcgcAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                                       <  4:100567/146‑1 (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                                          <  5:321681/134‑1 (MQ=255)                                                                                                                          aCGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCa                                          >  6:321681/1‑134 (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)                                                                                                                                         aGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCAt            <  2:5688/149‑1 (MQ=255)                                                                                                                                           tCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCAtt           >  8:322580/1‑148 (MQ=255)                                                                                                                                             cTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCa        >  2:50960/1‑149 (MQ=255)                                                                                                                                                 cTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACcgc    <  1:50960/149‑1 (MQ=255)                                                                                                                                                    ccGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACCGCGt  <  3:332620/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTG‑ATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACCGCGT  >  NC_002947/1618271‑1618564

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