Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,891,873	G→A	7.8%	V345V (GTG→GTA)	ttuD →	hydroxypyruvate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,891,873	0	G	A	7.8%	184.9 / 4.3	64	V345V (GTG→GTA)	ttuD	hydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTA  >  NC_002947/4891727‑4892016                                                                                                                                                    |                                                                                                                                                 gTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCgtgt                                                                                                                                                 <  2:307579/148‑1 (MQ=255)     aaCACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCAGGGCGTACACGc                                                                                                                                             >  2:69316/1‑149 (MQ=255)       cacGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCt                                                                                                                                          >  3:400923/1‑149 (MQ=255)        acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTg                                                                                                                                         <  2:230192/149‑1 (MQ=255)         cGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTTACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTgg                                                                                                                                        <  8:244656/149‑1 (MQ=255)         cGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTgg                                                                                                                                        <  5:48172/149‑1 (MQ=255)            ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGccc                                                                                                                                     >  1:325532/1‑148 (MQ=255)            ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCg                                                                                                                                     <  1:103545/149‑1 (MQ=255)            ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCg                                                                                                                                     >  4:180593/1‑149 (MQ=255)             cAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCAGGGCGTACACGcccggcagg                                                                                                                                     <  1:69316/149‑7 (MQ=255)                ccGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑a                                                                                                                                <  6:359904/149‑1 (MQ=255)                ccGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑a                                                                                                                                >  8:326506/1‑149 (MQ=255)                  gCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑aca                                                                                                                              <  4:329050/149‑1 (MQ=255)                     cAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCAGGGCGTACACGcccggca                                                                                                                                       <  3:145024/139‑5 (MQ=255)                        aGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACg                                                                                                                        <  4:187918/149‑1 (MQ=255)                              cTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATc                                                                                                                   >  4:353652/1‑148 (MQ=255)                               tGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGa                                                                                                                 >  6:56338/1‑149 (MQ=255)                               tGCGGGATCCTGTCGGGTTGCGAAACCACGGTGACCGGGCGGGGCAATTGCCGTGGCGGGCGTATCGCCGAGTGCCTGCTCAGCCTCACCGATAGCCTGTATGGCCTGCC‑GTGCGTGTACGCCCTGGCCGGTG‑TCACCGACGGCATCGa                                                                                                                 <  6:938/149‑1 (MQ=255)                                 cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGAt                                                                                                                <  7:374417/148‑1 (MQ=255)                                 cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATg                                                                                                               >  5:148832/1‑149 (MQ=255)                                      tCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCg                                                                                                          <  6:221954/149‑1 (MQ=255)                                       ccTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCg                                                                                                          >  8:189629/1‑148 (MQ=255)                                           tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGGCGGTGCACACCGACGGCATCGATGGCTCGgag                                                                                                       <  1:4671/148‑2 (MQ=255)                                           tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGgaa                                                                                                       <  3:279963/147‑1 (MQ=255)                                           tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGgaa                                                                                                       <  1:6463/147‑1 (MQ=255)                                             gggTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGt                                                                                                                                   >  8:344057/1‑118 (MQ=255)                                             gggTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGt                                                                                                                                   <  7:344057/118‑1 (MQ=255)                                                    gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGAc                                                                                                                         >  2:99307/1‑120 (MQ=255)                                                     aaaCCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                           <  3:353652/149‑1 (MQ=255)                                                      aaCCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                           <  4:125976/148‑1 (MQ=255)                                                      aaCCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                           >  3:125976/1‑148 (MQ=255)                                                      aaCCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTg                                                                                          <  7:189629/149‑1 (MQ=255)                                                            cgtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCa                                                                                    >  7:254789/1‑149 (MQ=255)                                                                  cgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAccc                                                                              >  1:165699/1‑149 (MQ=255)                                                                      cgcgGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACTGACGGCATCGATGGCTCGGAAGAAAACGCGGGTGCCTTCATGACCCCgg                                                                           <  2:325532/148‑1 (MQ=255)                                                                                gCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                 >  4:5092/1‑148 (MQ=255)                                                                                        gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCg                                                         >  7:174743/1‑148 (MQ=255)                                                                                        gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCg                                                         >  7:180050/1‑148 (MQ=255)                                                                                        gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGa                                                        <  3:180593/149‑1 (MQ=255)                                                                                         ggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGaa                                                       <  2:112378/149‑1 (MQ=255)                                                                                            gTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGcc                                                    <  7:326506/149‑1 (MQ=255)                                                                                                cGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTgg                                                <  8:254789/149‑1 (MQ=255)                                                                                                 gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                               <  5:56338/149‑1 (MQ=255)                                                                                                 gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                               >  3:201962/1‑149 (MQ=255)                                                                                                 gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                               <  4:400923/149‑1 (MQ=255)                                                                                                  ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                 >  2:44812/1‑130 (MQ=255)                                                                                                  ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                 <  1:44812/130‑1 (MQ=255)                                                                                                  ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                              <  6:148832/149‑1 (MQ=255)                                                                                                     aGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTg                                           >  3:117193/1‑149 (MQ=255)                                                                                                             cTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAg                                   <  4:201962/149‑1 (MQ=255)                                                                                                                    tCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGcc                                                                         >  1:11684/1‑104 (MQ=255)                                                                                                                    tCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGcc                                                                         <  2:11684/104‑1 (MQ=255)                                                                                                                     cACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGaa                                                       <  2:278865/121‑1 (MQ=255)                                                                                                                     cACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGaa                                                       >  1:278865/1‑121 (MQ=255)                                                                                                                       ccGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                         >  6:280643/1‑149 (MQ=255)                                                                                                                        cGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGGCGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                         >  6:226249/1‑148 (MQ=255)                                                                                                                          aaaGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGa                       >  6:404480/1‑148 (MQ=255)                                                                                                                                  aaaGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGAcaacaacaa              <  8:180050/149‑1 (MQ=255)                                                                                                                                   aaGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGAcaacaacaa              <  8:174743/148‑1 (MQ=255)                                                                                                                                   aaGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACaacaacaac             <  5:280643/149‑1 (MQ=255)                                                                                                                                        cggccagggCGTACACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                         >  2:139238/7‑133 (MQ=255)                                                                                                                                        cggccagggCGTACACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                         <  1:139238/127‑1 (MQ=255)                                                                                                                                           cc‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggctacgg     <  5:404480/149‑1 (MQ=255)                                                                                                                                            c‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggctacggc    >  5:348801/1‑149 (MQ=255)                                                                                                                                               ggCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAAcggctacggcta  >  8:271247/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCC‑GGGCGTGTACGCCCTGGCCGGTG‑ACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTA  >  NC_002947/4891727‑4892016

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