Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,609,406	T→C	5.9%	F307S (TTC→TCC)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,609,406	0	T	C	5.9%	212.6 / 2.9	68	F307S (TTC→TCC)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/35)
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.71e-01

TTGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCG  >  NC_000913/2609259‑2609552                                                                                                                                                    |                                                                                                                                                   ttGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTc                                                                                                                                                   >  3:111681/1‑149 (MQ=255)      aaTGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              >  1:359143/1‑147 (MQ=255)       aTGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTg                                                                                                                                             >  7:178580/1‑149 (MQ=255)           aCTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCgg                                                                                                                                          >  8:199680/1‑148 (MQ=255)               gggATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGaa                                                                                                                                       <  4:3322/147‑1 (MQ=255)                 gATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCCGCACGCGGAaaaatt                                                                                                                                   <  2:359143/149‑3 (MQ=255)                       ggTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGa                                                                                                                             <  7:181567/149‑1 (MQ=255)                        gTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGGCTAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAc                                                                                                                            <  6:302024/149‑1 (MQ=255)                        gTGCCGAAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAc                                                                                                                            <  8:170117/149‑1 (MQ=255)                              cAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGc                                                                                                                      >  5:239816/1‑149 (MQ=255)                                 cGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCgg                                                                                                                    >  8:354332/1‑148 (MQ=255)                                    gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCcgc                                                                                                                                                >  4:123411/1‑117 (MQ=255)                                    gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCcgc                                                                                                                                                <  3:123411/117‑1 (MQ=255)                                    gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAtt                                                                                                                >  6:345766/1‑149 (MQ=255)                                     cacaTACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTg                                                                                                               <  1:308525/149‑1 (MQ=255)                                         tACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAt                                                                                                                 <  7:10769/143‑1 (MQ=255)                                         tACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAt                                                                                                                 >  8:10769/1‑143 (MQ=255)                                              ttCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCATGCCATTTCGTTTGTTTTGATACAGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCgg                                                                                                       <  5:319853/148‑1 (MQ=255)                                                 gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcg                                                                                                   >  1:351898/1‑149 (MQ=255)                                                 gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcg                                                                                                   >  7:125326/1‑149 (MQ=255)                                                 gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTc                                                                                                    <  4:56933/148‑1 (MQ=255)                                                  cTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCTGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  <  7:199680/149‑1 (MQ=255)                                                     gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  >  7:173899/1‑146 (MQ=255)                                                     gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  <  8:173899/146‑1 (MQ=255)                                                       ttCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              <  7:296426/100‑3 (MQ=255)                                                       ttCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              >  8:296426/1‑98 (MQ=255)                                                          tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAAcc                                                                                          >  1:7393/1‑149 (MQ=255)                                                          tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACATTTGCCTCTCCTGCCATTTCGTTGGTTTTGATAACGGCTTTATGCAATTTTTCCGCGTGAGGGAAAAGTCGATGACGATGGCGGAATTGCTGATAGGGTAGCGTAAAAc                                                                                           >  4:203165/1‑148 (MQ=255)                                                            ctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTa                                                                                        <  4:111681/149‑1 (MQ=255)                                                             tgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTAc                                                                                       >  2:88006/1‑149 (MQ=255)                                                              tgAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCGTCTCCTGCCTTTTCGTTGGGTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACg                                                                                      <  6:49520/148‑1 (MQ=255)                                                                aaGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGt                                                                                    <  6:268071/149‑1 (MQ=255)                                                                  gTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCt                                                                                  <  7:350239/149‑1 (MQ=255)                                                                     gAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGa                                                                                <  8:178580/148‑1 (MQ=255)                                                                      aaCGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGaa                                                                               >  8:2035/1‑148 (MQ=255)                                                                             acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAGGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtg                                                                       >  2:155951/1‑149 (MQ=255)                                                                             acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtg                                                                       >  2:40524/1‑149 (MQ=255)                                                                                caGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtggtg                                                                    <  7:354332/149‑1 (MQ=255)                                                                                 aGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCgg                                                                                                                                          >  1:283098/1‑78 (MQ=255)                                                                                 aGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTAATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATttgtgg                                                                   >  3:192630/1‑149 (MQ=255)                                                                                   tCATTTGATGAACCTTGGCCTCTCCTGCCATTTCGTAGGATTAGATACCCTCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGCTGGCCGAATAGCTGATCGGGTCCCGTAAAACCTACGGTCTGGATCGGATTGgtggtg                                                                 >  4:57079/1‑149 (MQ=255)                                                                                      tttGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                              <  1:226924/149‑1 (MQ=255)                                                                                          cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCTGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGTTGGTGTTcgcc                                                           <  6:239816/148‑1 (MQ=255)                                                                                          cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCt                                                                            >  6:321071/1‑131 (MQ=255)                                                                                          cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCt                                                                            <  5:321071/131‑1 (MQ=255)                                                                                          cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccg                                                          >  8:325776/1‑149 (MQ=255)                                                                                          cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccg                                                          >  7:142531/1‑149 (MQ=255)                                                                                           ttAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAGAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                         >  2:91096/1‑149 (MQ=255)                                                                                             aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                         <  2:351898/147‑1 (MQ=255)                                                                                                    gCCTCTCCTGCCATTTCGTTGGTTTTGCTACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAATTCGATGACGATGGCGTATTTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                <  2:71380/149‑1 (MQ=255)                                                                                                    gCCTCTCCTGCCATTTCGTTGGTGTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                <  8:390227/149‑1 (MQ=255)                                                                                                      ctctCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACTGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                              <  8:125326/149‑1 (MQ=255)                                                                                                          ccTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAagag                                          >  3:349552/1‑149 (MQ=255)                                                                                                            tGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGc                                         >  6:93791/1‑148 (MQ=255)                                                                                                              ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAAc                                      <  5:345766/149‑1 (MQ=255)                                                                                                               cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGAtt                                                                        >  6:287653/1‑114 (MQ=255)                                                                                                               cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGAtt                                                                        <  5:287653/114‑1 (MQ=255)                                                                                                               cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACg                                     >  3:168486/1‑149 (MQ=255)                                                                                                                   tCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCTGGTAGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                  <  3:203165/148‑1 (MQ=255)                                                                                                                   tCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                  <  1:68238/148‑1 (MQ=255)                                                                                                                    cGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTg                                 <  1:40524/148‑1 (MQ=255)                                                                                                                     gTTGGTTTTGATACCGGCTTTATGCAATTTATCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGc                                <  1:88006/148‑1 (MQ=255)                                                                                                                              gATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGa                       <  8:142531/148‑1 (MQ=255)                                                                                                                                    ggCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGt                <  7:325776/149‑1 (MQ=255)                                                                                                                                        ttATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGGCTGCAAACGCTGAAACTGGTgcgc            >  2:219376/1‑149 (MQ=255)                                                                                                                                            gCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgcgc            >  3:57137/1‑145 (MQ=255)                                                                                                                                             cAATTTTTCCGCGTGCGGGAAGAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGAt       <  1:91096/149‑1 (MQ=255)                                                                                                                                             cAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGAt       >  3:323751/1‑149 (MQ=255)                                                                                                                                                   ttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgcg  <  1:219376/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TTGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCG  >  NC_000913/2609259‑2609552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑