Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,725	C→T	6.1%	Q49Q (CAG→CAA)	hmgC ←	maleylacetoacetate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,725	0	C	T	6.1%	296.0 / 5.5	97	Q49Q (CAG→CAA)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base C (40/51);  new base T (3/3);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCA  >  NC_002947/5241577‑5241870                                                                                                                                                      |                                                                                                                                                      tgGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGccc                                                                                                                                                       >  2:132559/1‑149 (MQ=255)    ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                                    >  7:254779/1‑149 (MQ=255)    ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                                    <  8:254779/149‑1 (MQ=255)       ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                                 <  7:176491/149‑3 (MQ=255)       ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                                 <  8:285685/149‑3 (MQ=255)       ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCggg                                                                                                                                                 <  7:224105/149‑1 (MQ=255)         aCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGtt                                                                                                                                               >  4:148833/1‑149 (MQ=255)          ccccccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTg                                                                                                                                              >  3:122609/1‑149 (MQ=255)           cccccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGa                                                                                                                                             <  4:84902/149‑1 (MQ=255)             cccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACg                                                                                                                                           <  4:94243/149‑1 (MQ=255)             cccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACg                                                                                                                                           >  8:41498/1‑149 (MQ=255)             cccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACg                                                                                                                                           >  2:121186/1‑149 (MQ=255)             cccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGAAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCgggaactcg                                                                                                                                           >  1:332028/1‑143 (MQ=255)                 aaCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                                       <  3:314173/149‑1 (MQ=255)                  aCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                                       >  5:39053/1‑148 (MQ=255)                    ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                                     >  7:130386/1‑148 (MQ=255)                          ccGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACAAGCACCTCACCACCGTCGGTGCGCAAGGCGG‑GTACGCGGCCCTGCGGGTTGCCGGCGACATAGccgt                                                                                                                              >  4:36486/1‑147 (MQ=255)                           cGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTc                                                                                                                             >  6:235941/1‑149 (MQ=255)                            gTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCg                                                                                                                            <  4:215857/149‑1 (MQ=255)                                cTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCAc                                                                                                                        >  7:7080/1‑149 (MQ=255)                                   ggcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgc                                                                                                                      <  4:194444/148‑1 (MQ=255)                                    gcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctg                                                                                                                    <  2:184785/149‑1 (MQ=255)                                    gcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctg                                                                                                                    <  2:21057/149‑1 (MQ=255)                                    gcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctg                                                                                                                    <  8:127522/149‑1 (MQ=255)                                     cggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctgc                                                                                                                   <  6:39053/149‑1 (MQ=255)                                      ggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctgct                                                                                                                  >  2:289783/1‑149 (MQ=255)                                            cGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑t                                                                                                        >  8:214594/1‑149 (MQ=255)                                                  caACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑Tgcag                                                                                                    >  4:36654/1‑147 (MQ=255)                                                    aCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAgg                                                                                                <  1:132559/149‑1 (MQ=255)                                                    aCGCAGGCTGTGGATAAAACTCTTCCAGGTACTCGATGATCGCAGGCGACTGCACCAGAACCTCACAACCGTCGGTGCGCAAGACCG‑CTACGCGGCCCACAGGGATAACGCCTACATAGCCATCGCCACGCTGCTAA‑CCC‑‑‑TGAACCAgg                                                                                                >  7:194567/1‑149 (MQ=255)                                                              tgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGc                                                                                       <  5:318724/148‑1 (MQ=255)                                                                gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAgg                                                                                    <  8:130386/149‑1 (MQ=255)                                                                gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAgg                                                                                    >  7:219555/1‑149 (MQ=255)                                                                gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAgg                                                                                    <  5:200538/149‑1 (MQ=255)                                                                gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAgg                                                                                    >  5:170942/1‑149 (MQ=255)                                                                 aTAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAgg                                                                                    >  1:21151/1‑148 (MQ=255)                                                                  tAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGa                                                                                  >  6:135843/1‑149 (MQ=255)                                                                   aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGAc                                                                                 <  3:36654/149‑1 (MQ=255)                                                                    aaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGAc                                                                                 <  4:300889/148‑1 (MQ=255)                                                                      ccTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTgg                                                                              >  4:367233/1‑149 (MQ=255)                                                                        tctTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGTGACTGGt                                                                             >  1:326234/1‑148 (MQ=255)                                                                            ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                           <  7:214594/146‑1 (MQ=255)                                                                            ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                           <  1:177727/146‑1 (MQ=255)                                                                            ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACAACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTAA‑CCC‑‑‑TGAAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                         >  8:264276/1‑148 (MQ=255)                                                                             cAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAagcc                                                                        <  3:4721/148‑1 (MQ=255)                                                                                gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAagccagcc                                                                    <  1:370708/149‑1 (MQ=255)                                                                                gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAagccagcc                                                                    >  2:209404/1‑149 (MQ=255)                                                                                  aCTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTg                                                                                             >  7:195943/1‑122 (MQ=255)                                                                                  aCTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTg                                                                                             <  8:195943/122‑1 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                                <  4:179996/147‑1 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                                >  4:310021/1‑147 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTc                                                               <  5:235941/148‑1 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTc                                                               <  4:122609/148‑1 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCa                                                              <  5:135843/149‑1 (MQ=255)                                                                                      gatgatCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCa                                                              <  1:289783/149‑1 (MQ=255)                                                                                          atCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCCGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGc                                                          <  7:155370/149‑1 (MQ=255)                                                                                           tCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGcc                                                         <  3:36486/149‑1 (MQ=255)                                                                                             gCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAg                                                       <  1:207690/149‑1 (MQ=255)                                                                                              ccGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGc                                                          >  6:210103/1‑145 (MQ=255)                                                                                              ccGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGc                                                          <  5:210103/145‑1 (MQ=255)                                                                                              ccGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAgg                                                      >  5:52285/1‑149 (MQ=255)                                                                                               cGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAgg                                                      <  1:102283/149‑1 (MQ=255)                                                                                                 gCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGc                                                    >  1:70402/1‑148 (MQ=255)                                                                                                 gCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGc                                                    >  4:263325/1‑148 (MQ=255)                                                                                                    aCTGCCCCCCCACCTCACCCCCGCCGGTGCGCAAGGCCG‑CTACGCGGCCCCGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑CGCAGCAGCTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATg                                                <  4:289015/149‑1 (MQ=255)                                                                                                     cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCACGGCAATgc                                               <  1:198675/149‑1 (MQ=255)                                                                                                      tGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGa                                                                                  <  2:214247/113‑1 (MQ=255)                                                                                                        cACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAggg                                                     <  8:234746/140‑1 (MQ=255)                                                                                                        cACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAggg                                                     >  7:234746/1‑140 (MQ=255)                                                                                                        cACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCa                                            >  3:90428/1‑149 (MQ=255)                                                                                                         aCCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCAc                                           >  7:153032/1‑149 (MQ=255)                                                                                                             gCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCTTTCAGTGCCAGGGCAATGCGCACCCgg                                       <  1:209404/149‑1 (MQ=255)                                                                                                                ccTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                                    <  2:326234/149‑1 (MQ=255)                                                                                                                   caccacCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGAc                                 <  7:176540/149‑1 (MQ=255)                                                                                                                    accaccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGcc                                                         >  2:181443/1‑124 (MQ=255)                                                                                                                    accaccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGcc                                                         <  1:181443/124‑1 (MQ=255)                                                                                                                     ccaccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCCTCTTACAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGa                                  <  1:121186/149‑1 (MQ=255)                                                                                                                      caccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTAGATAGCCCTCGCCACGCTTCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGGCATCCCCTTGAGTGCCAGGGCAATGCGCACGCGGTAGGACGAg                              <  8:194567/149‑1 (MQ=255)                                                                                                                      caccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAg                              <  3:310021/149‑1 (MQ=255)                                                                                                                        ccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCGCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGt                            <  4:116616/149‑1 (MQ=255)                                                                                                                         cGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTg                           <  2:332028/149‑1 (MQ=255)                                                                                                                          gTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTgg                          <  2:70402/149‑1 (MQ=255)                                                                                                                             ggTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGaa                        <  2:298554/148‑1 (MQ=255)                                                                                                                               tGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  3:325784/1‑149 (MQ=255)                                                                                                                               cccgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  4:366255/3‑149 (MQ=255)                                                                                                                                ccgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  5:19237/2‑148 (MQ=255)                                                                                                                                 cgcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGt                    >  6:35632/1‑148 (MQ=255)                                                                                                                                  gcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑CGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     <  6:209728/146‑1 (MQ=255)                                                                                                                                  gcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑CGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  5:209728/1‑146 (MQ=255)                                                                                                                                    aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa                  >  8:51494/1‑147 (MQ=255)                                                                                                                                    aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa                  >  3:338239/1‑147 (MQ=255)                                                                                                                                    aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATa                <  5:35632/149‑1 (MQ=255)                                                                                                                                       gCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt             <  3:148833/149‑1 (MQ=255)                                                                                                                                        ccG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgtg            >  5:135146/1‑149 (MQ=255)                                                                                                                                         cG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  7:267125/1‑139 (MQ=255)                                                                                                                                         cG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     <  8:267125/139‑1 (MQ=255)                                                                                                                                         cG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgtg            >  3:154746/1‑148 (MQ=255)                                                                                                                                             tACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACa        <  4:90428/149‑1 (MQ=255)                                                                                                                                                   gcgcTGCGGGTTGACGGCTACATTGCCCTCTCCTCGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGTACTGGTAAGCCAGCCCCTTCTGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCa  <  5:223285/146‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCA‑CCC‑‑‑TGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCA  >  NC_002947/5241577‑5241870

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