Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,387,265	G→T	6.1%	R40L (CGC→CTC)	PP_1207 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,387,265	0	G	T	6.1%	210.1 / 3.9	66	R40L (CGC→CTC)	PP_1207	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base T (2/2);  total (30/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.55e-01

GGATGTATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTG  >  NC_002947/1387120‑1387413                                                                                                                                                  |                                                                                                                                                     ggATGTATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAgcgcgc                                                                                                                                                   >  8:277662/1‑149 (MQ=255)  gATGTATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCt                                                                                                                                                  <  5:304315/149‑1 (MQ=255)  gATGTATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCt                                                                                                                                                  >  7:205970/1‑149 (MQ=255)      tATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTa                                                                                                                                              <  7:277662/149‑1 (MQ=255)          gCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCa                                                                                                                                           <  2:24054/148‑1 (MQ=255)           ccTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGGTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCTATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGt                                                                                                                                         <  5:51422/149‑1 (MQ=255)           ccTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTCTGCAAGTTGCTGATGAGCAAGGATGCCAATTACCCATGGTTCATTATGGTGCCCAAGCGCGCTGGTATCAGt                                                                                                                                         >  7:62162/1‑149 (MQ=255)            cTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAgcg                                                                                                                                        >  3:15839/1‑147 (MQ=255)                 ggAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTg                                                                                                                                   >  4:348433/1‑149 (MQ=255)                   agaTATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgg                                                                                                                                 <  4:15839/149‑4 (MQ=255)                         gCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCt                                                                                                                           >  1:197829/1‑149 (MQ=255)                          ccGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCt                                                                                                                           <  8:205970/148‑1 (MQ=255)                                  gTCCTCGATTCGCGGTTGCTGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   <  3:23952/148‑1 (MQ=255)                                   tCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGACCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCa                                                                                                                 >  7:338169/1‑149 (MQ=255)                                     cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGa                                                                                                                              >  6:143751/1‑134 (MQ=255)                                     cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGa                                                                                                                              <  5:143751/134‑1 (MQ=255)                                     cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   <  3:196384/145‑1 (MQ=255)                                     cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   >  4:196384/1‑145 (MQ=255)                                     cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCaga                                                                                                               <  1:227918/149‑1 (MQ=255)                                        gATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGc                                                                                                             >  2:114724/1‑148 (MQ=255)                                                tttGCAGCAGGATTCGCGGGTGTTGGGAGAGTTTCCGGTGTGCCAGTTGCTGGTGAGCAAGGATGCCAATTACCCATGCTTCATTCTGGTGCCGAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGt                                                                                                    <  5:109604/149‑1 (MQ=255)                                                   gcagcaGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTtg                                                                                                  <  3:348433/148‑1 (MQ=255)                                                    cagcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgggc                                                                                                                               <  6:359261/118‑6 (MQ=255)                                                    cagcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgggc                                                                                                                               >  5:359261/1‑113 (MQ=255)                                                      gcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTgg                                                                                               >  7:134017/1‑148 (MQ=255)                                                        agGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAg                                                                                            >  8:54694/1‑149 (MQ=255)                                                          gATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCag                                                                                                                <  8:275952/127‑1 (MQ=255)                                                          gATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCag                                                                                                                >  7:275952/1‑127 (MQ=255)                                                              cGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAac                                                                                       <  4:288226/148‑1 (MQ=255)                                                               gCTGGTGTTGGTAGATTTTCCGCTGTGGCAGTTGCTGCTGAGCAAGGATGCCTATTTCCCTTGGTTCTTTCTTGTGGCCAAGCGCTGTTGTATCAGTGAGCTGTTCTATCTGGTTTCGGGAGAGCAACTTCAGTTGTGGCTGGATacca                                                                                     <  8:62162/149‑1 (MQ=255)                                                               gCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAacca                                                                                     <  8:338169/149‑1 (MQ=255)                                                                                   cGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGtcg                                                                 <  8:188455/149‑3 (MQ=255)                                                                                   cGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTg                                                                 <  4:79017/149‑1 (MQ=255)                                                                                          ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGaaa                                                              >  2:249234/1‑145 (MQ=255)                                                                                          ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGaaa                                                              <  1:249234/145‑1 (MQ=255)                                                                                          ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           <  8:134017/148‑1 (MQ=255)                                                                                                 ctgctgAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATgg                                                   <  6:103121/149‑1 (MQ=255)                                                                                                      gAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGGGCTGGTATCAGTGAGCTGTTCGAGCTGGAt                                                                                                                       <  1:31685/76‑1 (MQ=255)                                                                                                      gAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGAt                                                                                                                       >  2:31685/1‑76 (MQ=255)                                                                                                        gCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGAc                                            <  2:197829/149‑1 (MQ=255)                                                                                                            ggATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGAc                                            >  4:60891/1‑145 (MQ=255)                                                                                                                 ccAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAAc                                   <  3:118757/149‑1 (MQ=255)                                                                                                                   aaTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGt                                 >  8:12369/1‑149 (MQ=255)                                                                                                                   aaTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGt                                 >  4:39233/1‑149 (MQ=255)                                                                                                                    aTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTGGAAGGCCAGCTATGGGGGCGACAAGATGAACGTg                                >  8:322083/1‑149 (MQ=255)                                                                                                                     ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTAt                                                     <  6:300203/127‑1 (MQ=255)                                                                                                                     ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTAt                                                     >  5:300203/1‑127 (MQ=255)                                                                                                                     ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTgg                               <  7:12369/149‑1 (MQ=255)                                                                                                                         ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGTAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           >  4:319169/1‑117 (MQ=255)                                                                                                                         ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGTAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           <  3:319169/117‑1 (MQ=255)                                                                                                                         ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCACCTATGGGGCCGACAAGATGAACGTGGCTa                            >  7:239580/1‑148 (MQ=255)                                                                                                                           aTGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGACCGTGGCGACCt                         >  8:235326/1‑149 (MQ=255)                                                                                                                             ggTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGa                                        >  7:66358/1‑132 (MQ=255)                                                                                                                             ggTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGa                                        <  8:66358/132‑1 (MQ=255)                                                                                                                                tCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTggg                     <  7:658/148‑1 (MQ=255)                                                                                                                                tCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTggg                     >  8:658/1‑148 (MQ=255)                                                                                                                                   ttCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACaa                                          <  6:39952/124‑1 (MQ=255)                                                                                                                                   ttCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACaa                                          >  5:39952/1‑124 (MQ=255)                                                                                                                                      tGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTg              >  3:149278/1‑149 (MQ=255)                                                                                                                                        gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCATAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTg              <  4:324502/147‑1 (MQ=255)                                                                                                                                        gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCTACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            >  8:70360/1‑149 (MQ=255)                                                                                                                                        gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGACCGTGGCTACCTTGGGCAACGTg              <  7:54694/147‑1 (MQ=255)                                                                                                                                        gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            <  4:149278/149‑1 (MQ=255)                                                                                                                                         tGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTc           <  3:60891/149‑1 (MQ=255)                                                                                                                                          gCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            >  3:155888/1‑147 (MQ=255)                                                                                                                                            ccAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGTCAGGAAACGTCCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGc        <  1:114724/149‑1 (MQ=255)                                                                                                                                            ccAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGATGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGc        >  6:105353/1‑149 (MQ=255)                                                                                                                                                  gcgcTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTg  <  4:155888/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGATGTATTGCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTG  >  NC_002947/1387120‑1387413

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