Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,709	A→G	8.3%	L55L (TTG→CTG)	hmgC ←	maleylacetoacetate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,709	0	A	G	8.3%	203.9 / 7.3	72	L55L (TTG→CTG)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base G (2/4);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATA  >  NC_002947/5241566‑5241856                                                                                                                                                  |                                                                                                                                                     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)          gatgGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGc                                                                                                                                                   <  4:263056/138‑1 (MQ=255)          gatgGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGc                                                                                                                                                   >  3:263056/1‑138 (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                                                                                                                                  >  4:130997/1‑149 (MQ=255)               ccgccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGc                                                                                                                                  <  6:291407/149‑1 (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                                                                                                                               <  1:150709/149‑3 (MQ=255)                  ccACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                               <  2:173001/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:327970/148‑1 (MQ=255)                                                      gCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCAc                                                                                             <  5:77258/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                                                                   <  3:29955/148‑1 (MQ=255)                                                                                aaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGAc                                                                   <  2:216389/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                                                           >  3:177979/1‑148 (MQ=255)                                                                                        ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                           >  6:343167/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    >  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)                                                                                                                                                  |                                                                                                                                                     GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCC‑TCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCAC‑CGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATA  >  NC_002947/5241566‑5241856

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