Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,838,196	C→G	6.4%	K125N (AAG→AAC)	PP_1643 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,838,196	0	C	G	6.4%	290.7 / 6.2	91	K125N (AAG→AAC)	PP_1643	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base G (3/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA  >  NC_002947/1838048‑1838340                                                                                                                                                      |                                                                                                                                                 cATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAc                                                                                                                                                  >  5:312110/1‑149 (MQ=255)  aTTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACt                                                                                                                                                 >  3:48336/1‑149 (MQ=255)  aTTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTACCACA‑GTAATATCACCGCACTTTTGGCCAGTGCAGTACCTGTGGCAGCGGGCGTGCCAGCGATAAAAAAAACTGCCTTAGAACGCCATGCGCCGCTACTGCTGGGACt                                                                                                                                                 >  8:122698/1‑149 (MQ=255)     gAGCCTTGTGTTGCGCCCATTCGCGGGCATGCCCGCTCGCACA‑GTAATCCCACCGCACTTAAGGCCTGTGCAGTACCTGCGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTc                                                                                                                                               <  6:18147/148‑1 (MQ=255)      aGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCa                                                                                                                                             <  1:254434/149‑1 (MQ=255)          ttGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAgcgc                                                                                                                                         >  1:96410/1‑149 (MQ=255)             tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAGTTCCAgcatga                                                                                                                                       <  2:263098/148‑5 (MQ=255)             tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCAc                                                                                                                                       <  8:221015/148‑1 (MQ=255)             tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACg                                                                                                                                      >  4:118631/1‑149 (MQ=255)             tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGACTGTGCAGTACCTGTGGGAGCGTGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCAc                                                                                                                                       <  7:152285/148‑1 (MQ=255)              gtTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACg                                                                                                                                      >  5:310492/1‑148 (MQ=255)                 gCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTa                                                                                                                                    >  6:342143/1‑147 (MQ=255)                     ccATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAg                                                                                                                              >  7:102600/1‑149 (MQ=255)                     ccATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAg                                                                                                                              >  4:163713/1‑149 (MQ=255)                           gcgGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGc                                                                                                                        <  8:59691/149‑1 (MQ=255)                           gcgGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGc                                                                                                                        <  6:179946/149‑1 (MQ=255)                           gcgGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGc                                                                                                                        <  6:312110/149‑1 (MQ=255)                           gcgGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGACCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCgg                                                                                                                         >  6:15791/1‑148 (MQ=255)                           gcgGGCACGCCCGCTCCCACAGGTACTC‑CACTGGCCTTACGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCgg                                                                                                                         <  5:62888/148‑1 (MQ=17)                            cgGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCa                                                                                                                       <  8:247761/149‑1 (MQ=255)                              ggCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAAt                                                                                                                     >  4:261298/1‑149 (MQ=255)                              ggCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAAt                                                                                                                     <  2:114285/149‑1 (MQ=255)                                cATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAGTTCCAgcat                                                                                                                                         <  2:297380/127‑3 (MQ=255)                                cATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAGTTCCAgcat                                                                                                                                         >  1:297380/1‑125 (MQ=255)                                   gCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGaa                                                                                                                  >  4:52011/1‑147 (MQ=255)                                    cccGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAGTTCCAgcat                                                                                                                                         >  1:215644/1‑121 (MQ=255)                                    cccGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAAGTTCCAgcat                                                                                                                                         <  2:215644/123‑3 (MQ=255)                                       gCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGccc                                                                                                            >  5:299547/1‑149 (MQ=255)                                           ccACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCa                                                                                                          <  4:48336/147‑1 (MQ=255)                                           ccACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAg                                                                                                         <  3:52011/148‑1 (MQ=255)                                            caca‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGc                                                                                                       >  6:290633/1‑149 (MQ=255)                                             aca‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGaa                                                                                                                  <  6:119069/137‑1 (MQ=255)                                             aca‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGaa                                                                                                                  >  5:119069/1‑137 (MQ=255)                                             aca‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                                                                                                      >  5:366613/1‑149 (MQ=255)                                               a‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCa                                                                                                    >  4:307149/1‑149 (MQ=255)                                                 gTAATCGCACCGCAGTTAAGGCCTGAGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTCGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAc                                                                                                   >  5:306882/1‑149 (MQ=255)                                                 gTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAc                                                                                                   <  8:102600/149‑1 (MQ=255)                                                  tAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAcc                                                                                                  >  8:200632/1‑149 (MQ=255)                                                   aaTCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCa                                                                                                 >  8:14446/1‑149 (MQ=255)                                                     tCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAg                                                                                                >  7:172121/1‑148 (MQ=255)                                                        caccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACt                                                                                            >  1:201011/1‑149 (MQ=255)                                                          ccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                                                                                                      <  6:18486/121‑1 (MQ=255)                                                          ccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                                                                                                      >  5:18486/1‑121 (MQ=255)                                                           cgcacTTAAGGCCTGTGCAGTACCTGGGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCa                                                                                         <  2:257907/149‑1 (MQ=255)                                                            gcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCaa                                                                                        <  8:172121/149‑1 (MQ=255)                                                              acTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAGCAGCCCCAGGCTCACCAGCACTTCAagc                                                                                      >  7:198505/1‑149 (MQ=255)                                                              acTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagc                                                                                      <  3:261298/149‑1 (MQ=255)                                                                   aGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccagc                                                                                  <  6:310492/148‑1 (MQ=255)                                                                    ggCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCg                                                                                <  3:118631/149‑1 (MQ=255)                                                                     gCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCg                                                                                >  6:299564/1‑148 (MQ=255)                                                                         gtgCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAc                                                                           >  3:67379/1‑149 (MQ=255)                                                                         gtgCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAc                                                                           >  2:122430/1‑149 (MQ=255)                                                                           gCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccg                                                                         <  3:11618/149‑1 (MQ=255)                                                                               tACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAcc                                                                          >  4:263058/1‑144 (MQ=255)                                                                               tACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAcc                                                                          <  3:263058/144‑1 (MQ=255)                                                                                     tgGGAGCTGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCgg                                                               <  7:200632/149‑1 (MQ=255)                                                                                     tgGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCTCCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCgg                                                               <  3:163713/149‑1 (MQ=255)                                                                                       ggAGCGTGCGTGCCCTCGAAAAATCTTACTGACTTCGCCCGACATGCGCCTCTCATGCTTGAACTTCGAGCGCTCGTACAGTAGCCCGGCACTGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCACGCGCCGGGt                                                             <  6:299547/149‑1 (MQ=255)                                                                                                gTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAAcgccg                                                    <  6:366613/149‑1 (MQ=255)                                                                                                    ccGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccag                                                                                   <  5:60475/114‑1 (MQ=255)                                                                                                    ccGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccag                                                                                   >  6:60475/1‑114 (MQ=255)                                                                                                     cgcgAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGc                                               <  5:15791/149‑1 (MQ=255)                                                                                                      gcgAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCTAACGCCGCCAGCa                                              >  4:68283/1‑149 (MQ=255)                                                                                                              tCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCtt                                       >  6:139133/1‑148 (MQ=255)                                                                                                                  aCTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAGCAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTtgat                                    <  8:198505/147‑1 (MQ=255)                                                                                                                  aCTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTtga                                     <  7:14446/146‑1 (MQ=255)                                                                                                                  aCTGCCTTCGCCCGCCATGCGCAGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCAGGCAATGCACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTtgatg                                   <  4:297030/148‑1 (MQ=255)                                                                                                                    tGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCTCGTCGAACGCCGCCAGCACGCCCTTGATGaaa                                >  4:375034/1‑149 (MQ=255)                                                                                                                      ccTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGaaa                                <  3:307149/147‑1 (MQ=255)                                                                                                                       cTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGaaa                                >  1:3878/1‑146 (MQ=255)                                                                                                                       cTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGaaa                                <  2:3878/146‑1 (MQ=255)                                                                                                                          cgcccgccATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGata                          >  6:215017/1‑149 (MQ=255)                                                                                                                          cgcccgccATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGat                           <  5:290633/148‑1 (MQ=255)                                                                                                                           gcccgccATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACGAGCAATTCTAGCCTGCAGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGata                          <  5:272267/148‑1 (MQ=255)                                                                                                                              cgccATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCAAACGCCGCCAGCACGCCCTTGATGAAATAGATATTca                      >  5:112812/1‑149 (MQ=255)                                                                                                                                ccATGCGCCGCTCATGCTGGAACATCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTcacc                    <  4:219147/149‑1 (MQ=255)                                                                                                                                 cATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTcacca                   >  1:23407/1‑149 (MQ=255)                                                                                                                                      gccgcTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGaa               <  7:52693/148‑1 (MQ=255)                                                                                                                                        cgcTCATGCTGGTACTTCCAGCTCACGTACAGTAGCCCTGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCTAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCa            <  7:122698/149‑1 (MQ=255)                                                                                                                                        cgcTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGc             >  3:105022/1‑148 (MQ=255)                                                                                                                                        cgcTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGc             <  4:105022/148‑1 (MQ=255)                                                                                                                                         gcTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAg           <  5:215017/149‑1 (MQ=255)                                                                                                                                          cTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAgg          >  1:361670/1‑149 (MQ=255)                                                                                                                                          cTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAgg          >  1:287670/1‑149 (MQ=255)                                                                                                                                           gtgcgCTGGAAGTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCCCGCCCTTGATGAAATAGATATTCACCACGAAGCAGGt         >  5:33251/5‑149 (MQ=255)                                                                                                                                              tGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTc        >  1:254495/1‑147 (MQ=255)                                                                                                                                                cTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCATGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAg     <  2:287670/148‑1 (MQ=255)                                                                                                                                                cTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAgg    >  5:306636/1‑149 (MQ=255)                                                                                                                                                cTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAgg    <  2:96410/149‑1 (MQ=255)                                                                                                                                                cTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAg     <  2:361670/148‑1 (MQ=255)                                                                                                                                                cTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAg     <  2:23407/148‑1 (MQ=255)                                                                                                                                                  ggAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCa  >  6:147949/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 CATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACA‑GTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA  >  NC_002947/1838048‑1838340

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