Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,709	A→G	7.6%	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	7.6%	186.2 / 5.1	66	L55L (TTG→CTG)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base G (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 = 8.18e-01

ATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAG  >  NC_002947/5241561‑5241857                                                                                                                                                     |                                                                                                                                                      aTGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCaa                                                                                                                                                       <  3:303544/149‑1 (MQ=255)    tCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAgg                                                                                                                                                     <  7:247926/148‑1 (MQ=255)      gCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGGACCAGCACCTCACCACCGTCGGTGCGCAAGGCCg                                                                                                                                                  <  2:267749/149‑1 (MQ=255)      gCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCg                                                                                                                                                  >  1:142683/1‑149 (MQ=255)            aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                           >  5:351103/1‑149 (MQ=255)            aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                           <  3:193174/149‑1 (MQ=255)            aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTAcgc                                                                                                                                           >  6:151352/1‑149 (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                                                                                                                    >  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                                                                                                                      <  4:180591/133‑1 (MQ=255)                                                 tCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                      >  3:180591/1‑133 (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                                                                                                                            <  1:731/100‑1 (MQ=255)                                                                            tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGa                                                                                                                            >  2:731/1‑100 (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                                                                        <  5:68631/148‑1 (MQ=255)                                                                                gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                        <  6:351103/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                                                                  <  5:62503/149‑1 (MQ=255)                                                                                     aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                  <  4:145580/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                       >  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  <  1:360251/149‑1 (MQ=255)                                                                                                                                                     aGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg  <  4:34108/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAG  >  NC_002947/5241561‑5241857

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