Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,108,335	G→C	7.2%	W23C (TGG→TGC)	PP_5358 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,108,335	0	G	C	7.2%	216.8 / 6.2	69	W23C (TGG→TGC)	PP_5358	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base C (2/3);  total (33/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 = 1.49e-01

GCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGT  >  NC_002947/6108188‑6108481                                                                                                                                                     |                                                                                                                                                   gcagcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGc                                                                                                                                                   <  2:210563/149‑1 (MQ=255)   agcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGca                                                                                                                                                  >  1:241136/1‑148 (MQ=255)     cgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACaa                                                                                                                                               <  8:75139/149‑1 (MQ=255)      gcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAAc                                                                                                                                              >  3:221981/1‑149 (MQ=255)         gTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                           <  4:334851/149‑2 (MQ=255)         gTATGAACTTGATTGCTACCTGTACGACAGGCCACGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                           <  3:213002/149‑2 (MQ=255)          tATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACCGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                           <  1:20417/149‑2 (MQ=255)             gAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGt                                                                                                                                       >  7:35501/1‑149 (MQ=255)                   gATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGc                                                                                                                                    <  3:335414/146‑1 (MQ=255)                   gATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGc                                                                                                                                    >  4:335414/1‑146 (MQ=255)                   gATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                  <  7:141603/148‑1 (MQ=255)                           ccTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAgg                                                                                                                           <  8:35501/147‑1 (MQ=255)                           ccTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGcc                                                                                                                         >  5:59527/1‑149 (MQ=255)                             tGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCt                                                                                                                       >  2:181911/1‑149 (MQ=255)                             agTACGACAGGCAAGGAACGCAGTTAAGCTGGATCTAACTGTTTGGAGCAAATGCAAATGTAGACCCAGGAGAACCTCTAAGTCCAGATCCTGGAATCCA‑CGGAAGGCGCCTGGTTGTGGCACAACCAAGGGTTGCGTGAAAAggacgt                                                                                                                       >  4:256150/2‑145 (MQ=255)                                  gACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                      <  2:174028/145‑1 (MQ=255)                                  gACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                      >  1:174028/1‑145 (MQ=255)                                    cAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACgg                                                                                                                <  1:160454/149‑1 (MQ=255)                                    cAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACgg                                                                                                                <  6:54394/149‑1 (MQ=255)                                           ggAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAAc                                                                                                         <  1:181911/149‑1 (MQ=255)                                                 cAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCAtt                                                                                                   >  5:142689/1‑149 (MQ=255)                                                  aGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGc                                                                                                                                                   >  7:236615/1‑100 (MQ=255)                                                  aGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGc                                                                                                                                                   <  8:236615/100‑1 (MQ=255)                                                   gTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                 <  4:315899/149‑1 (MQ=255)                                                     tAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTg                                                                                                  >  7:823/1‑146 (MQ=255)                                                     tAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTg                                                                                                  <  8:823/146‑1 (MQ=255)                                                         cTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCaaa                                                                                           <  8:105359/149‑1 (MQ=255)                                                            gATCTACCTGTTTGGAGGATATGCAAATGTTGACCCTGGAGAACCTCTACGTCCTGATTGTGGTAGCCA‑CGGCAGTCGCCTGGTTGTTGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCaaa                                                                                           <  3:256150/146‑1 (MQ=255)                                                                  ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTg                                                                                  >  1:241398/1‑149 (MQ=255)                                                                   cTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTg                                                                                  >  3:34090/1‑148 (MQ=255)                                                                     gTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCaaa                                                                                                          >  4:118341/1‑122 (MQ=255)                                                                     gTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCaaa                                                                                                          <  3:118341/122‑1 (MQ=255)                                                                     gTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGACCTGATGCTGGTAGCCA‑CCGCAGGCGCCTGGTTGTGGCACAACCATCGGTTGCGTGAAAACGCCCTGGAACGCGTCAAACAGCATTGCGCCAAACTCGACCTGGAg                                                                               >  2:142863/1‑149 (MQ=255)                                                                          gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCt                                                                          >  7:271226/1‑149 (MQ=255)                                                                          gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCt                                                                          <  2:215981/149‑1 (MQ=255)                                                                             cATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCg                                                                                       <  1:148949/133‑1 (MQ=255)                                                                             cATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCg                                                                                       >  2:148949/1‑133 (MQ=255)                                                                                   cAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGc                                                                 >  3:118137/1‑149 (MQ=255)                                                                                    aaaTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGcc                                                                >  2:127506/1‑149 (MQ=255)                                                                                         ttGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAg                                                            <  2:241398/148‑1 (MQ=255)                                                                                              ccTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTc                                                       >  4:271441/1‑148 (MQ=255)                                                                                               cTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                     <  2:241136/149‑1 (MQ=255)                                                                                                  gagaACCTCTTCGTCCTGATGCTGGTAGTCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGt                                                                              >  3:143018/1‑121 (MQ=255)                                                                                                  gagaACCTCTTCGTCCTGATGCTGGTAGTCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGt                                                                              <  4:143018/121‑1 (MQ=255)                                                                                                       ccTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGTAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCg                                             <  3:94784/149‑1 (MQ=255)                                                                                                       ccTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                     <  8:273922/141‑1 (MQ=255)                                                                                                       ccTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                     >  7:273922/1‑141 (MQ=255)                                                                                                       ccTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCg                                             >  5:150512/1‑149 (MQ=255)                                                                                                        ctctTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGc                                            >  7:201055/1‑149 (MQ=255)                                                                                                          ctTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                          <  4:221981/149‑1 (MQ=255)                                                                                                           ttCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACAATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCTCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTAACGCCtt                                         >  6:175832/1‑149 (MQ=255)                                                                                                             cGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                          <  8:201055/146‑1 (MQ=255)                                                                                                                      gCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATg                              >  8:84475/1‑149 (MQ=255)                                                                                                                               cA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                     <  6:169508/149‑1 (MQ=255)                                                                                                                               cA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                     <  6:142689/149‑1 (MQ=255)                                                                                                                               cA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                     >  5:169508/1‑149 (MQ=255)                                                                                                                               cA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                     >  2:187054/1‑149 (MQ=255)                                                                                                                                  cGGCAGGCGCCTGGTTGTGGCACTACCATGGGTTGCGTGAATAGGCCCTGTAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                   <  1:142863/149‑1 (MQ=255)                                                                                                                                  cGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                   <  1:187054/149‑1 (MQ=255)                                                                                                                                          gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt           <  6:281121/149‑1 (MQ=255)                                                                                                                                          gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt           >  5:90390/1‑149 (MQ=255)                                                                                                                                          gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt           <  8:271226/149‑1 (MQ=255)                                                                                                                                           catGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg          >  4:310875/3‑149 (MQ=255)                                                                                                                                            cTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg         >  4:301395/1‑149 (MQ=255)                                                                                                                                            atGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg         >  5:341642/2‑149 (MQ=255)                                                                                                                                               gTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAc      <  6:150512/149‑1 (MQ=255)                                                                                                                                                  gtgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc    <  7:84475/148‑1 (MQ=255)                                                                                                                                                  gtgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc    <  1:137508/148‑1 (MQ=255)                                                                                                                                                   tgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGt  >  2:213961/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCA‑CGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGT  >  NC_002947/6108188‑6108481

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