Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,339	T→C	6.3%	D178G (GAC→GGC)	hmgC ←	maleylacetoacetate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,339	0	T	C	6.3%	186.7 / 2.3	63	D178G (GAC→GGC)	hmgC	maleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base T (24/35);  new base C (2/2);  total (26/37)
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.66e-01

GAAGCTGTCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATC  >  NC_002947/5241191‑5241485                                                                                                                                                     |                                                                                                                                                   gAAGCTGTCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGt                                                                                                                                                    >  8:12326/1‑149 (MQ=255)       gTCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGt                                                                                                                                              >  7:124712/1‑149 (MQ=255)        tCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAgg                                                                                                                                               <  5:307314/147‑1 (MQ=255)        tCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGt                                                                                                                                              <  6:263324/148‑1 (MQ=255)           gAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGAt                                                                                                                                          <  7:307511/149‑1 (MQ=255)                 aTCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTg                                                                                                                                      >  8:126047/1‑147 (MQ=255)                 aTCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGaa                                                                                                                                    <  6:78254/149‑1 (MQ=255)                    ggCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCg                                                                                                                                 >  7:12829/1‑149 (MQ=255)                    ggCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCg                                                                                                                                 <  2:99666/149‑1 (MQ=255)                       aaGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTc                                                                                                                              <  3:34994/149‑1 (MQ=255)                              cggcggCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCcgc                                                                                                                        <  6:272689/148‑1 (MQ=255)                                 cggcgAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTa                                                                                                                     >  8:162191/1‑148 (MQ=255)                                  ggcgAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGaa                                                                                                                                    <  3:284880/132‑1 (MQ=255)                                  ggcgAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGaa                                                                                                                                    >  4:284880/1‑132 (MQ=255)                                        ccGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCg                                                                                                             >  3:7691/1‑149 (MQ=255)                                        ccGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCg                                                                                                             <  2:98541/149‑1 (MQ=255)                                           ttcacttcacTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGa                                                                                                          <  8:124712/149‑1 (MQ=255)                                              acttcactGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCa                                                                                                       <  7:265704/149‑1 (MQ=255)                                                  cactGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCcgc                                                                                                                        >  7:201406/1‑128 (MQ=255)                                                  cactGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCcgc                                                                                                                        <  8:201406/128‑1 (MQ=255)                                                   actGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGa                                                                                                          <  7:119867/141‑1 (MQ=255)                                                   actGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGa                                                                                                          >  8:119867/1‑141 (MQ=255)                                                   actGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGCCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTa                                                                                                   >  4:270634/1‑148 (MQ=255)                                                   actGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGCCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTa                                                                                                   >  4:265974/1‑148 (MQ=255)                                                     tGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCGGTGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCTTTCGGCCGCGTACAACTGCGGGATCAGGTac                                                                                                  <  5:352705/147‑1 (MQ=255)                                                         ccGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATcc                                                                                            <  7:162191/149‑1 (MQ=255)                                                           gggCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGc                                                                                          <  8:12829/149‑1 (MQ=255)                                                                    ggCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAAc                                                                                                                 >  4:65351/1‑117 (MQ=255)                                                                    ggCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAAc                                                                                                                 <  3:65351/117‑1 (MQ=255)                                                                     gctgctGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGc                                                                                <  5:225406/149‑1 (MQ=255)                                                                            ggCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGCCGTGTCGc                                                                         >  8:248349/1‑149 (MQ=255)                                                                                            gCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCg                                                           >  4:278808/1‑147 (MQ=255)                                                                                            gCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCg                                                           <  3:278808/147‑1 (MQ=255)                                                                                            gCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGCCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGt                                                          <  3:265974/148‑1 (MQ=255)                                                                                            gCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGCCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGt                                                          <  3:270634/148‑1 (MQ=255)                                                                                             cGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGa                                                        <  7:21079/149‑1 (MQ=255)                                                                                                 cGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGc                                                    <  3:203845/149‑1 (MQ=255)                                                                                                      ggggcgcGGCCGCCTGGGCGGCAACCCGCAGGATGCGCGGGTAGCTGTCGAGGTCGAGGTTGAAGCGGTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGGCGCCAAAGCAGAAACCGTGATCGCCAATc                                               <  6:294110/146‑1 (MQ=255)                                                                                                           gcgGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTa                                                                                                   <  8:211322/92‑1 (MQ=255)                                                                                                           gcgGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTg                                          >  8:293847/1‑149 (MQ=255)                                                                                                             ggCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCt                                        >  6:30307/1‑149 (MQ=255)                                                                                                              gccgccAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCGAAAGCAGAGACCGTGATCGCCAATCAGTTGCTc                                       >  7:321923/1‑149 (MQ=255)                                                                                                                   cAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTg                                  >  1:212785/1‑149 (MQ=255)                                                                                                                    aGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGc                                 >  5:216015/1‑149 (MQ=255)                                                                                                                     gTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGcc                                >  8:118180/1‑149 (MQ=255)                                                                                                                     gTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGcc                                >  8:113248/1‑149 (MQ=255)                                                                                                                         ggCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTg                                          <  7:192874/135‑1 (MQ=255)                                                                                                                         ggCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTg                                          >  8:192874/1‑135 (MQ=255)                                                                                                                            aaCCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGcc                         <  2:212785/149‑1 (MQ=255)                                                                                                                               ccGCAGGATGCGCGGGAAGCTGTCGTGGTCGATGTTGTAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCATCTGCTCCACTGCCGCGAGCCCtt                      <  7:51761/149‑1 (MQ=255)                                                                                                                                cGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTg                     >  3:220668/1‑149 (MQ=255)                                                                                                                                 gCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTgg                    <  5:30307/149‑1 (MQ=255)                                                                                                                                  cAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGc                   <  3:266793/149‑1 (MQ=255)                                                                                                                                   aGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTTTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCt                  <  5:203126/149‑1 (MQ=255)                                                                                                                                     gATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGa                >  8:7115/1‑149 (MQ=255)                                                                                                                                       tGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATc              <  7:7115/149‑1 (MQ=255)                                                                                                                                            gggAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTg         >  8:202477/1‑149 (MQ=255)                                                                                                                                              gAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGc       <  7:118180/149‑1 (MQ=255)                                                                                                                                              gAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGc       <  7:113248/149‑1 (MQ=255)                                                                                                                                                 gCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTgg                    >  6:286851/1‑133 (MQ=255)                                                                                                                                                 gCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTgg                    <  5:286851/133‑1 (MQ=255)                                                                                                                                                  cTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGAt   <  4:220668/149‑1 (MQ=255)                                                                                                                                                   tGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGa    >  4:81325/1‑147 (MQ=255)                                                                                                                                                   tGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGa    <  3:81325/147‑1 (MQ=255)                                                                                                                                                     tCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATc  <  7:248349/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GAAGCTGTCCGAGCGAATCGGCAAGCTTCCGGCGGCGAACCGTTCACTTCACTGAGCCGGGCTGTCCGGCTGCTGGGCGGGATGGGCCTTGGCGAACGCCGGGTGCGCGGCCGCCAGTGCGGCAACCCGCAGGATGCGCGGGAAGCTGTCGAGGTCGATGTTGAAGCGTTCGGCCGCGTACAACTGCGGGATCAGGTACACATCCGCCAGCCCAGGCGTGTCGCCAAAGCAGAAACCGTGATCGCCAATCAGTTGCTCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATC  >  NC_002947/5241191‑5241485

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