Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,054,462	A→C	6.8%	M1G (GTG→GGG) †‡	PP_0912 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,054,462	0	A	C	6.8%	265.7 / 8.4	89	M1G (GTG→GGG) †‡	PP_0912	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base C (2/4);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

TGGACGATTCGGCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCG  >  NC_002947/1054315‑1054607                                                                                                                                                    |                                                                                                                                                  tGGACGATTCGGCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGca                                                                                                                                                   >  7:70073/1‑148 (MQ=255)   gACGATTCGGCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACaa                                                                                                                                                <  5:151214/149‑1 (MQ=255)        ttCGGCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCg                                                                                                                                           >  3:353328/1‑149 (MQ=255)           ggCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGc                                                                                                                                        >  8:170113/1‑149 (MQ=255)            tccTTCGCACCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCACCAAACGAAAACGCCCGGCACAAGGCCGGGCg                                                                                                                                       >  2:39261/2‑149 (MQ=255)                 cGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGGAAACGAAAACGCCCGGCCTTGTGCCGGGCGttt                                                                                                                                    <  1:241826/147‑1 (MQ=255)                 cGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTg                                                                                                                                   <  7:299737/148‑1 (MQ=255)                  gCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCCTTGTGCCGGGCGttttc                                                                                                                                  <  4:172800/148‑3 (MQ=255)                  gCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCg                                                                                                                                 <  5:110636/149‑1 (MQ=255)                   ccccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCCTTGTGCCGGGCGTTTTCgt                                                                                                                                <  5:332071/149‑2 (MQ=255)                   ccccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCg                                                                                                                                 >  8:122757/1‑148 (MQ=255)                   ccccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCg                                                                                                                                 <  7:122757/148‑1 (MQ=255)                    cccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGt                                                                                                                               >  4:359487/1‑149 (MQ=255)                    cccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGt                                                                                                                               >  6:4687/1‑149 (MQ=255)                    cccAGACATCAAGAATCTGCAGCACTATGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAg                                                                                                                                               >  2:327445/1‑133 (MQ=255)                    cccAGACATCAAGAATCTGCAGCACTATGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAg                                                                                                                                               <  1:327445/133‑1 (MQ=255)                      cAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTg                                                                                                                              >  2:21627/1‑148 (MQ=255)                       aGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGaa                                                                                                                            <  1:160019/149‑1 (MQ=255)                           aTCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAgcgc                                                                                                                        >  6:22381/1‑149 (MQ=255)                             cAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGTAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                      >  6:23010/1‑149 (MQ=255)                              aaGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGACCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                      <  3:273218/148‑1 (MQ=255)                                gAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAcag                                                                                                                   <  4:211963/149‑1 (MQ=255)                                gAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGAATCTGCCTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAcag                                                                                                                   >  1:289462/1‑149 (MQ=255)                                  aTCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAcagca                                                                                                                 <  4:130609/149‑1 (MQ=255)                                   tCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAcagcag                                                                                                                >  5:233235/1‑149 (MQ=255)                                       cagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGCAAACGCCCGGCCTTGTGCCGGGCGTTTTCgtt                                                                                                                               <  6:68680/130‑3 (MQ=38)                                       cagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGCAAACGCCCGGCCTTGTGCCGGGCGTTTTCgtt                                                                                                                               >  5:68680/1‑128 (MQ=38)                                       cagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAgg                                                                                                               <  5:200882/146‑1 (MQ=255)                                       cagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGAt                                                                                                             <  6:233235/148‑1 (MQ=255)                                           aCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGaa                                                                                                        <  4:122454/149‑1 (MQ=255)                                              gtgCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTc                                                                                                     >  7:219308/1‑149 (MQ=255)                                               tgCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTc                                                                                                     <  7:170113/148‑1 (MQ=255)                                               tgCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCg                                                                                                    <  8:283917/149‑1 (MQ=255)                                                   tcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGAt                                                                                                >  4:4918/1‑149 (MQ=255)                                                    cctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATc                                                                                               >  4:207256/1‑149 (MQ=255)                                                       gTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAca                                                                                                                    <  7:308347/125‑1 (MQ=255)                                                       gTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAca                                                                                                                    >  8:308347/1‑125 (MQ=255)                                                        ttAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCg                                                                                             >  8:120471/1‑147 (MQ=255)                                                        ttAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCg                                                                                             >  4:750/1‑147 (MQ=255)                                                        ttAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCa                                                                                           <  8:70073/149‑1 (MQ=255)                                                             tttGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAca                                                                                                                    >  8:300129/1‑119 (MQ=255)                                                             tttGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTAca                                                                                                                    <  7:300129/119‑1 (MQ=255)                                                             tttGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTc                                                                                      >  4:145657/1‑149 (MQ=255)                                                                 aTCGAGGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCggg                                                                                  >  8:260680/1‑149 (MQ=255)                                                                   cGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                   <  2:261048/146‑1 (MQ=255)                                                                   cGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                   >  1:261048/1‑146 (MQ=255)                                                                     aaGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCa                                                                              <  5:22381/149‑1 (MQ=255)                                                                      aGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAg                                                                             <  7:120471/149‑1 (MQ=255)                                                                         ccTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGtt                                                                          >  8:311812/1‑149 (MQ=255)                                                                          cTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTg                                                                         <  3:145657/149‑1 (MQ=255)                                                                           tACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                   <  2:19043/138‑1 (MQ=255)                                                                           tACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                   >  1:19043/1‑138 (MQ=255)                                                                             cAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGc                                                                       <  2:289462/148‑1 (MQ=255)                                                                              aGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCaa                                                                     >  5:25899/1‑149 (MQ=255)                                                                                gATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCCTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATg                                                                   >  4:286287/1‑149 (MQ=255)                                                                                    aGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCa                                                               <  7:311812/149‑1 (MQ=255)                                                                                    aGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCa                                                               <  7:108827/149‑1 (MQ=255)                                                                                              tGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGc                                                     <  1:173826/149‑1 (MQ=255)                                                                                              tGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGc                                                     <  5:131691/149‑1 (MQ=255)                                                                                                cACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCa                                                    <  5:4687/148‑1 (MQ=255)                                                                                                cACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAg                                                   <  3:207256/149‑1 (MQ=255)                                                                                                    ccAGGCAGCGGGAAGCGTTTTACCCGGCATCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTg                                               >  1:24730/1‑149 (MQ=255)                                                                                                        gCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGt                                                                                                   <  5:79698/93‑1 (MQ=255)                                                                                                        gCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGt                                                                                                   >  6:79698/1‑93 (MQ=255)                                                                                                          aGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgat                                         >  1:75705/1‑149 (MQ=255)                                                                                                          aGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgat                                         <  1:39261/149‑1 (MQ=255)                                                                                                              ggAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGAc                                     <  1:21627/149‑1 (MQ=255)                                                                                                                 aGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCg                                  >  3:331819/1‑149 (MQ=255)                                                                                                                          cccGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatc                                        <  5:96970/134‑1 (MQ=255)                                                                                                                          cccGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatc                                        >  6:96970/1‑134 (MQ=255)                                                                                                                           ccGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGa                        >  3:138877/1‑149 (MQ=255)                                                                                                                            cGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGAc                       >  8:227239/1‑149 (MQ=255)                                                                                                                              gcagcaAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCa                     <  3:4918/149‑1 (MQ=255)                                                                                                                                agcaAACGAACACGCCCGGAACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAgg                   <  4:205083/149‑1 (MQ=255)                                                                                                                                agcaAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAgg                   >  3:195626/1‑149 (MQ=255)                                                                                                                                 gcaAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAggg                  <  8:219308/149‑1 (MQ=255)                                                                                                                                  caAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtc                                            >  7:314921/1‑122 (MQ=255)                                                                                                                                  caAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtc                                            <  8:314921/122‑1 (MQ=255)                                                                                                                                  caAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggc                 >  1:38258/1‑149 (MQ=255)                                                                                                                                   aaaCGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcg                <  3:359487/149‑1 (MQ=255)                                                                                                                                    aaCGAAACCGCCCGGCACAAGGCCGGGGGTTTGCGGTGAAGCGCTACAGCAGCCTCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcgg               <  3:750/149‑1 (MQ=255)                                                                                                                                      cgcaaaCGCCCGGCCTTGTGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcggc              >  7:86690/4‑148 (MQ=255)                                                                                                                                       gAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGGTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGa            <  7:260680/149‑1 (MQ=255)                                                                                                                                         aaaCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATg          >  4:206985/1‑149 (MQ=255)                                                                                                                                           aCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGc        >  3:100753/1‑149 (MQ=255)                                                                                                                                           aCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGc        >  1:209557/1‑149 (MQ=255)                                                                                                                                              cccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTc     <  4:331819/149‑1 (MQ=255)                                                                                                                                               ccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTcc    >  2:259972/1‑149 (MQ=255)                                                                                                                                               ccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTcc    <  3:286287/149‑1 (MQ=255)                                                                                                                                                 ggCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCg  <  4:59451/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TGGACGATTCGGCCTTCGCCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCG  >  NC_002947/1054315‑1054607

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