Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,719	Δ1 bp	7.7%	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	.	7.7%	210.6 / 10.5	65	coding (153/633 nt)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base G (25/35);  new base . (2/3);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCT  >  NC_002947/5241572‑5241867                                                                                                                                                     |                                                                                                                                                     aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                                     >  5:351103/1‑149 (MQ=255) aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                                     >  6:151352/1‑149 (MQ=255) aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                                     <  3:193174/149‑1 (MQ=255) aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCCCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                                     >  6:361577/1‑149 (MQ=255) aaTGATGGCCGCCACCCCCCGAACCCTGACCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                                     <  3:4680/149‑1 (MQ=255)          cgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                            <  3:347622/149‑1 (MQ=255)          cgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCg                                                                                                                                            <  2:142683/149‑1 (MQ=255)            ctcccccccGACCCCTGGCCCGTGCCGCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACGGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCggg                                                                                                                                          <  1:318519/146‑1 (MQ=255)            ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgca                                                                                                                                          <  5:362709/149‑3 (MQ=255)             cACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGAGCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGCgcac                                                                                                                                         <  5:66586/149‑4 (MQ=255)              aCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGtt                                                                                                                                        <  7:27638/149‑1 (MQ=255)                     gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGtt                                                                                                                                        <  4:255645/142‑1 (MQ=255)                     gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGtt                                                                                                                                        >  3:255645/1‑142 (MQ=255)                       aCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTAc                                                                                                                               <  5:129273/149‑1 (MQ=255)                         ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                              >  2:360251/1‑148 (MQ=255)                         ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                              >  6:349313/1‑148 (MQ=255)                             gCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAg                                                                                                                           <  5:30880/147‑1 (MQ=255)                                   gCCTCGGCGGTGGCCGGCAGCAACGCAGGCAGTGGATAAACCTCTTCCAGGTACTCGATGCTCGCCCGCGCCTGCACCGGCCCCTCACCCCCGTCGGTGCGCCAGGCCG‑GTACGCGGCCCGGCGGGGTGACGGCGCCAACGCCCTCGcc                                                                                                                   >  8:150721/1‑149 (MQ=255)                                     cTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCAc                                                                                                                 >  1:113944/1‑149 (MQ=255)                                      tCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                                >  3:180591/1‑133 (MQ=255)                                      tCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                                <  4:180591/133‑1 (MQ=255)                                                      gcaACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTgcagc                                                                                                <  5:349313/149‑1 (MQ=255)                                                         aCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                             >  2:351867/1‑149 (MQ=255)                                                           gCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACg                                                                                                                <  7:318904/128‑1 (MQ=255)                                                           gCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACg                                                                                                                >  8:318904/1‑128 (MQ=255)                                                               gCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTg                                                                                                    >  6:108823/1‑136 (MQ=255)                                                               gCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTg                                                                                                    <  5:108823/136‑1 (MQ=255)                                                                 tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGa                                                                                                                                      >  2:731/1‑100 (MQ=255)                                                                 tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGa                                                                                                                                      <  1:731/100‑1 (MQ=255)                                                                   tgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGc                                                                                    <  1:351867/148‑1 (MQ=255)                                                                     gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACgg                                                                                                                                   >  7:334072/1‑99 (MQ=255)                                                                     gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACgg                                                                                                                                   <  8:334072/99‑1 (MQ=255)                                                                     gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  6:351103/148‑1 (MQ=255)                                                                     gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  5:68631/148‑1 (MQ=255)                                                                        aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  4:289171/145‑1 (MQ=255)                                                                        aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  >  3:289171/1‑145 (MQ=255)                                                                          aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                             >  6:367877/1‑132 (MQ=255)                                                                          aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                             <  5:367877/132‑1 (MQ=255)                                                                          aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                            <  4:145580/149‑1 (MQ=255)                                                                          aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                            <  5:62503/149‑1 (MQ=255)                                                                               ttCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAag                                                                       >  4:196318/1‑149 (MQ=255)                                                                                 ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                        <  6:278199/146‑1 (MQ=255)                                                                                 ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCCCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                        <  5:361577/146‑1 (MQ=255)                                                                                    ggTACTCGATGAGCGCCGGCGACTGCACCCGCACCTCACCACCGGCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGGTGACGGCTACATAGCCCTCGCCCCGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagccagc                                                                  <  7:150721/149‑1 (MQ=255)                                                                                     gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCgg                                                                                     <  5:22216/129‑1 (MQ=255)                                                                                     gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCgg                                                                                     >  6:22216/1‑129 (MQ=255)                                                                                        cTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCt                                                              >  1:136787/1‑149 (MQ=255)                                                                                         tCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                             >  1:161728/1‑149 (MQ=255)                                                                                                        gACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAAt                                              <  2:113944/149‑1 (MQ=255)                                                                                                         aCTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATg                                             >  2:221079/1‑149 (MQ=255)                                                                                                          cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                      <  4:232429/123‑1 (MQ=255)                                                                                                          cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                      >  3:232429/1‑123 (MQ=255)                                                                                                                    aCCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                                  <  2:136787/149‑1 (MQ=255)                                                                                                                    aCCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                                  <  3:196318/149‑1 (MQ=255)                                                                                                                       tCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGa                               >  3:34108/1‑149 (MQ=255)                                                                                                                         accaccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGAc                              <  2:161728/148‑1 (MQ=255)                                                                                                                           caccGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAg                           >  1:173195/1‑149 (MQ=255)                                                                                                                                  ggTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGaa                     <  5:151352/148‑1 (MQ=255)                                                                                                                                    tGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                                 >  3:107540/1‑134 (MQ=255)                                                                                                                                    tGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                                 <  4:107540/134‑1 (MQ=255)                                                                                                                                     ccgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg                  >  1:310489/2‑148 (MQ=255)                                                                                                                                      cgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTg                        >  2:249721/1‑141 (MQ=255)                                                                                                                                      cgcAGGGCCGCGTAC‑CGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTg                        <  1:249721/141‑1 (MQ=255)                                                                                                                                         aaGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa               >  7:225361/1‑147 (MQ=255)                                                                                                                                          aGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg            <  4:34108/149‑1 (MQ=255)                                                                                                                                          aGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg            <  1:360251/149‑1 (MQ=255)                                                                                                                                           ggCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg           >  5:37346/1‑149 (MQ=255)                                                                                                                                           ggCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg           >  8:154742/1‑149 (MQ=255)                                                                                                                                            gCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt          <  6:37346/149‑1 (MQ=255)                                                                                                                                            gCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt          <  2:173195/149‑1 (MQ=255)                                                                                                                                                     gcgGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCt  >  5:303677/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCT  >  NC_002947/5241572‑5241867

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