Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,609,899	C→T	5.3%	S471S (AGC→AGT)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,609,899	0	C	T	5.3%	244.8 / 2.5	75	S471S (AGC→AGT)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (2/2);  total (39/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 = 9.76e-01

ACGTTTGATTCGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTG  >  NC_000913/2609756‑2610045                                                                                                                                                |                                                                                                                                                   aCGTTTGATTCGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCTCTGGAGCATGc                                                                                                                                                >  1:204545/1‑148 (MQ=255)   gTTTGATTCGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCTCTGGAGCATGCTc                                                                                                                                              <  2:204545/148‑1 (MQ=255)        aTTCGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGa                                                                                                                                         >  6:128000/1‑148 (MQ=255)           cGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCGTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATa                                                                                                                                       <  4:143658/147‑1 (MQ=255)                    tGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAAtt                                                                                                                             >  3:209385/1‑148 (MQ=255)                     gCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTg                                                                                                                            >  2:211371/1‑148 (MQ=255)                      cTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATaa                                                                                                                                      >  8:428515/1‑137 (MQ=255)                      cTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGATTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATaa                                                                                                                                      <  7:428515/137‑1 (MQ=255)                         gAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTCGCTATAAACCTCACGCATTCGCGCTGGGCGTTGTTGACGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTc                                                                                                                         <  3:390548/147‑1 (MQ=255)                                cccTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGgcg                                                                                                                 >  4:376144/1‑148 (MQ=255)                                   tcaacaacaTGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtg                                                                                                              >  4:267895/4‑148 (MQ=255)                                    cGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgc                                                                                                             >  8:336251/1‑148 (MQ=255)                                    cGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgc                                                                                                             >  7:282454/1‑148 (MQ=255)                                      acaacaTGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCaaaa                                                                                                                                >  6:59765/1‑127 (MQ=255)                                      acaacaTGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCaaaa                                                                                                                                <  5:59765/127‑1 (MQ=255)                                         acaTGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgccgtgc                                                                                                        <  7:363481/148‑1 (MQ=255)                                            tGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCTCCAATCGGTGAACACGTGCACTGGAGCATGCTCGTTGATAACCAAAAATTGTTCCGCTGGCTCTGCCGTgccgc                                                                                                     >  2:246018/1‑148 (MQ=255)                                                       ccACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGcc                                                                                           >  6:46154/1‑147 (MQ=255)                                                       ccACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACATACGcc                                                                                           >  7:433698/1‑147 (MQ=255)                                                        cACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCaa                                                                                         >  7:196080/1‑148 (MQ=255)                                                        cACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCa                                                                                          >  6:252606/1‑147 (MQ=255)                                                         aCTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgccg                                                                                                           <  5:192581/129‑1 (MQ=255)                                                         aCTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgccg                                                                                                           >  6:192581/1‑129 (MQ=255)                                                           tgctgACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTg                                                                                      >  1:130684/1‑148 (MQ=255)                                                             ctgACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTgg                                                                                     <  3:45363/147‑1 (MQ=255)                                                               gACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCg                                                                                  >  3:130189/1‑148 (MQ=255)                                                                   gAAGGCTTTAGCTATAAACGTCACGCATTCGCGCTGGGCTTTGTTTAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAAACAAAAATTGTTCCGCTGTCGCTGCCGTTCCGCCACCTACGCCAACTGGCCGgtgt                                                                              <  1:246018/148‑1 (MQ=255)                                                                   gAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGgtgt                                                                              <  4:167056/148‑1 (MQ=255)                                                                    aaGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGtt                                                                             >  5:167397/1‑148 (MQ=255)                                                                     aGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTg                                                                            >  7:382534/1‑148 (MQ=255)                                                                      ggCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTGGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGc                                                                           <  5:128000/148‑1 (MQ=255)                                                                        cTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGt                                                                         >  2:197111/1‑148 (MQ=255)                                                                          ttAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTa                                                                       >  2:201861/1‑148 (MQ=255)                                                                           tAGCTTTAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTa                                                                       <  4:209385/147‑1 (MQ=255)                                                                            aGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTAc                                                                      <  4:130189/147‑1 (MQ=255)                                                                            aGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTAc                                                                      <  3:267895/147‑1 (MQ=255)                                                                              cTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACAt                                                                    <  8:196080/147‑1 (MQ=255)                                                                                ataAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACAtgct                                                                 >  6:148389/1‑148 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGTATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGc                                                          >  6:227785/1‑147 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGTATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGc                                                          >  2:368784/1‑147 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgccgtg                                                                                                         >  6:339680/1‑100 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCtgccgtg                                                                                                         <  5:339680/100‑1 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGt                                                                         >  7:340921/1‑132 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGt                                                                         <  8:340921/132‑1 (MQ=255)                                                                                        cACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGc                                                          >  3:168671/1‑147 (MQ=255)                                                                                           gCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGTATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATa                                                      <  5:227785/148‑1 (MQ=255)                                                                                           gCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGTATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATa                                                      <  1:368784/148‑1 (MQ=255)                                                                                              ttCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCg                                                   >  3:404097/1‑148 (MQ=255)                                                                                              ttCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCg                                                   <  1:201861/148‑1 (MQ=255)                                                                                               tCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAATTGGCCGGTGTTGCGTTACATGCTGCGCGGCa                                                         >  6:75315/1‑141 (MQ=255)                                                                                               tCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAATTGGCCGGTGTTGCGTTACATGCTGCGCGGCa                                                         <  5:75315/141‑1 (MQ=255)                                                                                               tCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGt                                                  >  4:406495/1‑148 (MQ=255)                                                                                                 gcgcTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGttt                                                <  2:443127/148‑1 (MQ=255)                                                                                                 gcgcTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGtt                                                 <  5:46154/147‑1 (MQ=255)                                                                                                   gcTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCt                                              <  3:376144/148‑1 (MQ=255)                                                                                                    cTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTg                                             <  4:404097/148‑1 (MQ=255)                                                                                                         cTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCa                                        <  2:130684/148‑1 (MQ=255)                                                                                                          tttGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCAc                                       >  8:415894/1‑148 (MQ=255)                                                                                                                gAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTg                                  <  5:148389/147‑1 (MQ=255)                                                                                                                     gcCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATc                            <  8:282454/148‑1 (MQ=255)                                                                                                                      ccACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCg                           >  4:114698/1‑148 (MQ=255)                                                                                                                       cACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACa                                                                     >  8:381923/1‑105 (MQ=255)                                                                                                                       cACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACa                                                                     <  7:381923/105‑1 (MQ=255)                                                                                                                        aCGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCgg                          >  3:393720/1‑147 (MQ=255)                                                                                                                                 aGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGAt                >  1:255388/1‑148 (MQ=255)                                                                                                                                  gACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATc               <  4:393720/148‑1 (MQ=255)                                                                                                                                  gACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAAAACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATc               <  7:336251/148‑1 (MQ=255)                                                                                                                                    cGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATccc             <  1:197111/148‑1 (MQ=255)                                                                                                                                      tGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTg           <  8:382534/148‑1 (MQ=255)                                                                                                                                      tGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTg           <  1:211371/148‑1 (MQ=255)                                                                                                                                          cTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCgg                          >  3:242193/1‑129 (MQ=255)                                                                                                                                          cTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCgg                          <  4:242193/129‑1 (MQ=255)                                                                                                                                            ggAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCt                    >  1:165524/1‑133 (MQ=255)                                                                                                                                            ggAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCt                    <  2:165524/133‑1 (MQ=255)                                                                                                                                            ggAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTc     >  8:323826/1‑148 (MQ=255)                                                                                                                                                cATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTg  >  6:221966/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                                   ACGTTTGATTCGCTGGCAATGCTGGAATTTGCCCTCGACAACATGCCGGATACCCCACTGCTGACCGAAGGCTTTAGCTATAAACCTCACGCATTCGCGCTGGGCTTTGTTGAAGCGCCACGCGGTGAAGACGTGCACTGGAGCATGCTCGGTGATAACCAAAAATTGTTCCGCTGGCGCTGCCGTGCCGCCACCTACGCCAACTGGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTG  >  NC_000913/2609756‑2610045

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