Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	262,746	Δ1 bp	8.2%	coding (263/282 nt)	PP_0211 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	262,746	0	G	.	8.2%	190.5 / 10.3	61	coding (263/282 nt)	PP_0211	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base . (3/2);  total (32/29)
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.38e-01

TCACAATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGC  >  NC_002947/262598‑262894                                                                                                                                                     |                                                                                                                                                      tCACAATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCg                                                                                                                                                       >  6:231968/1‑149 (MQ=255)   acaATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCg                                                                                                                                                       <  3:325506/147‑1 (MQ=255)   acaATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGt                                                                                                                                                      >  1:35236/1‑148 (MQ=255)    caATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGa                                                                                                                                                    >  4:169823/1‑149 (MQ=255)         cagTGCGGCGCGCCTGTCCGGAGTGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATTGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTg                                                                                                                                                <  8:171337/146‑1 (MQ=255)           gTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGt                                                                                                                                              >  3:111897/1‑148 (MQ=255)                    gcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGa                                                                                                                                   <  2:203271/149‑1 (MQ=255)                      cTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAg                                                                                                                                 >  6:80587/1‑149 (MQ=255)                       tGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCGAGGGTATGGCGTACAACTGGTGTTCCGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGc                                                                                                                                >  7:212117/1‑149 (MQ=255)                         tCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGc                                                                                                                                <  5:231968/147‑1 (MQ=255)                              aGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATccc                                                                                                                         <  3:169823/149‑1 (MQ=255)                               gggATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                           <  6:318916/146‑1 (MQ=255)                               gggATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                           >  5:318916/1‑146 (MQ=255)                                       ccTGGTGCGCGACGATGTGCGCGTGGAGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACACTAGGGCTTACTAGCGTGACCTGGACA‑GCAAGTGAAGCTGATCCCTTCAGTAGc                                                                                                                >  2:269869/1‑149 (MQ=255)                                        cTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGcc                                                                                                               <  4:111897/149‑1 (MQ=255)                                         tGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGAc                                                                                                                                                   >  8:164547/1‑113 (MQ=255)                                         tGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGAc                                                                                                                                                   <  7:164547/113‑1 (MQ=255)                                             gcgcGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑TGACCAGGTCacgcaa                                                                                                                                       >  4:245551/1‑115 (MQ=255)                                             gcgcGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑TGACCAGGTCacgcaa                                                                                                                                       <  3:245551/120‑6 (MQ=255)                                                 gacgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGt                                                                                                      >  5:269637/1‑149 (MQ=255)                                                  acgaTGTGCGCGTGGTGCAAATCGAGTCGCGAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTg                                                                                                     <  8:212117/149‑1 (MQ=255)                                                   cgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTg                                                                                                     >  3:89653/1‑148 (MQ=255)                                                        tgCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGaa                                                                                                <  6:33438/148‑1 (MQ=255)                                                            cgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                            >  7:117208/1‑148 (MQ=255)                                                             gtggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGga                                                                                          <  6:125550/149‑1 (MQ=255)                                                             gtggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGga                                                                                          >  4:64632/1‑149 (MQ=255)                                                                  gCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGAcc                                                                                     >  2:123866/1‑149 (MQ=255)                                                                   cAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCg                                                                                    >  8:185293/1‑149 (MQ=255)                                                                   cAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCg                                                                                    <  6:269637/149‑1 (MQ=255)                                                                    aaaTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGc                                                                                   <  2:185105/149‑1 (MQ=255)                                                                      aTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                            >  6:97831/1‑106 (MQ=255)                                                                      aTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                            <  5:97831/106‑1 (MQ=255)                                                                       tCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTa                                                                                <  7:185293/149‑1 (MQ=255)                                                                          aGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGc                                                                              >  6:178061/1‑148 (MQ=255)                                                                              gCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAAc                                                                          <  1:123866/148‑1 (MQ=255)                                                                                         gCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                               <  2:294525/148‑1 (MQ=255)                                                                                                 cTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCt                                                      >  3:110373/1‑149 (MQ=255)                                                                                                            cGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCg                                           <  8:117208/149‑1 (MQ=255)                                                                                                               tGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGaca                                        >  6:347657/1‑149 (MQ=255)                                                                                                                   cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACAcccc                                    >  5:117885/1‑149 (MQ=255)                                                                                                                   cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTAGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGAcccccc                                    >  1:292505/1‑149 (MQ=255)                                                                                                                     tGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGa                                  <  3:64632/149‑1 (MQ=255)                                                                                                                      gAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCa                                                                                                                     >  4:120182/1‑65 (MQ=255)                                                                                                                      gAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCa                                                                                                                     <  3:120182/65‑1 (MQ=255)                                                                                                                       aGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAgt                                <  5:178061/149‑1 (MQ=255)                                                                                                                        gcgcGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAgtg                               >  5:78367/1‑149 (MQ=255)                                                                                                                           cgGGATCTACCCTTGGGCTTTCTTGCGTGAGCTGGTCA‑GGAAGTGAAGCTGATCCCTTCAGTTGCCGGACGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACGTTCTCTTTGGCGACACCCCGAGTGTgg                            <  1:269869/149‑1 (MQ=255)                                                                                                                           cgGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                            <  3:164738/149‑1 (MQ=255)                                                                                                                           cgGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                            <  4:110373/149‑1 (MQ=255)                                                                                                                              gATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTa                         >  8:121012/1‑149 (MQ=255)                                                                                                                                   aCCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCAGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAggg                    >  2:322732/1‑149 (MQ=255)                                                                                                                                      cTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgcg                  >  4:144894/1‑148 (MQ=255)                                                                                                                                      cTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCAGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgcg                  >  2:226701/1‑148 (MQ=255)                                                                                                                                         gggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCaa               <  6:78367/148‑1 (MQ=255)                                                                                                                                          ggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGctc                                                  >  6:158682/1‑112 (MQ=255)                                                                                                                                          ggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAAc              <  3:144894/148‑1 (MQ=255)                                                                                                                                           cttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg            >  6:326456/5‑149 (MQ=255)                                                                                                                                            ttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg           >  5:286978/4‑149 (MQ=255)                                                                                                                                            ttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg           <  5:326456/146‑1 (MQ=255)                                                                                                                                            cTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGACCTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg           <  2:35236/149‑1 (MQ=255)                                                                                                                                            cTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg           <  5:347657/149‑1 (MQ=255)                                                                                                                                                 ttGCGTGACCTGGTCA‑GTAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAAc      >  5:265893/1‑149 (MQ=255)                                                                                                                                                    cGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTg   >  7:287136/1‑149 (MQ=255)                                                                                                                                                     gTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc  >  4:151565/1‑149 (MQ=255)                                                                                                                                                     gTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc  >  4:57153/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCACAATGCCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGC  >  NC_002947/262598‑262894

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