Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,599,297	T→G	11.8%	Q171P (CAG→CCG) ‡	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,599,297	0	T	G	11.8%	253.0 / 24.5	93	Q171P (CAG→CCG) ‡	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base G (5/6);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGG  >  NC_002947/4599150‑4599444                                                                                                                                                     |                                                                                                                                                    tGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTg                                                                                                                                                    >  4:421283/1‑149 (MQ=255)      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:2671/1‑149 (MQ=255)                       aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                              >  8:8978/1‑149 (MQ=255)                       aCTCAACGTCAGGGGAGGCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                              <  3:151650/149‑1 (MQ=255)                        cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcg                                                                                                                             <  2:35839/149‑1 (MQ=255)                         tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                            >  7:20697/1‑149 (MQ=255)                         tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                            >  5:201125/1‑149 (MQ=255)                         tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                      <  8:191986/139‑1 (MQ=38)                         tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                      >  7:191986/1‑139 (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                                                                                                                          >  5:332642/1‑149 (MQ=255)                           aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                          >  7:105538/1‑149 (MQ=255)                                     gAGTCACTGCCCAGCGAATCTGTGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                                >  6:12035/1‑149 (MQ=255)                                     gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                                 <  6:81427/148‑1 (MQ=255)                                        tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                             >  8:51009/1‑149 (MQ=255)                                        tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                             <  5:54930/149‑1 (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                                                                                       <  6:52232/149‑1 (MQ=255)                                                              ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                       <  5:12035/149‑1 (MQ=255)                                                              ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                       >  2:69170/1‑149 (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                                                                                            >  5:234841/1‑130 (MQ=255)                                                                            aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                            <  6:234841/130‑1 (MQ=255)                                                                               aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                           <  4:321757/80‑1 (MQ=255)                                                                               aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                           >  3:321757/1‑80 (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                                                         <  5:374378/149‑1 (MQ=255)                                                                                            ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                         <  6:325981/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                                                                                                     <  1:44751/95‑1 (MQ=255)                                                                                                      cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                     >  2:44751/1‑95 (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                         >  4:430497/1‑149 (MQ=255)                                                                                                                            cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                         >  7:352522/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)                                                                                                                                                    cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCgg  >  5:169797/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGG  >  NC_002947/4599150‑4599444

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