Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,199,077	C→T	6.1%	pseudogene (2604/2847 nt)	yehH →	DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,199,077	0	C	T	6.1%	207.4 / 2.7	65	pseudogene (2604/2847 nt)	yehH	DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

AGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACG  >  NC_000913/2198930‑2199213                                                                                                                                                    |                                                                                                                                         tgtggCTCTGTATCCTGGATTATTGCATGTGAAAGACGTCTGTTAAGCCGACTCTCTGGCGGGATTTGCCTGGGATCTGTGGACCGGCTGGCAGACCGCTGGCGCGCCGTCTAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCg                                                                                                                                         <  2:115524/146‑1 (MQ=255)  gagGCTCTGTATACGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGGCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCgg                                                                                                                                        <  6:300253/149‑1 (MQ=255)      ctctGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAAc                                                                                                                                    <  7:12396/149‑1 (MQ=255)           tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGa                                                                                                                                  <  8:213501/146‑1 (MQ=255)           tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATg                                                                                                                                <  5:73765/148‑1 (MQ=255)            aTCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGac                                                                                                                              >  7:260858/1‑149 (MQ=255)              ccGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGacac                                                                                                                            >  8:63774/1‑149 (MQ=255)                gggAGTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTACCCTCGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCCGTGCGTTAGGCGTTCTCGGTAACGATGACAccg                                                                                                                          >  4:112369/1‑149 (MQ=255)                 ggATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgc                                                                                                                         >  7:328403/1‑149 (MQ=255)                        tGCAGGTGAAAGCCGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCaaa                                                                                                                  >  3:261442/1‑149 (MQ=255)                                      gTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAAt                                                                                                    <  2:315371/149‑1 (MQ=255)                                              cgccgACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcct                                                                                            >  5:106491/1‑149 (MQ=255)                                               gccgACTCACTGGCGGGTTTTGCCTGGGATCTGTTTACCGCCTGGCAGACAGCTGGCGCGCCGTCCAAAGAGAGTAGGGCGTTCACTGCGTTAGGCGTTCTCGGTATCGACGCCCCCGCCCGCACACTGACGCCATTAATACGCgcctg                                                                                           >  2:75933/1‑149 (MQ=255)                                               gccgACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctg                                                                                           <  8:328403/149‑1 (MQ=255)                                                       aCTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    <  6:137633/148‑1 (MQ=255)                                                       aCTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    <  6:137640/148‑1 (MQ=255)                                                       aCTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    <  7:63774/148‑1 (MQ=255)                                                       aCTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    >  7:161419/1‑148 (MQ=255)                                                       aCTGGCGGGATTTGCCTGCGATCTGTTTACCGCCTCGCACACCGCTGGCGCGCCGTCGAAACAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    >  2:384290/1‑148 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTTGGTAACGAt                                                                                                                                 <  5:366548/99‑1 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTTGGTAACGAt                                                                                                                                 <  5:372268/99‑1 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTTGGTAACGAt                                                                                                                                 >  6:366548/1‑99 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTTGGTAACGAt                                                                                                                                 >  6:372268/1‑99 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATc                                                                               >  5:325940/1‑149 (MQ=255)                                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCAGTCGAAAGAGAGTTGTTCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATc                                                                               >  8:58490/1‑149 (MQ=255)                                                                      cTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAc                                                                    >  6:96457/1‑149 (MQ=255)                                                                       tGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAcg                                                                   >  1:76991/1‑149 (MQ=255)                                                                         ggATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAcgct                                                                 >  2:37480/1‑148 (MQ=255)                                                                          gATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACgcgc                                                                <  8:318645/149‑1 (MQ=255)                                                                           aTCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAcgcg                                                                 <  2:6587/147‑1 (MQ=255)                                                                           aTCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAcgcg                                                                 >  1:6587/1‑147 (MQ=255)                                                                               gTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCCTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGCCGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCg                                                           >  1:109634/1‑149 (MQ=255)                                                                                   aCCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgg                                                       >  8:137501/1‑149 (MQ=255)                                                                                      gCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgggtt                                                    >  6:305643/1‑149 (MQ=255)                                                                                      gCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgggtt                                                    >  1:11032/1‑149 (MQ=255)                                                                                           gCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcc                                                                                        >  2:34787/1‑108 (MQ=255)                                                                                           gCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcc                                                                                        <  1:34787/108‑1 (MQ=255)                                                                                           gCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGa                                                 <  1:221107/147‑1 (MQ=255)                                                                                                ccGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATtctc                                          <  3:30427/149‑1 (MQ=255)                                                                                                 cGCTGGCGCGCCGGCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCg                                         <  4:248619/149‑1 (MQ=255)                                                                                                  gCTGGCGCGCCGTCGACAGAGAGTTGGGCGTTCACTGCGATAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCgc                                        <  3:112369/149‑1 (MQ=255)                                                                                                       cgcgcCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    <  4:90486/100‑1 (MQ=255)                                                                                                       cgcgcCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt                                                                                    >  3:90486/1‑100 (MQ=255)                                                                                                          gcCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCg                                <  2:11032/149‑1 (MQ=255)                                                                                                           ccGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCgg                               >  5:128099/1‑149 (MQ=255)                                                                                                             gTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTa                             <  8:260858/149‑1 (MQ=255)                                                                                                                gAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTg                          <  5:96457/149‑1 (MQ=255)                                                                                                                 aaaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGCCACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGa                         >  8:341957/1‑149 (MQ=255)                                                                                                                  aaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                                                                                                   <  2:66381/74‑1 (MQ=255)                                                                                                                  aaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                                                                                                   >  1:66381/1‑74 (MQ=255)                                                                                                                  aaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGat                        >  2:303571/1‑149 (MQ=255)                                                                                                                    gagagTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGatat                      <  6:325940/149‑1 (MQ=255)                                                                                                                         ttGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGccc                 <  8:58159/149‑1 (MQ=255)                                                                                                                             gCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGCCGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTAt             >  5:105645/1‑149 (MQ=255)                                                                                                                             gCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTAt             >  5:219239/1‑149 (MQ=255)                                                                                                                              cGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATg            <  1:303571/149‑1 (MQ=255)                                                                                                                                ttCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCa          <  5:305643/149‑1 (MQ=255)                                                                                                                                 tCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAg         <  7:137501/149‑1 (MQ=255)                                                                                                                                   aCTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCt       >  3:327369/1‑149 (MQ=255)                                                                                                                                   aCTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCt       >  4:77978/1‑149 (MQ=255)                                                                                                                                     tGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCtt      >  1:97011/1‑148 (MQ=255)                                                                                                                                       cGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAAc   <  8:161419/149‑1 (MQ=255)                                                                                                                                       cGTTAGGCGTGCGCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAAc   <  1:316257/149‑1 (MQ=255)                                                                                                                                        gTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACg  <  2:76991/149‑1 (MQ=255)                                                                                                                                        gTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACg  <  1:384290/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                         AGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACG  >  NC_000913/2198930‑2199213

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