Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,714:1	+C	8.3%	coding (158/633 nt)	hmgC ←	maleylacetoacetate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,714	1	.	C	8.3%	276.7 / 17.0	72	coding (158/633 nt)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base . (31/35);  new base C (2/4);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTG  >  NC_002947/5241566‑5241860                                                                                                                                                        |                                                                                                                                                   gCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCg                                                                                                                                                     >  4:96087/1‑149 (MQ=255)   aaCCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCg‑g                                                                                                                                                   >  4:29955/1‑148 (MQ=255)       aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                              >  5:62417/1‑149 (MQ=255)        atgacgGCCGCCACCCCCCCAACCCTGG‑CCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCgg                                                                                                                                            <  4:224567/149‑1 (MQ=255)         tgatgGCCGCCACCCCCCGAACCCTGGCCCGTGCCGTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACAAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTAcgcg                                                                                                                                             >  1:55473/1‑149 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                      >  1:173001/1‑149 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc                                                                                                                                      <  6:47825/149‑1 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                      <  6:291407/149‑1 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                      >  4:130997/1‑149 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                      <  3:271638/149‑1 (MQ=255)                cgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                     >  6:187190/1‑149 (MQ=255)                 gccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCgg                                                                                                                                    >  1:366546/1‑149 (MQ=255)                  ccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                   <  2:173001/149‑3 (MQ=255)                  ccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                   <  1:150709/149‑3 (MQ=255)                   cACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGt                                                                                                                                  <  2:268718/149‑1 (MQ=255)                   cACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGt                                                                                                                                  >  2:93934/1‑149 (MQ=255)                    aCCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgcacc                                                                                                                                 <  3:34764/149‑5 (MQ=255)                     ccccccGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTg                                                                                                                                <  6:88083/149‑1 (MQ=255)                         ccGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACgg                                                                                                                            <  5:187190/149‑1 (MQ=255)                         ccGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCAGCGACTGCCCCAGCACCTCACCAC‑CGTCGGTGCGCAAGGACG‑GTACGCCGCCCTGCGGGTTGACgg                                                                                                                            >  8:142485/1‑149 (MQ=255)                           gAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCt                                                                                                                          <  3:130997/149‑1 (MQ=255)                                 tGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAg                                                                                                                    >  6:77258/1‑149 (MQ=255)                                      cGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTc                                                                                                               >  2:232708/1‑149 (MQ=255)                                           c‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑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                                                                                                                                     >  6:329083/1‑111 (MQ=255)                                                       ccGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                     <  5:329083/111‑1 (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                                                                     >  7:336332/1‑149 (MQ=255)                                                                                 aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                     <  1:93934/149‑1 (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                                                  <  8:336332/149‑1 (MQ=255)                                                                                                    tgatCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGt                                                  <  3:96087/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        >  3:69600/1‑149 (MQ=255)                                                                                                                                               gcAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa        >  5:174231/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)                                                                                                                                                        |                                                                                                                                                   GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTG  >  NC_002947/5241566‑5241860

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