Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,038,518:1	+A	9.8%	coding (64/255 nt)	PP_4444 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,038,518	1	.	A	9.8%	222.6 / 20.9	61	coding (64/255 nt)	PP_4444	transposase
Reads supporting (aligned to +/- strand):  ref base . (28/27);  new base A (4/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

TGCTGACGCGGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACAT  >  NC_002947/5038371‑5038666                                                                                                                                                      |                                                                                                                                                     tgCTGACGCGGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑a                                                                                                                                                     <  3:346955/149‑1 (MQ=255)     gACGCGGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑Agcgc                                                                                                                                                 >  2:141371/1‑149 (MQ=255)          ggTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTggg                                                                                                                                            >  3:291773/1‑149 (MQ=255)           gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTggg                                                                                                                                            >  2:284978/1‑148 (MQ=255)           gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                           <  4:291773/149‑1 (MQ=255)           gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                           <  4:350375/149‑1 (MQ=255)           gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                           <  7:281293/149‑1 (MQ=255)           gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                           <  1:141371/149‑1 (MQ=255)               ccGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGt                                                                                                                                       >  4:249497/1‑149 (MQ=255)                cGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTc                                                                                                                                      >  3:240880/1‑149 (MQ=255)                 gACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCCAAGCGC‑TGGGCaacatc                                                                                                                                     <  6:45228/149‑6 (MQ=18)                      aGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGt                                                                                                                                   >  8:187027/1‑146 (MQ=255)                      aGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGt                                                                                                                                   <  7:187027/146‑1 (MQ=255)                      aGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGca                                                                                                                                <  1:284978/149‑1 (MQ=255)                           gCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTc                                                                                                                           <  3:249497/149‑1 (MQ=255)                               gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGGCCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                       <  6:72760/149‑1 (MQ=255)                               gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                       >  8:208647/1‑149 (MQ=255)                               gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                       <  8:276163/149‑1 (MQ=255)                               gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                       >  7:276163/1‑149 (MQ=255)                                tGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTAc                                                                                                                      >  4:268893/1‑149 (MQ=255)                                tGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTAc                                                                                                                      <  2:191788/149‑1 (MQ=255)                                     tCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTg                                                                                                                 <  3:268893/149‑1 (MQ=255)                                         gCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGAt                                                                                                              >  4:32833/1‑148 (MQ=255)                                                           gTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGccc                                                                                           >  2:222164/1‑149 (MQ=255)                                                                tATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGa                                                                                      <  7:208647/149‑1 (MQ=255)                                                                tATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGa                                                                                      <  1:222164/149‑1 (MQ=255)                                                                    ggCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCGAAGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTGTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGag                                                                                     >  1:280922/1‑147 (MQ=255)                                                                     gCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGagcagc                                                                                 <  2:122728/149‑1 (MQ=255)                                                                      ccTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCg                                                                                >  2:125899/1‑149 (MQ=255)                                                                            ttACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑Agc                                                                                                                                                   >  3:72222/1‑76 (MQ=255)                                                                            ttACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑Agc                                                                                                                                                   <  4:72222/76‑1 (MQ=255)                                                                            ttACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGc                                                                                                 <  5:135378/126‑1 (MQ=255)                                                                            ttACGCAAGGTGCGTAGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGc                                                                                                 >  6:135378/1‑126 (MQ=255)                                                                              aCGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTcaac                                                                           >  5:100248/1‑146 (MQ=255)                                                                                    ggTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCATAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACgatg                                                                  <  1:105405/149‑1 (MQ=255)                                                                                             ggCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAgg                                                         >  8:74154/1‑149 (MQ=255)                                                                                                    gCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACt                                                  >  6:277400/1‑149 (MQ=255)                                                                                                      ttACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTgc                                                >  4:221055/1‑149 (MQ=255)                                                                                                        aCAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTgc                                                <  3:32833/147‑1 (MQ=255)                                                                                                          aagcaagTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTgcgc                                              >  6:165273/1‑147 (MQ=255)                                                                                                                   aCCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTgcgcgc                                   <  7:74154/149‑1 (MQ=255)                                                                                                                     ag‑agaAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATAAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCt                                  >  2:5750/2‑148 (MQ=255)                                                                                                                        agaAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGaa                               >  1:93018/1‑149 (MQ=255)                                                                                                                         gaAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAAc                              <  3:162293/149‑1 (MQ=255)                                                                                                                               cAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCTAACTGCGCAAGCTGCGCGCTGAACTCAAg                         <  1:5750/148‑1 (MQ=255)                                                                                                                                 aaCCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGTTCAAGCTATACAGCAAGCCCCAGGAGCAGCGGCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCgg                      <  1:286578/149‑1 (MQ=255)                                                                                                                                  aCCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCggg                     <  3:221055/149‑1 (MQ=255)                                                                                                                                     tGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGa                  >  3:156562/1‑149 (MQ=255)                                                                                                                                       ttGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTgg                                                                                                                <  2:141571/53‑1 (MQ=255)                                                                                                                                       ttGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTgg                                                                                                                >  1:141571/1‑53 (MQ=255)                                                                                                                                          ccGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACgg               >  6:345464/1‑147 (MQ=255)                                                                                                                                           cGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACgg               >  2:286806/1‑146 (MQ=255)                                                                                                                                            gATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGAc                 >  8:208745/1‑143 (MQ=255)                                                                                                                                            gATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGAc                 <  7:208745/143‑1 (MQ=255)                                                                                                                                              acatGCCCAAGCGC‑TGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGAc                                                                      >  8:22021/4‑88 (MQ=255)                                                                                                                                              acatGCCCAAGCGC‑TGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGAc                                                                      <  7:22021/85‑1 (MQ=255)                                                                                                                                              acatGCCCAAGCGC‑TGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCACGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGc         >  4:131591/4‑149 (MQ=255)                                                                                                                                              acatGCCCAAGCGC‑TGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGACCTCAAGCGGGTGACGGAAGAGc         >  4:126222/4‑149 (MQ=255)                                                                                                                                              acatGCCCAAGCGC‑TGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGc         >  4:29716/4‑149 (MQ=255)                                                                                                                                                ttGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCga       <  4:240880/149‑1 (MQ=255)                                                                                                                                                  gCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCgag      <  1:125899/148‑1 (MQ=255)                                                                                                                                                  gCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGACTGGATAAAGCTCTACAGCAAGCCCCAGGATCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCgag      >  7:69949/1‑148 (MQ=255)                                                                                                                                                    cc‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACa   >  5:62764/1‑149 (MQ=255)                                                                                                                                                    cc‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACa   >  4:174154/1‑149 (MQ=255)                                                                                                                                                     c‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACAt  >  5:153741/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGCTGACGCGGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACAT  >  NC_002947/5038371‑5038666

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