Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,719	Δ1 bp	8.6%	coding (153/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,719	0	G	.	8.6%	222.0 / 14.2	70	coding (153/633 nt)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base . (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 = 5.33e-01

AATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTG  >  NC_002947/5241572‑5241860                                                                                                                                                       |                                                                                                                                              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                                                                                                                                      <  3:271638/149‑1 (MQ=255)         ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                      <  6:291407/149‑1 (MQ=255)         ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                      >  4:130997/1‑149 (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                                                                                       >  4:63465/1‑149 (MQ=255)                                                         aaCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGcagcag                                                                                       >  5:156270/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                                                                       <  3:29955/148‑1 (MQ=255)                                                                          aaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGAc                                                                       <  2:216389/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)                                                                                                                                                       |                                                                                                                                              AATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTG  >  NC_002947/5241572‑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: ‑