Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,649,985	T→A	5.7%	K310M (AAG→ATG) ‡	PP_4113 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,649,985	0	T	A	5.7%	268.6 / 4.5	87	K310M (AAG→ATG) ‡	PP_4113	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base A (2/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

ACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4649838‑4650131                                                                                                                                                     |                                                                                                                                                     aCATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCtt                                                                                                                                                     >  7:323882/1‑149 (MQ=255)   aTCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGc                                                                                                                                                   <  4:54645/149‑1 (MQ=255)      gcgcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGccc                                                                                                                                                 <  8:307864/148‑1 (MQ=255)      gcgcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                  >  3:410343/1‑147 (MQ=255)           cgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCGGGGCTGGGCAAGCACAGcc‑ct                                                                                                                                          <  2:9454/149‑2 (MQ=255)             ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCgg                                                                                                                                        <  4:380932/149‑1 (MQ=255)             ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGGGCAAGCACAGCC‑CTgg                                                                                                                                        <  7:20345/149‑1 (MQ=255)                aCCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCa                                                                                                                                      >  7:367449/1‑148 (MQ=255)                 cccGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCaaa                                                                                                                                    <  1:323466/149‑1 (MQ=255)                   ccgTGAATGCCCCAGCTCCCGTTCCGCACCCGCACGCCATCTTCGACCAGGTGCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCaaaaa                                                                                                                                  <  5:327478/147‑1 (MQ=255)                   cGGGGAATGCCCCAGCTCCAGTTCAGCAGCGGCACGCCATCTTCGACGAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACG‑ATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCC‑CTGGcttcga                                                                                                                                  <  6:340467/149‑6 (MQ=255)                        aaTGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCAATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑cc                                                                                                                              <  1:95420/148‑1 (MQ=255)                          tGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCa                                                                                                                          >  4:355624/1‑149 (MQ=255)                             cccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCAccc                                                                                                                       >  2:76907/1‑149 (MQ=255)                             cccAGATACAGTTCAGCACCCGCAAGCCATCTTCGACCAGGTTCACCGAGGCAGAGATGTTGGAGGTGATACCGGCATCCGAGTTGAGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTTCTTGCCCAGCC‑ACGGCAAAAAG‑CCCTCAcc                                                                                                                        >  2:231031/1‑148 (MQ=255)                              ccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCAACGAGGAAGCGATGTTGGCGGTGATACCGGCATCCGATTTGCGCTCTACC‑ATCACCTCGAAGCCCGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑ACCTCACACt                                                                                                                      >  1:117543/1‑149 (MQ=255)                                    ccAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAg                                                                                                                >  7:191943/1‑149 (MQ=255)                                     cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGc                                                                                                               >  2:416155/1‑149 (MQ=255)                                      aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCg                                                                                                              <  8:191943/149‑1 (MQ=255)                                      aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCg                                                                                                              <  7:376606/149‑1 (MQ=255)                                          cAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                          >  3:364908/1‑149 (MQ=255)                                           aGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAAt                                                                                                         <  3:152778/149‑1 (MQ=255)                                           aGCACCCGCACGCCATCTTCGACCAGCTTCACCGAGTCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAC‑CCCTCACCCTGCGCAGCgagact                                                                                                         >  5:62293/1‑147 (MQ=255)                                             cACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACA‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGc                                                                                                       <  5:367779/149‑1 (MQ=255)                                                ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                    >  2:306037/1‑149 (MQ=255)                                                ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                    >  5:291365/1‑149 (MQ=255)                                                       ccATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                              <  3:175014/148‑1 (MQ=255)                                                       ccATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                              >  4:175014/1‑148 (MQ=255)                                                        cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATgg                                                                                            >  4:409273/1‑149 (MQ=255)                                                        cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATgg                                                                                            >  2:326656/1‑149 (MQ=255)                                                         aTCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGc                                                                                           >  7:280068/1‑149 (MQ=255)                                                            ttCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGc                                                                                        >  5:156240/1‑149 (MQ=255)                                                                aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTg                                                                                    <  4:4913/149‑1 (MQ=255)                                                                 ccAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGATGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                    <  8:32795/132‑1 (MQ=255)                                                                 ccAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGATGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                    >  7:32795/1‑132 (MQ=255)                                                                  cAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGc                                                                                  >  5:116838/1‑149 (MQ=255)                                                                     gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTg                                                                               >  3:292721/1‑149 (MQ=255)                                                                      ttCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGTGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGc                                                                              <  8:123857/149‑1 (MQ=255)                                                                          ccGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAGCCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCgg                                                                           <  1:326656/149‑1 (MQ=255)                                                                            gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  5:89724/149‑1 (MQ=255)                                                                            gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  1:76907/149‑1 (MQ=255)                                                                            gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                        <  3:355624/149‑1 (MQ=255)                                                                                cAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                    <  4:410343/149‑1 (MQ=255)                                                                                 aGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAc                                                                   <  6:158692/149‑1 (MQ=255)                                                                                 aGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAc                                                                   >  8:395881/1‑149 (MQ=255)                                                                                   cGATGTTGGCGGTGATACCCGCATCCGAGTTGCGCACTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGAGACATGGCTGCCGTGGCTTGCGGGGCGGTTACCg                                                                 >  4:258610/1‑149 (MQ=255)                                                                                      tGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGACGTGGCTTGCCGGCCGGTTACCGTcg                                                              >  6:425761/1‑149 (MQ=255)                                                                                           gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                          <  5:192069/100‑1 (MQ=255)                                                                                           gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                          >  6:192069/1‑100 (MQ=255)                                                                                           gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                         <  6:291365/149‑1 (MQ=255)                                                                                            cGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTTCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCg                                                        <  7:259409/149‑1 (MQ=255)                                                                                                 aTACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTaga                                                    <  1:416155/148‑1 (MQ=255)                                                                                                  tACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTagag                                                   <  6:286739/148‑1 (MQ=255)                                                                                                     cGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                               >  3:48131/1‑149 (MQ=255)                                                                                                       gCATCCGAGTTGGGCTCTACC‑ATCACATCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAT‑CCAGCACCCTGGGGAGCGAGCTTGCCTGCGACATGGCTGCCGTGGCTTGCGGGGCGGTGAGCGTCGCGGGCGGAGAGGCATGt                                             >  4:237798/1‑149 (MQ=255)                                                                                                        cATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTc                                            >  5:393559/1‑149 (MQ=255)                                                                                                          tCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                          <  6:174502/149‑1 (MQ=255)                                                                                                           ccGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTgg                                         >  1:134918/1‑149 (MQ=255)                                                                                                            cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCTGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                        >  1:332212/1‑149 (MQ=255)                                                                                                            cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                        <  5:160788/149‑1 (MQ=255)                                                                                                            cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                        <  4:292721/149‑1 (MQ=255)                                                                                                                   cgcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGtta                                 >  7:107449/1‑148 (MQ=255)                                                                                                                   cgcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTc                                 <  4:48131/149‑1 (MQ=255)                                                                                                                    gcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTCCCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCg                                <  2:316234/149‑1 (MQ=255)                                                                                                                    gcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCg                                <  6:156240/149‑1 (MQ=255)                                                                                                                         aCC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTg                                    >  8:412607/1‑140 (MQ=255)                                                                                                                         aCC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTg                                    <  7:412607/140‑1 (MQ=255)                                                                                                                          cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                          >  7:318892/1‑149 (MQ=255)                                                                                                                          cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                          >  1:363172/1‑149 (MQ=255)                                                                                                                          cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcc                           <  6:393559/148‑1 (MQ=255)                                                                                                                          cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGc                            <  7:298781/147‑1 (MQ=255)                                                                                                                           c‑ATCACCTCGAAGCGAGGGCTGTTCTTGCCGAGCG‑CCGTCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGACTGCTTCATGGCTGCCGTGTCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcgcc                         <  8:367449/149‑1 (MQ=255)                                                                                                                              tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCTGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCt                       <  1:231031/149‑1 (MQ=255)                                                                                                                              tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCt                       <  8:280068/149‑1 (MQ=255)                                                                                                                                aCCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                      <  3:409273/148‑1 (MQ=255)                                                                                                                                 gctttttgCCGGGGCTGGGCAAGCACAGCCTCTGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                      >  1:259491/8‑148 (MQ=255)                                                                                                                                     gAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAAt                <  5:425761/149‑1 (MQ=255)                                                                                                                                     gAAGCCAGGGATTTGCTTGCCCAGCC‑CCGGCAAATAG‑CCCTCAGCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTTGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTGGATaa                 <  6:62293/148‑1 (MQ=255)                                                                                                                                         ccggggCTGGGCAAGCACAGCC‑CTGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt            >  2:121712/4‑149 (MQ=255)                                                                                                                                         ccAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt            >  2:117593/1‑149 (MQ=255)                                                                                                                                         ccAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGAAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCTCAGCCCCCTCGATAATGGCt            >  4:372999/1‑149 (MQ=255)                                                                                                                                           aGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTg           >  7:53254/1‑148 (MQ=255)                                                                                                                                           aGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTg           >  1:190698/1‑148 (MQ=255)                                                                                                                                             ggCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAc        >  8:252758/1‑149 (MQ=255)                                                                                                                                                  tgCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg   >  5:231727/1‑149 (MQ=255)                                                                                                                                                  tgCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg   >  1:46362/1‑149 (MQ=255)                                                                                                                                                    cTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCggg  <  5:211283/148‑1 (MQ=255)                                                                                                                                                    cTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCggg  <  2:332212/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4649838‑4650131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑