Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,206,129	T→C	6.2%	F478S (TTC→TCC)	PP_2811 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,206,129	0	T	C	6.2%	187.5 / 2.2	64	F478S (TTC→TCC)	PP_2811	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCG  >  NC_002947/3205981‑3206275                                                                                                                                                     |                                                                                                                                                   cTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtctt                                                                                                                                                    <  7:16269/149‑1 (MQ=255)       ccAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCT‑GCAGGTCACGCCCGGCGTCTTCTACTcg                                                                                                                                             >  7:80696/1‑149 (MQ=255)         aGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTcgc                                                                                                                                            <  6:327576/149‑1 (MQ=255)         aGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTcgc                                                                                                                                            <  5:127404/149‑1 (MQ=255)                aaTGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTGACGCCCGGCGTCTTCTACTCGcgctcgc                                                                                                                                      <  6:336465/148‑1 (MQ=255)                aaTGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGcgctcgct                                                                                                                                     >  5:150241/1‑149 (MQ=255)                    ccTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTgg                                                                                                                                  <  1:162188/148‑1 (MQ=255)                       ggCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggcggg                                                                                                                              <  7:286369/149‑1 (MQ=255)                          aaCTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTCCTACTCGCGCTCGCTTggcgggcg                                                                                                                            >  5:25863/1‑148 (MQ=255)                          aaCTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTCCTACTCGCGCTCGCTTGGCGGGCGa                                                                                                                           >  5:20554/1‑149 (MQ=255)                           aCTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGaa                                                                                                                          >  5:130756/1‑149 (MQ=255)                               gATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGc                                                                                                                       <  6:9720/148‑1 (MQ=255)                                 tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtcttct                                                                                                                                                  >  3:92925/1‑119 (MQ=255)                                 tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtcttct                                                                                                                                                  <  4:92925/119‑1 (MQ=255)                                        aCAGCGCCCAAGTCCCTTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAg                                                                                                             <  2:307253/149‑1 (MQ=255)                                        aCAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAg                                                                                                             <  3:187392/149‑1 (MQ=255)                                        aCAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCT‑GCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAgc                                                                                                            <  8:80696/149‑1 (MQ=255)                                         cAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAgc                                                                                                            <  1:110138/149‑1 (MQ=255)                                         cAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAgc                                                                                                            <  1:105130/149‑1 (MQ=255)                                              cccAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCAc                                                                                                       >  3:354899/1‑149 (MQ=255)                                                 aaGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCtt                                                                                                    <  7:52481/149‑1 (MQ=255)                                                  aGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTCCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTc                                                                                                   <  6:25863/149‑1 (MQ=255)                                                  aGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTCCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTc                                                                                                   <  6:20554/149‑1 (MQ=255)                                                    tCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCAc                                                                                                 >  3:213755/1‑149 (MQ=255)                                                      ccGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggcgggc                                                                                                                             >  3:146/1‑119 (MQ=255)                                                      ccGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggcgggc                                                                                                                             <  4:146/119‑1 (MQ=255)                                                      ccGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCAcc                                                                                                >  6:157218/1‑148 (MQ=255)                                                           gCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGAc                                                                                          <  6:150241/149‑1 (MQ=255)                                                             gACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGAcgc                                                                                        >  8:348903/1‑149 (MQ=255)                                                               cAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAg                                                                                            <  8:286068/143‑1 (MQ=255)                                                               cAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAg                                                                                            >  7:286068/1‑143 (MQ=255)                                                                  gACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagca                                                                                   <  7:348903/149‑1 (MQ=255)                                                                   aCCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagcag                                                                                  <  4:354899/149‑1 (MQ=255)                                                                    ccTCGTCGGGCATCAACCTTGACTTCAGCCTGACCTATGACGGCAAGTTGATGCCGTGCTTGCAGGTCACGCCCTGCGGCTTCTACTCGCGCTCGCTTGGCGGTCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagcagc                                                                                 <  5:379226/149‑1 (MQ=255)                                                                    ccTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagcag                                                                                  >  8:324263/1‑148 (MQ=255)                                                                         tcgGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGaa                                                                            >  5:277908/1‑149 (MQ=255)                                                                            ggCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCt                                                                         <  2:180131/149‑1 (MQ=255)                                                                             gCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTg                                                                        >  8:257800/1‑149 (MQ=255)                                                                              cATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGCACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTg                                                                        >  7:331517/1‑148 (MQ=255)                                                                               aTCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAAc                                                                           >  2:263417/1‑144 (MQ=255)                                                                               aTCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAAc                                                                           <  1:263417/144‑1 (MQ=255)                                                                               aTCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTa                                                                      >  1:76206/1‑149 (MQ=255)                                                                               aTCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTa                                                                      >  5:68575/1‑149 (MQ=255)                                                                               aTCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGATCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTa                                                                      >  5:2418/1‑149 (MQ=255)                                                                                   aCCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACACCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTg                                                                  <  1:60812/149‑1 (MQ=255)                                                                                        gACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAAt                                                               <  7:79505/147‑1 (MQ=255)                                                                                            tCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCG‑AGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTAc                                                        >  2:293138/1‑149 (MQ=255)                                                                                             cAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTa                                                         <  7:317063/148‑1 (MQ=255)                                                                                              aGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTAc                                                        >  2:188493/1‑148 (MQ=255)                                                                                               gCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCG‑AGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGaa                                                     <  1:293138/149‑1 (MQ=255)                                                                                                ccTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGa                                                      >  2:26667/1‑148 (MQ=255)                                                                                                  tGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAAc                                                   >  2:229547/1‑149 (MQ=255)                                                                                                   gACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAAcc                                                  >  7:343277/1‑149 (MQ=255)                                                                                                      cTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCgg                                               <  6:277908/149‑1 (MQ=255)                                                                                                        aTGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGag                                                                                     >  2:71852/1‑109 (MQ=255)                                                                                                         tGACGGCACGTTGCTGCCGGGCTGGCAGATCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGc                                              <  6:2418/147‑1 (MQ=255)                                                                                                               cACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTgg                                       >  2:230921/1‑148 (MQ=255)                                                                                                                       tgcCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGACCCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTcc                              >  1:181959/1‑149 (MQ=255)                                                                                                                       tgcCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTcc                              >  2:118393/1‑149 (MQ=255)                                                                                                                       tgcCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTcc                              >  7:8301/1‑149 (MQ=255)                                                                                                                        gcCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTccc                             >  5:51043/1‑149 (MQ=255)                                                                                                                         ccGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTcc                              <  1:230921/147‑1 (MQ=255)                                                                                                                                   aGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTAcgcg                  <  1:26667/149‑1 (MQ=255)                                                                                                                                    ggTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTAcgcg                  >  5:159674/1‑148 (MQ=255)                                                                                                                                        aCGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGt              <  8:8301/148‑1 (MQ=255)                                                                                                                                            ccGGCGTCTTCTACTCGCGCTCGCTTGGCCGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCTGATTTCCCTCAACTACGCGAAGTTCATg         >  1:119219/1‑149 (MQ=255)                                                                                                                                                 gtcgtcTTCTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGggcgg    <  8:331517/149‑1 (MQ=255)                                                                                                                                                 gTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCAGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTAATGggcgg    >  1:382885/1‑149 (MQ=255)                                                                                                                                                    ttctACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGgcggcg  >  2:6867/1‑148 (MQ=255)                                                                                                                                                    ttctACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGgcggcg  >  8:76366/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTTCAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCG  >  NC_002947/3205981‑3206275

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