Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,725	C→T	7.5%	Q49Q (CAG→CAA)	hmgC ←	maleylacetoacetate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,725	0	C	T	7.5%	195.0 / 3.9	66	Q49Q (CAG→CAA)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base T (2/3);  total (28/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.78e-01

CGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCT  >  NC_002947/5241581‑5241867                                                                                                                                                  |                                                                                                                                               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                                                                                                                                        >  3:255645/1‑142 (MQ=255)            gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGtt                                                                                                                                        <  4:255645/142‑1 (MQ=255)              aCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTAc                                                                                                                               <  5:129273/149‑1 (MQ=255)                ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                              >  6:349313/1‑148 (MQ=255)                ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                              >  2:360251/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                                                                                                                >  8:318904/1‑128 (MQ=255)                                                  gCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACg                                                                                                                <  7:318904/128‑1 (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                                                                                                                                   <  8:334072/99‑1 (MQ=255)                                                            gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACgg                                                                                                                                   >  7:334072/1‑99 (MQ=255)                                                            gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  5:68631/148‑1 (MQ=255)                                                            gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  6:351103/148‑1 (MQ=255)                                                               aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  >  3:289171/1‑145 (MQ=255)                                                               aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                                  <  4:289171/145‑1 (MQ=255)                                                                 aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                             <  5:367877/132‑1 (MQ=255)                                                                 aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                             >  6:367877/1‑132 (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                                                                      >  3:232429/1‑123 (MQ=255)                                                                                                 cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                                      <  4:232429/123‑1 (MQ=255)                                                                                                           aCCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                                  <  3:196318/149‑1 (MQ=255)                                                                                                           aCCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                                  <  2:136787/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                                 <  4:107540/134‑1 (MQ=255)                                                                                                                           tGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                                 >  3:107540/1‑134 (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            <  1:360251/149‑1 (MQ=255)                                                                                                                                 aGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg            <  4:34108/149‑1 (MQ=255)                                                                                                                                  ggCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg           >  5:37346/1‑149 (MQ=255)                                                                                                                                  ggCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg           >  8:154742/1‑149 (MQ=255)                                                                                                                                   gCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt          <  2:173195/149‑1 (MQ=255)                                                                                                                                   gCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt          <  6:37346/149‑1 (MQ=255)                                                                                                                                            gcgGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCt  >  5:303677/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               CGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGTGTACAGCT  >  NC_002947/5241581‑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: ‑