Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,599,303	G→C	9.2%	P169R (CCA→CGA)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,599,303	0	G	C	9.2%	261.3 / 10.9	84	P169R (CCA→CGA)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base C (3/5);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATT  >  NC_002947/4599155‑4599448                                                                                                                                                      |                                                                                                                                                  tCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACAGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTg                                                                                                                                                   <  5:350673/149‑1 (MQ=255)   cAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggc                                                                                                                                                 <  6:346826/149‑1 (MQ=255)   cAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTgg                                                                                                                                                  <  6:50436/148‑1 (MQ=255)    agccGGAACCGAAGACTCAACGTCTGAGGAGTCACTGCCCAGCGACGCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                                                <  8:343117/146‑1 (MQ=255)    aGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                                                >  6:270398/1‑149 (MQ=255)    aGACGGTACCGAAGACTCAACGACAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcc                                                                                                                                                >  8:35316/1‑149 (MQ=255)     gACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                               <  3:422620/149‑1 (MQ=255)       cGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctac                                                                                                                                             >  3:77446/1‑146 (MQ=255)        ggTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctacc                                                                                                                                            <  6:411848/149‑5 (MQ=255)         gTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctaccc                                                                                                                                           <  6:168473/149‑6 (MQ=255)          tACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGcccc                                                                                                                                          >  2:358462/1‑149 (MQ=255)            ccGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGt                                                                                                                                        >  6:404695/1‑149 (MQ=255)             cGAAGACTCAACGTTAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                       >  8:83633/1‑149 (MQ=255)             cGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                       >  1:411301/1‑149 (MQ=255)              gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGCTAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgg                                                                                                                                      <  4:282044/149‑2 (MQ=37)              gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGc                                                                                                                                      <  4:391401/149‑1 (MQ=255)               aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                                     >  6:374378/1‑149 (MQ=255)               aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTggg                                                                                                                                     <  4:77446/149‑3 (MQ=38)                aGACTCAACGCCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgggg                                                                                                                                    <  7:35316/149‑4 (MQ=37)                  aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                  >  8:8978/1‑149 (MQ=255)                  aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                  >  8:2671/1‑149 (MQ=255)                  aCTCAACGTCAGGGGAGGCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                  <  3:151650/149‑1 (MQ=255)                   cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcg                                                                                                                                 <  2:35839/149‑1 (MQ=255)                    tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                >  5:201125/1‑149 (MQ=255)                    tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                >  7:20697/1‑149 (MQ=255)                    tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                          >  7:191986/1‑139 (MQ=38)                    tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                          <  8:191986/139‑1 (MQ=38)                     cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                               <  7:2671/149‑1 (MQ=255)                     cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                               <  7:8978/149‑1 (MQ=255)                      aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                <  6:378962/147‑1 (MQ=255)                      aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                              >  5:52232/1‑149 (MQ=255)                      aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                              >  7:105538/1‑149 (MQ=255)                      aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                              >  5:332642/1‑149 (MQ=255)                                gAGTCACTGCCCAGCGAATCTGTGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                                    >  6:12035/1‑149 (MQ=255)                                gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                                     <  6:81427/148‑1 (MQ=255)                                   tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                 <  5:54930/149‑1 (MQ=255)                                   tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                 >  8:51009/1‑149 (MQ=255)                                       tGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                             >  7:214131/1‑149 (MQ=255)                                         cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                             >  7:377068/1‑147 (MQ=255)                                           cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                          >  7:132207/1‑148 (MQ=255)                                           cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACAAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                          >  7:339325/1‑148 (MQ=255)                                            aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                         >  5:122338/1‑148 (MQ=255)                                             gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                       <  6:122338/149‑1 (MQ=255)                                                 aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                                   <  2:411301/149‑1 (MQ=255)                                                      cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                                                  <  5:441286/145‑1 (MQ=255)                                                      cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATc                                                                                               <  6:201125/148‑1 (MQ=255)                                                         ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                           >  2:69170/1‑149 (MQ=255)                                                         ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                           <  6:52232/149‑1 (MQ=255)                                                         ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                           <  5:12035/149‑1 (MQ=255)                                                          gCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                           >  8:31085/1‑148 (MQ=255)                                                           ccGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCa                                                                                         <  6:332642/149‑1 (MQ=255)                                                            cGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAg                                                                                        <  7:51009/149‑1 (MQ=255)                                                              ttCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTc                                                                                      <  7:83633/149‑1 (MQ=255)                                                                  tgcagATAAACGCACCGGGGGGGTCGATGATCCCCGCATCATTGGTGATACACCCCTCCCTCGCGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                                    <  1:358462/144‑1 (MQ=255)                                                                       aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                <  6:234841/130‑1 (MQ=255)                                                                       aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                >  5:234841/1‑130 (MQ=255)                                                                          aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                               >  3:321757/1‑80 (MQ=255)                                                                          aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                               <  4:321757/80‑1 (MQ=255)                                                                          aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgctgct                                                                          <  3:212510/149‑1 (MQ=255)                                                                           aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGctc                                                                         <  1:231808/149‑1 (MQ=255)                                                                                 cgggcgggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                                   <  7:31085/149‑1 (MQ=255)                                                                                       ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                             <  6:325981/149‑1 (MQ=255)                                                                                       ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                             <  5:374378/149‑1 (MQ=255)                                                                                       ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                             >  4:425101/1‑149 (MQ=255)                                                                                             tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                       <  8:20697/149‑1 (MQ=255)                                                                                             tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                       <  2:6542/149‑1 (MQ=255)                                                                                              gatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                                      >  6:313327/1‑149 (MQ=255)                                                                                                tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                          >  1:401858/1‑143 (MQ=255)                                                                                                tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                          <  2:401858/143‑1 (MQ=255)                                                                                                 cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                         >  2:44751/1‑95 (MQ=255)                                                                                                 cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                         <  1:44751/95‑1 (MQ=255)                                                                                                  aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTATGCGTGCGCCTTCGCTGGGCGGGCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                                   <  2:299253/148‑1 (MQ=255)                                                                                                  aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAg                                                  >  4:247087/1‑149 (MQ=255)                                                                                                        tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATGAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                            <  8:48266/149‑1 (MQ=255)                                                                                                        tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                            >  7:257281/1‑149 (MQ=255)                                                                                                         catTGCTGATCCACCAC‑CCCAGGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGGCAACGAGAACCTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGCgg                                           >  8:386002/1‑149 (MQ=255)                                                                                                          atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAgg                                           <  3:421283/148‑1 (MQ=255)                                                                                                          atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGt                                          <  8:257281/149‑1 (MQ=255)                                                                                                              cTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                                      >  7:300901/1‑149 (MQ=255)                                                                                                                      ccac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc                              <  5:270398/149‑1 (MQ=255)                                                                                                                       cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                             >  7:352522/1‑149 (MQ=255)                                                                                                                       cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                             >  4:430497/1‑149 (MQ=255)                                                                                                                        ac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCgg                            >  6:135943/1‑149 (MQ=255)                                                                                                                         c‑cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATATAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGt                           <  4:107898/149‑1 (MQ=255)                                                                                                                           cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                          <  8:132207/149‑1 (MQ=255)                                                                                                                            gcACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                          >  1:99909/2‑148 (MQ=38)                                                                                                                                ggggccggCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                     >  4:214933/7‑149 (MQ=38)                                                                                                                                ggggTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                     <  8:339325/149‑1 (MQ=255)                                                                                                                                      ggCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGAc               <  3:247087/149‑1 (MQ=255)                                                                                                                                              ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg       >  5:286600/1‑149 (MQ=255)                                                                                                                                              ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg       >  2:91759/1‑149 (MQ=255)                                                                                                                                              caggcgcTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg       <  2:99909/146‑1 (MQ=255)                                                                                                                                               cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCgg      >  5:169797/1‑149 (MQ=255)                                                                                                                                                   gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTTCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt  <  8:214131/149‑1 (MQ=255)                                                                                                                                                   gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGGCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt  >  7:274543/1‑149 (MQ=255)                                                                                                                                                   gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt  >  5:167623/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  TCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATT  >  NC_002947/4599155‑4599448

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