Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,898,333	Δ1 bp	5.6%	coding (211/957 nt)	PP_4306 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,898,333	0	G	.	5.6%	235.9 / 6.2	72	coding (211/957 nt)	PP_4306	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base . (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCACTTCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGG  >  NC_002947/4898187‑4898481                                                                                                                                                   |                                                                                                                                                     gATCACTTCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGc                                                                                                                                                    <  3:266499/149‑1 (MQ=255)     aCTTCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTc                                                                                                                                                >  5:192860/1‑149 (MQ=255)       ttCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCgg                                                                                                                                              >  1:205400/1‑149 (MQ=255)       ttCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCgg                                                                                                                                              >  4:98962/1‑149 (MQ=255)           cgTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGggcgg                                                                                                                                          <  6:244182/149‑1 (MQ=255)                 tGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTg                                                                                                                                                      <  7:56466/131‑1 (MQ=255)                 tGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTg                                                                                                                                                      >  8:56466/1‑131 (MQ=255)                 tGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACg                                                                                                                                    >  4:56931/1‑149 (MQ=255)                  ggCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACg                                                                                                                                    >  2:40295/1‑148 (MQ=255)                        gCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAg                                                                                                                             <  6:192860/149‑1 (MQ=255)                        gCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCACAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAg                                                                                                                             >  3:129230/1‑149 (MQ=255)                         cTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAgg                                                                                                                            <  1:15570/149‑1 (MQ=255)                             cccAGCAGCCACTGGCCAGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCGACTGCTGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGctac                                                                                                                        <  6:95569/149‑1 (MQ=255)                              ccAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGcgacc                                                                                                                       <  5:244961/149‑1 (MQ=255)                                 gcagcCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGcgaccgac                                                                                                                    <  7:40731/149‑1 (MQ=255)                                     ccttgagcaacacctgcttgccgccctcgaccctgccccacaggaaacccgccgcctgttccacggtcgcggccgctgctggcagggcctggagcaggtcaccgtcgaCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCa                                                                                                                <  3:245698‑M1/41‑1 (MQ=255)                                       ttgagcaacacctgcttgccgccctcgaccctgccccacaggaaacccgccgcctgttccacggtcgcggccgctgctggcagggcctggagcaggtcaccgtcgaCTgg                                                                                                                                                     <  3:154907‑M1/4‑1 (MQ=255)                                       ttgagcaacacctgcttgccgccctcgaccctgccccacaggaaacccgccgcctgttccacggtcgcggccgctgctggcagggcctggagcaggtcaccgtcgaCTgg                                                                                                                                                     >  4:154907‑M1/107‑110 (MQ=255)                                       aCTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACg                                                                                                              >  5:18120/1‑149 (MQ=255)                                        cTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGc                                                                                                             >  6:328097/1‑149 (MQ=255)                                           gCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCACAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCt                                                                                                          <  4:129230/149‑1 (MQ=255)                                               cacctgcttgccgccctcgaccctgccccacaggaaacccgccgcctgttccacggtcgcggccgctgctggcagggcctggagcaggtcaccgtcgaCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAg                                                                                                      >  4:245698‑M1/99‑149 (MQ=255)                                                  ctttCGGGCAGGTAGCGGTGCTGGATGTTGACTGCCTGACCGGTCCACTGCGGGCGCTCGGCCATGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCa                                                                                                   <  7:146851/146‑1 (MQ=255)                                                     cttgccgccctcgaccctgccccacaggaaacccgccgcctgttccacggtcgcggccgctgctggcagggcctggagcaggtcaccgtcgaCTGGCCTTCGGGAGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGt                                                                                                 >  6:180100‑M1/93‑148 (MQ=255)                                                        ggCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGAc                                                                                             <  3:56931/149‑1 (MQ=255)                                                             gTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGcga                                                                                                                         >  6:65754/1‑116 (MQ=255)                                                             gTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGa                                                                                        >  2:27951/1‑149 (MQ=255)                                                              tAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGa                                                                                        <  6:33278/148‑1 (MQ=255)                                                              tAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGAc                                                                                       <  1:40295/149‑1 (MQ=255)                                                                         gatgatGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGc                                                                            >  3:8191/1‑149 (MQ=255)                                                                          atgatgACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGcc                                                                           >  2:162652/1‑149 (MQ=255)                                                                            gatgaCTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCt                                                                         >  4:177938/1‑149 (MQ=255)                                                                             atgaCTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCt                                                                         <  4:262110/148‑1 (MQ=255)                                                                               gaCTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGc                                                                       >  2:29224/1‑148 (MQ=255)                                                                                  tGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcag                                                                    >  4:221822/1‑148 (MQ=255)                                                                                  tGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcag                                                                    >  8:77901/1‑148 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGGTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcagca                                                                  <  4:8191/148‑1 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCa                                                                                                       >  1:43916/1‑111 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCa                                                                                                       <  2:43916/111‑1 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcagcag                                                                 >  7:231215/1‑149 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcagca                                                                  <  5:328097/148‑1 (MQ=255)                                                                                    ccTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCcagca                                                                  <  3:221822/148‑1 (MQ=255)                                                                                      tGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCg                                                               >  4:160467/1‑149 (MQ=255)                                                                                        aCCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGc                                                             <  3:16435/149‑1 (MQ=255)                                                                                        aCCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGc                                                             >  8:105059/1‑149 (MQ=255)                                                                                          ccGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAg                                                                                                      <  8:187705/106‑1 (MQ=255)                                                                                          ccGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAg                                                                                                      >  7:187705/1‑106 (MQ=255)                                                                                          ccGTCCACTGCGAGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCcg                                                           >  4:231840/1‑149 (MQ=255)                                                                                           cGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCcgc                                                          <  8:238027/149‑1 (MQ=255)                                                                                            gTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCcgc                                                          >  4:160287/1‑148 (MQ=255)                                                                                             tCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGcc                                                                                                            <  1:231870/97‑1 (MQ=255)                                                                                             tCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGcc                                                                                                            >  2:231870/1‑97 (MQ=255)                                                                                                 cTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGt                                                    <  8:231215/149‑1 (MQ=255)                                                                                                    cGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGACCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTgg                                                  >  6:89054/1‑148 (MQ=255)                                                                                                       gcgcTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGaa                                                >  4:57159/1‑147 (MQ=255)                                                                                                         gcTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAgg                                            <  6:18120/149‑1 (MQ=255)                                                                                                         gcTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCT‑GCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCTACGGTGACCTGCTCCAGGCCATGCCAGCAGCGGCCGAGACCGTGGAACAgtc                                           >  8:117732/1‑147 (MQ=255)                                                                                                         gcTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCT‑GCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAggc                                           >  8:122205/1‑149 (MQ=255)                                                                                                              gCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCT‑GCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTTGAACAGGcggcgg                                      <  7:117732/149‑1 (MQ=255)                                                                                                              gCCAAGCTGCGCAGCATGGCTTCCAGTGCGGCGAGCT‑GCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGcggcgg                                      <  7:122205/149‑1 (MQ=255)                                                                                                               ccAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGTTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGcggcgg                                      >  8:216976/1‑149 (MQ=255)                                                                                                               ccAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGcggcgg                                      <  7:39715/149‑1 (MQ=255)                                                                                                                cAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCggg                                     >  3:99856/1‑149 (MQ=255)                                                                                                                    cTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTc                                 <  1:162652/149‑1 (MQ=255)                                                                                                                    cTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTc                                 <  3:231840/149‑1 (MQ=255)                                                                                                                      gcgcAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGCCCGTGGAACAGGCGGCGGGTTTCCt                               >  6:153479/1‑149 (MQ=255)                                                                                                                      gcgcAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCt                               >  8:239301/1‑149 (MQ=255)                                                                                                                         cagcaTGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCtgtg                            >  5:326086/1‑149 (MQ=255)                                                                                                                            caTGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTgggg                         >  1:84142/1‑149 (MQ=255)                                                                                                                                gCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAgg                     <  7:105059/149‑1 (MQ=255)                                                                                                                                gCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAgg                     <  3:177938/149‑1 (MQ=255)                                                                                                                                     cAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCg                 <  7:239301/148‑1 (MQ=255)                                                                                                                                     cAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGa                <  1:27951/149‑1 (MQ=255)                                                                                                                                        ttCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAggg             <  4:99856/149‑1 (MQ=255)                                                                                                                                            gCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGTTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGggcggc         <  7:216976/149‑1 (MQ=255)                                                                                                                                                   GGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAgg  <  7:143910‑M1/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GATCACTTCGCGTTGCTGGCAGGGCTCACCCAGCAGCCACTGGCCCGGGCTGTCGGGCAGGTAGCGGTGCTGGATGATGACTGCCTGACCCGTCCACTGCGGGCGCTCGGCCAAGCTGCGCAGCATGGCTTCCAGTTCGGCGAGCTGGCCTTCGGGCGGCTCACGGAACAAGGCGACCGACAGCACGCCCTGCAGCCAGTCGACGGTGACCTGCTCCAGGCCCTGCCAGCAGCGGCCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGG  >  NC_002947/4898187‑4898481

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