Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,761,301	T→A	6.4%	L118Q (CTG→CAG)	gpmI →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,761,301	0	T	A	6.4%	291.4 / 7.0	94	L118Q (CTG→CAG)	gpmI	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base A (3/3);  total (49/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 = 8.89e-01

ATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAA  >  NC_002947/5761153‑5761449                                                                                                                                                     |                                                                                                                                                     atGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGAct                                                                                                                                                      >  8:288503/1‑149 (MQ=255) atGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACAGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGAct                                                                                                                                                      >  5:274865/1‑149 (MQ=255)   gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGGctgc                                                                                                                                                    >  5:435837/1‑149 (MQ=255)      ccTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGt                                                                                                                                                 <  6:116272/149‑1 (MQ=255)       cTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTc                                                                                                                                                >  8:13048/1‑149 (MQ=255)           gcgcAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGGCTGCTGTCCGAt                                                                                                                                            <  6:435837/149‑1 (MQ=255)            cgcAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGagg                                                                                                                                           <  7:159659/149‑3 (MQ=255)               aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                        <  3:371465/149‑1 (MQ=255)               aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                        <  2:284822/149‑1 (MQ=255)               aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGaggatg                                                                                                                                        <  4:408833/149‑6 (MQ=255)                  tCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGgtggtg                                                                                                                                     >  8:370510/1‑149 (MQ=255)                  tCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGgtggtg                                                                                                                                     <  2:365027/149‑1 (MQ=255)                  gcgcgTCGTGTATCAGGACTTCACGCGGGTCACCAAGGCCATCCGCGACGGCGAGTGCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGgtggtg                                                                                                                                     <  4:311933/148‑1 (MQ=255)                   cgcgTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGgtgt                                                                                                                                    <  3:312511/149‑1 (MQ=255)                    gcgTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGgtgt                                                                                                                                    <  7:135741/148‑1 (MQ=255)                    gcgTCGTGCATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGgtgt                                                                                                                                    <  3:184452/148‑1 (MQ=255)                     cgtcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCAGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTa                                                                                                                                   >  6:86079/1‑148 (MQ=255)                      gtcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTaca                                                                                                                                 <  5:426002/149‑1 (MQ=255)                        cgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                                                                                                                               <  3:293002/149‑1 (MQ=255)                         gtgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                                                                                                                               <  7:289725/148‑1 (MQ=255)                             aTCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                                                                                                                          >  2:26976/1‑149 (MQ=255)                              tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggt                                                                                                                                         >  5:250909/1‑133 (MQ=255)                              tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggt                                                                                                                                         <  6:250909/133‑1 (MQ=255)                              tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCATCGGACAGCAGTCCGaggat                                                                                                                                         <  6:285596/133‑5 (MQ=255)                              tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCATCGGACAGCAGTCCGaggat                                                                                                                                         >  5:285596/1‑129 (MQ=255)                               cAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                                                                                                                          >  7:68878/1‑147 (MQ=255)                                 ggACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAt                                                                                                                                            >  8:192834/1‑127 (MQ=255)                                 ggACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAt                                                                                                                                            <  7:192834/127‑1 (MQ=255)                                     ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCa                                                                                                                  >  3:69597/1‑149 (MQ=255)                                      cgtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCAc                                                                                                                 >  8:409266/4‑149 (MQ=255)                                       cACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACt                                                                                                                >  4:127718/1‑149 (MQ=255)                                         cgcgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtg                                                                                                              <  7:153122/149‑1 (MQ=255)                                          gcgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtgg                                                                                                             >  6:84076/1‑149 (MQ=255)                                           cgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggt                                                                                                            <  2:65366/149‑1 (MQ=255)                                            gggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtg                                                                                                           <  7:333652/149‑1 (MQ=255)                                            gggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtg                                                                                                           >  8:202936/1‑149 (MQ=255)                                              gTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGc                                                                                                         <  5:86079/149‑1 (MQ=255)                                                  ccAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATg                                                                                                     >  7:199305/1‑149 (MQ=255)                                                   cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATgg                                                                                                    >  2:162001/1‑149 (MQ=255)                                                      ggCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTAGCCATGGcc                                                                                                  >  1:365373/1‑148 (MQ=255)                                                        ccATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGaa                                                                                               <  2:116200/149‑1 (MQ=255)                                                           tCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTg                                                                                            >  8:277756/1‑149 (MQ=255)                                                              gcgACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTgg                                                                                           >  1:158401/1‑147 (MQ=255)                                                               cgACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccg                                                                                        <  7:409266/149‑1 (MQ=255)                                                                   ggCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCa                                                                                    >  2:165996/1‑149 (MQ=255)                                                                        gTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCTGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTg                                                                               <  6:274865/149‑1 (MQ=255)                                                                        gTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTg                                                                               <  1:26976/149‑1 (MQ=255)                                                                           cttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccg                                                                                        <  8:274826/137‑1 (MQ=255)                                                                           cttaGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGgccg                                                                                        >  7:274826/1‑137 (MQ=255)                                                                            ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTg                                                                               <  1:35580/145‑1 (MQ=255)                                                                            ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTg                                                                               >  2:35580/1‑145 (MQ=255)                                                                            ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGgcgc                                                                           >  5:309297/1‑149 (MQ=255)                                                                              cGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGTTGGTTGTGTACTCAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGgcg                                                                            <  5:301058/146‑1 (MQ=255)                                                                              cGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGgcg                                                                            <  7:202936/146‑1 (MQ=255)                                                                              cGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCt                                                                          >  2:14977/1‑148 (MQ=255)                                                                                     ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGaaaa                                                                     >  3:203617/1‑146 (MQ=255)                                                                                     ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGaaaa                                                                     <  4:203617/146‑1 (MQ=255)                                                                                     ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAGAGAt                                                                  <  7:13048/149‑1 (MQ=255)                                                                                       ggTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTgg                                                                                           <  3:246136/122‑1 (MQ=255)                                                                                       ggTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTgg                                                                                           >  4:246136/1‑122 (MQ=255)                                                                                             gACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTg                                                           >  3:311977/1‑148 (MQ=255)                                                                                             gACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTg                                                           >  3:116762/1‑148 (MQ=255)                                                                                                      tgtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTcc                                                 <  8:199305/149‑1 (MQ=255)                                                                                                       gtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCt                                                >  8:300120/1‑149 (MQ=255)                                                                                                          gACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGACGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGa                                             >  2:391413/1‑149 (MQ=255)                                                                                                           aCAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGAt                                            <  1:162001/149‑1 (MQ=255)                                                                                                           aCAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGAt                                            <  7:300120/149‑1 (MQ=255)                                                                                                               ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGcc                                        >  4:402833/1‑149 (MQ=255)                                                                                                               ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACACCTTCCt                                                <  6:296044/141‑1 (MQ=255)                                                                                                               ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACACCTTCCt                                                >  5:296044/1‑141 (MQ=255)                                                                                                                gCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCcg                                       >  7:248384/1‑149 (MQ=255)                                                                                                                 cAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCcgc                                      <  7:90253/149‑1 (MQ=255)                                                                                                                      ggTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGacac                                 >  4:415766/1‑149 (MQ=255)                                                                                                                       gTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACAcc                                >  4:191908/1‑149 (MQ=255)                                                                                                                        tGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACAccc                               >  6:85787/1‑149 (MQ=255)                                                                                                                            ggCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGcc                           >  3:65446/1‑149 (MQ=255)                                                                                                                            ggCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGc                            >  1:244811/1‑148 (MQ=255)                                                                                                                             gCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCa                          <  4:116762/149‑1 (MQ=255)                                                                                                                             gCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCa                          >  5:308694/1‑149 (MQ=255)                                                                                                                                 ggCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGCCCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGAACTACCTGCACGCCTTCCTCGATGGCCGCGCCACCCCGCCACGCa                      >  7:362928/1‑149 (MQ=255)                                                                                                                                 ggCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCcgc                                      >  7:166151/1‑133 (MQ=255)                                                                                                                                 ggCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCcgc                                      <  8:166151/133‑1 (MQ=255)                                                                                                                                  gCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAg                     >  3:58737/1‑149 (MQ=255)                                                                                                                                   ccGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAgc                    >  6:168052/1‑149 (MQ=255)                                                                                                                                     gTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAgcgc                  <  7:277756/149‑1 (MQ=255)                                                                                                                                        cacaTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCa               <  3:67115/149‑1 (MQ=255)                                                                                                                                          caccatCGGACAGCAGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgt             >  3:142006/6‑149 (MQ=255)                                                                                                                                           aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtc            <  5:262617/149‑1 (MQ=255)                                                                                                                                           aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGAAgtc            >  8:78943/1‑149 (MQ=255)                                                                                                                                             ccTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgt          >  7:258044/1‑149 (MQ=255)                                                                                                                                             catCGGACAGCAGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgt          >  6:23452/3‑149 (MQ=255)                                                                                                                                              cTCGGACTGCTGTCCGATGGTGGTGTACACCGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgtc         >  7:58640/1‑149 (MQ=255)                                                                                                                                              cTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgtc         >  1:12670/1‑149 (MQ=255)                                                                                                                                              cTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgt          <  4:69597/148‑1 (MQ=255)                                                                                                                                                cGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCa       >  2:435591/1‑149 (MQ=255)                                                                                                                                                  gACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATc     <  7:288503/149‑1 (MQ=255)                                                                                                                                                  gACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATc     <  5:84076/149‑1 (MQ=255)                                                                                                                                                     tgctgTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGaa  >  3:163949/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAA  >  NC_002947/5761153‑5761449

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