Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,945,617	Δ1 bp	12.5%	coding (394/531 nt)	PP_2578 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,945,617	0	G	.	12.5%	194.0 / 21.8	64	coding (394/531 nt)	PP_2578	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/34);  new base . (5/3);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

CCGCAGCGCCCTGCAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAACCTTGGCCT  >  NC_002947/2945472‑2945764                                                                                                                                                  |                                                                                                                                                     ccGCAGCGCCCTGCAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAg                                                                                                                                                    <  4:40996/148‑1 (MQ=255)    cAGCGCCCTGCAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑Ga                                                                                                                                               <  7:209874/149‑1 (MQ=255)        gcCCTGCAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTg                                                                                                                                           <  4:396593/149‑1 (MQ=255)              cAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGt                                                                                                                                       <  5:306187/147‑1 (MQ=255)              cAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTAc                                                                                                                                     <  7:25782/149‑1 (MQ=255)                ggACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGt                                                                                                                                   <  5:48027/149‑1 (MQ=255)                 gACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGt                                                                                                                                   >  6:5473/1‑148 (MQ=255)                    cTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGt                                                                                                                               >  3:54682/1‑149 (MQ=255)                               cccGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGActgctg                                                                                                                    <  5:10377/149‑1 (MQ=255)                                ccGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACtgctgc                                                                                                                   >  5:343843/1‑149 (MQ=255)                                 cggcggGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCg                                                                                                                  <  6:32921/149‑1 (MQ=255)                                   gcgggcggTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGa                                                                                                                <  2:185886/149‑1 (MQ=255)                                    cgggcggTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGac                                                                                                               <  5:109260/149‑1 (MQ=255)                                      ggcggTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGacaa                                                                                                             <  7:11658/149‑1 (MQ=255)                                       gcggTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGacaac                                                                                                            <  6:343843/149‑1 (MQ=255)                                                    gTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGcc                                                                                               >  1:182584/1‑149 (MQ=255)                                                      cGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTg                                                                                             >  4:7712/1‑149 (MQ=255)                                                       gCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTgg                                                                                            >  5:16785/1‑149 (MQ=255)                                                        cTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGAGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGt                                                                                           <  1:179842/149‑1 (MQ=255)                                                           cGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGACGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGaa                                                                                        >  5:224523/1‑149 (MQ=255)                                                            gggTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAAt                                                                                       <  3:19292/149‑1 (MQ=255)                                                            gggTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAAt                                                                                       <  5:5473/149‑1 (MQ=255)                                                              gTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGACGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATgc                                                                                     <  6:224523/149‑1 (MQ=255)                                                                          cGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCaa                                                                         <  4:141001/149‑1 (MQ=255)                                                                           gCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAAc                                                                        <  2:182584/149‑1 (MQ=255)                                                                             ttCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCt                                                                      >  7:297434/1‑149 (MQ=255)                                                                                gACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGGAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTgc                                                                    <  5:440700/148‑1 (MQ=255)                                                                                    cGGTACCAGGGCGCTACGCCTCGCTGGACGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAg                                                               <  5:259325/149‑1 (MQ=255)                                                                                       tACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTg                                                            <  3:396794/149‑1 (MQ=255)                                                                                         ccAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTgg                                                           <  3:33014/148‑1 (MQ=255)                                                                                           aGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGcc                                                         >  6:154958/1‑148 (MQ=255)                                                                                            gggCGCTACTCCTCGCGGGAGGCCGTGGGGCGCCTGAAGCAGCAAGGGCAGGGCGAGCT‑GACCTGCCGGACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACTTGCGCAGGTTGGcctg                                                       <  4:15919/149‑3 (MQ=255)                                                                                                       cTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCggg                                            >  3:178508/1‑149 (MQ=255)                                                                                                       cTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCggg                                            >  2:145759/1‑149 (MQ=255)                                                                                                       cTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCggg                                            >  7:26832/1‑149 (MQ=255)                                                                                                       cTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAACTTGGCCGAGGCGGTGCcgg                                            >  8:208435/1‑149 (MQ=255)                                                                                                        tCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCgggg                                           >  8:101634/1‑149 (MQ=255)                                                                                                          gCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGa                                         <  7:101634/149‑1 (MQ=255)                                                                                                            tGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGt                                       <  8:297434/149‑1 (MQ=255)                                                                                                               aGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GATCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGcc                                    <  1:145759/149‑1 (MQ=255)                                                                                                               aGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGcc                                    >  3:77852/1‑149 (MQ=255)                                                                                                                 gCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGTGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCg                                  <  5:246720/149‑1 (MQ=255)                                                                                                                        ggCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAggg                            >  8:349870/1‑148 (MQ=255)                                                                                                                         gcgcCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAggg                            >  1:257131/1‑147 (MQ=255)                                                                                                                         gcgcCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTa                          <  4:77852/149‑1 (MQ=255)                                                                                                                          cgcCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCCAGGCGGTGCGGGGAAGTCCCCGGCAGGGTAg                         >  3:356730/1‑149 (MQ=255)                                                                                                                           gcCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGc                        <  3:7712/149‑1 (MQ=255)                                                                                                                               gAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGcc                    >  3:292892/1‑149 (MQ=255)                                                                                                                                aaGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGc                     <  7:349870/147‑1 (MQ=255)                                                                                                                                 agcagcAAGGGCACGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCGGGTGAATGCGCACCAGGCCAACCTGCGCAAGTGGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCt                  <  3:334916/149‑1 (MQ=255)                                                                                                                                   cagcaAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGc                <  2:257131/149‑1 (MQ=255)                                                                                                                                   cagcaAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGCCCGAGGCGGTGCGGGCGAGTGCCCGGCAGGGTAGCTGCCCTg                 >  5:228869/1‑148 (MQ=255)                                                                                                                                    agcaAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGATAGCTGCCCTGCa               <  6:16785/149‑1 (MQ=255)                                                                                                                                     gcaAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCAt              >  1:21640/1‑149 (MQ=255)                                                                                                                                      gtacggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGt                           >  3:156802/5‑135 (MQ=255)                                                                                                                                      gtacggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGt                           <  4:156802/131‑1 (MQ=255)                                                                                                                                         cggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGc                <  2:116616/142‑1 (MQ=255)                                                                                                                                         cggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGc                >  1:116616/2‑143 (MQ=255)                                                                                                                                         cggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCa               <  7:184329/143‑1 (MQ=255)                                                                                                                                         cggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCa               >  8:184329/2‑144 (MQ=255)                                                                                                                                         cggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAAc          >  1:90542/2‑149 (MQ=255)                                                                                                                                          ggCAGGTC‑AGCTCGACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAAcc         >  7:134135/1‑149 (MQ=255)                                                                                                                                              ggTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAACCTTgg     >  8:251313/1‑149 (MQ=255)                                                                                                                                               gTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAACCTTGGc    <  2:21640/149‑1 (MQ=255)                                                                                                                                                  gAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTg                                                                                             >  7:258308/1‑57 (MQ=255)                                                                                                                                                  gAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTg                                                                                             <  8:258308/57‑1 (MQ=255)                                                                                                                                                  gAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAACCTTGGCCt  <  4:292892/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCGCAGCGCCCTGCAGGACCTGTTGGGCGGCCCGGCGGGCGGTGGCAAGGCGTCGCTGCGGGTGGATGTGCAACGCTTCGACACGGTACCAGGGCGCTACGCCTCGCTGGAGGCCGTGTGGCGCCTGAAGCAGCAAGGGCAGGTCGAGCT‑GACCTGCCGTACGTCAGTGCAGACTGCTGCGGACAACAGCATCACCAGCCTGGTGAATGCGCACCAGGCCAACCTGCGCAAGTTGGCCGAGGCGGTGCGGGGGAGTGCCCGGCAGGGTAGCTGCCCTGCATGAACCTTGGCCT  >  NC_002947/2945472‑2945764

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