Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,038,524	Δ1 bp	9.7%	coding (70/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,524	0	T	.	9.7%	189.6 / 14.4	62	coding (70/255 nt)	PP_4444	transposase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base . (4/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACATCTTA  >  NC_002947/5038380‑5038670                                                                                                                                                   |                                                                                                                                                   ggTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTggg                                                                                                                                                >  3:291773/1‑149 (MQ=255)  gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTggg                                                                                                                                                >  2:284978/1‑148 (MQ=255)  gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                               <  4:350375/149‑1 (MQ=255)  gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                               <  1:141371/149‑1 (MQ=255)  gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                               <  7:281293/149‑1 (MQ=255)  gTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGc                                                                                                                                               <  4:291773/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                                                                                                                                       <  7:187027/146‑1 (MQ=255)             aGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGt                                                                                                                                       >  8:187027/1‑146 (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                                                                                                                           >  7:276163/1‑149 (MQ=255)                      gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                           >  8:208647/1‑149 (MQ=255)                      gTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTa                                                                                                                           <  8:276163/149‑1 (MQ=255)                       tGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTAc                                                                                                                          <  2:191788/149‑1 (MQ=255)                       tGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTAc                                                                                                                          >  4:268893/1‑149 (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‑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                                                                                                                    >  1:141571/1‑53 (MQ=255)                                                                                                                              ttGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTgg                                                                                                                    <  2:141571/53‑1 (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)                                                                                                                                               gcgcTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACATCt    >  3:102587/1‑149 (MQ=255)                                                                                                                                                   ttGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACATCTTa  >  6:313593/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GGTTGCCGACACAGCGGGCCAGTGAGATCGGGCAACTGCTGCCGCATCAGTGGGTATCGGCCTGACTTACGCAAGGTGCGTTGGGCGGACGCTTACAAGCAAGTGACCG‑AGAAAGGCAAACCTGTTGCCGATGTTGCCC‑AGCGCTTGGGCATGTCCGTGCACAGTCTTTACGCCTGGATCAAGCTCTACAGCAAGCCCCAGGAGCAGCGTCAACAAGACGATGATCAGCAGGCCGAACTGCGCAAGCTGCGCGCTGAACTCAAGCGGGTGACGGAAGAGCGAGACATCTTA  >  NC_002947/5038380‑5038670

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