Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,725	C→T	8.5%	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	8.5%	197.6 / 7.0	70	Q49Q (CAG→CAA)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/4);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCA  >  NC_002947/5241580‑5241870                                                                                                                                                    |                                                                                                                                                  ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                                >  1:173001/1‑149 (MQ=255) ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                                <  6:47825/149‑1 (MQ=255) ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                                >  4:130997/1‑149 (MQ=255) ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                                <  6:291407/149‑1 (MQ=255) ccgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                                <  3:271638/149‑1 (MQ=255)  cgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                               >  6:187190/1‑149 (MQ=255)   gccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCgg                                                                                                                                              >  1:366546/1‑149 (MQ=255)    ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                             <  2:173001/149‑3 (MQ=255)    ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                             <  1:150709/149‑3 (MQ=255)     cACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGt                                                                                                                                            <  2:268718/149‑1 (MQ=255)     cACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGt                                                                                                                                            >  2:93934/1‑149 (MQ=255)      aCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgcacc                                                                                                                                           <  3:34764/149‑5 (MQ=255)       ccccccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTg                                                                                                                                          <  6:88083/149‑1 (MQ=255)           ccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACgg                                                                                                                                      <  5:187190/149‑1 (MQ=255)           ccGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCAGCGACTGCCCCAGCACCTCACCAC‑CGTCGGTGCGCAAGGACG‑GTACGCCGCCCTGCGGGTTGACgg                                                                                                                                      >  8:142485/1‑149 (MQ=255)             gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCt                                                                                                                                    <  3:130997/149‑1 (MQ=255)                   tGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAg                                                                                                                              >  6:77258/1‑149 (MQ=255)                        cGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTc                                                                                                                         >  2:232708/1‑149 (MQ=255)                             cTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCAc                                                                                                                    >  2:287436/1‑149 (MQ=255)                               cggcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAG‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgc                                                                                                                  >  2:90122/1‑149 (MQ=255)                                ggcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgct                                                                                                                 <  2:366546/149‑1 (MQ=255)                                 gcggTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACgctg                                                                                                                <  1:248202/149‑1 (MQ=255)                                    gTGGCCGGCAGCAACGCGGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTc                                                                                                             <  5:68110/149‑1 (MQ=255)                                     tGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCa                                                                                                            >  1:216389/1‑149 (MQ=255)                                       gCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCAc                                                                                                           <  5:77258/148‑1 (MQ=255)                                       gCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCAc                                                                                                           <  5:327970/148‑1 (MQ=255)                                        ccGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                               <  5:329083/111‑1 (MQ=255)                                        ccGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                               >  6:329083/1‑111 (MQ=255)                                        ccGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCAccc                                                                                                         >  5:54161/1‑149 (MQ=255)                                            cagcaACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTgca                                                                                                     <  7:155087/149‑1 (MQ=255)                                                aaCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGcagcag                                                                                                 >  5:156270/1‑149 (MQ=255)                                                aaCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGcagcag                                                                                                 >  4:63465/1‑149 (MQ=255)                                                     aGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGa                                                                                            >  7:344352/1‑149 (MQ=255)                                                      ggCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGAc                                                                                           <  2:55473/149‑1 (MQ=255)                                                         tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCgg                                                                                        >  4:119089/1‑149 (MQ=255)                                                                 aaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGAc                                                                                 <  2:216389/148‑1 (MQ=255)                                                                 aaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGAc                                                                                 <  3:29955/148‑1 (MQ=255)                                                                  aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                               <  1:93934/149‑1 (MQ=255)                                                                  aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                               >  7:336332/1‑149 (MQ=255)                                                                   ccTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTgg                                                                              <  1:232708/149‑1 (MQ=255)                                                                         ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                         >  6:343167/1‑148 (MQ=255)                                                                         ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                         >  3:177979/1‑148 (MQ=255)                                                                         ccAGGTACTCGATGATCGCCAGCGACTGCACCAGCACCTCACCAC‑CGTCGTTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                           <  7:142485/146‑1 (MQ=255)                                                                               aCTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGcccc                                                                  >  3:28406/1‑149 (MQ=255)                                                                                  cGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGCCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGAACGGCAGGGACTGGTAAGCCAGCCCCtt                                                                >  7:100638/1‑148 (MQ=255)                                                                                   gatgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                                >  7:300099/1‑147 (MQ=255)                                                                                   gatgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTc                                                               <  5:219053/148‑1 (MQ=255)                                                                                   gatgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCa                                                              <  3:63465/149‑1 (MQ=255)                                                                                   gatgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGCTAGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                                >  8:109615/1‑147 (MQ=255)                                                                                    atgatCGCCGGCGACTGCACCAGCACCTCACCACTCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                                >  1:50507/1‑147 (MQ=255)                                                                                     tgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGt                                                            <  3:96087/149‑1 (MQ=255)                                                                                     tgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGt                                                            <  8:336332/149‑1 (MQ=255)                                                                                        tCGCCGGCGACTGCACCAACACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCa                                                              >  4:303494/1‑144 (MQ=255)                                                                                        tCGCCGGCGACTGCACCAACACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCa                                                              <  3:303494/144‑1 (MQ=255)                                                                                            cGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAggg                                                     <  5:343167/149‑1 (MQ=255)                                                                                                       ccAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCAcc                                          <  6:54161/149‑1 (MQ=255)                                                                                                          gCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCgg                                       <  1:287436/149‑1 (MQ=255)                                                                                                          gCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCGTCAGTGCCAGGGCAATGCGCACCCgg                                       <  1:201730/149‑1 (MQ=255)                                                                                                           cACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGt                                      <  3:119089/149‑1 (MQ=255)                                                                                                                 accac‑cGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGa                                  <  8:100638/147‑1 (MQ=255)                                                                                                                 accac‑cGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGAc                                 >  3:96718/1‑148 (MQ=255)                                                                                                                   cac‑cGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAg                              <  6:62417/149‑1 (MQ=255)                                                                                                                             cccgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                     >  4:36836/3‑149 (MQ=255)                                                                                                                                gcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa                  >  5:174231/1‑149 (MQ=255)                                                                                                                                gcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa                  >  3:69600/1‑149 (MQ=255)                                                                                                                                  aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa                  >  5:82900/1‑147 (MQ=255)                                                                                                                                  aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATa                <  1:90122/149‑1 (MQ=255)                                                                                                                                    ggCCT‑GTACGCGGCCCTGCGGGTGGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg              <  4:243047/149‑1 (MQ=255)                                                                                                                                      ccG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgtg            <  8:344352/149‑1 (MQ=255)                                                                                                                                                 gCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCa  <  4:96718/149‑1 (MQ=255)                                                                                                                                                   ccTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAgg                             >  4:176055/1‑120 (MQ=255)                                                                                                                                                   ccTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAgg                             <  3:176055/120‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCTCCA  >  NC_002947/5241580‑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: ‑