Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,851,654	0	T	G	39.8%	133.2 / 29.7	84	G218G (GGT→GGG)	hypB	GTP hydrolase involved in nickel liganding into hydrogenases
Reads supporting (aligned to +/- strand):  ref base T (7/43);  new base G (33/0);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.25e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCT  >  NC_000913/2851505‑2851803                                                                                                                                                        |                                                                                                                                                      gACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGt                                                                                                                                                       <  2:538805/150‑1 (MQ=255)    gCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGa                                                                                                                                                     <  2:932271/149‑1 (MQ=255)    gCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAg                                                                                                                                                   <  2:237289/151‑1 (MQ=255)        aGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACaa                                                                                                                                               <  1:928234/151‑1 (MQ=255)         gatgatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGCAGGGGAAGACaa                                                                                                                                               >  1:366490/1‑150 (MQ=255)         gatgatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACaaa                                                                                                                                              <  1:198703/151‑1 (MQ=255)           tgatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGA‑GCTTT‑CCGTTTCCGAAGGGGGAAACaaaac                                                                                                                                            >  1:644418/1‑149 (MQ=255)            gatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCa                                                                                                                                           <  2:422382/151‑1 (MQ=255)                gCCGACGCCGCCCCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACt                                                                                                                                         <  1:221615/149‑1 (MQ=255)                gCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGa                                                                                                                                       <  2:820584/151‑1 (MQ=255)                gCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGG‑GCGTT‑TCCTTACCCGAGGGGGAGAAAAAacac                                                                                                                                          >  2:289934/1‑148 (MQ=255)                 ccGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGaa                                                                                                                                      <  2:79052/151‑1 (MQ=255)                   gACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGaaa                                                                                                                                     <  2:586266/150‑1 (MQ=255)                     cgccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAGGGGAAAACAAACCACTGAAata                                                                                                                                   >  1:1028875/1‑150 (MQ=255)                      gccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAatat                                                                                                                                  <  2:208917/150‑1 (MQ=255)                       ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAGGGGAAGACAAAACACTGAAata                                                                                                                                   >  1:1382025/1‑148 (MQ=255)                       ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGAAGGGGGaaacaaaacaacggaata                                                                                                                                   >  2:489307/1‑140 (MQ=255)                       ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATcc                                                                                                                                <  1:894169/151‑1 (MQ=255)                       ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATcc                                                                                                                                <  1:1008991/151‑1 (MQ=255)                        cgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCg                                                                                                                               <  1:910801/151‑1 (MQ=255)                         gcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGc                                                                                                                              <  2:498516/151‑1 (MQ=255)                            ccGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GGCTT‑TCGTTTTCGAATGGGGAGAAAAACCCCTGAAATAACCg                                                                                                                               >  2:343569/1‑147 (MQ=255)                              gcgTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCatat                                                                                                                          <  1:1247728/150‑1 (MQ=255)                               cgTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCCTTACCGAAGGGGAAGACAAACCACCGAAAAAACCCCATATGt                                                                                                                        >  1:953455/1‑151 (MQ=255)                               cgTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATg                                                                                                                         <  1:813562/150‑1 (MQ=255)                                gTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGtt                                                                                                                       <  1:378259/151‑1 (MQ=255)                                        cTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGGAGGGGAAGAAAAAACACTGAAAAATCCGCATATGTTTgccgcc                                                                                                                >  1:564367/1‑150 (MQ=255)                                        cTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTTCCGGAGGTGAAGAAAAACCACCGAAAAATCCGCCTATGTTTtccgcc                                                                                                                >  1:178433/1‑150 (MQ=255)                                          ggACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGAAGGGGAAGACAAACCACCGAAATATCCGCATATGTTTgccgccgcc                                                                                                             >  1:100786/1‑151 (MQ=255)                                          ggACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTgccgccgcc                                                                                                             <  2:387885/151‑1 (MQ=255)                                          ggACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTgccgccgc                                                                                                              <  2:144819/150‑1 (MQ=255)                                          ggACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGG‑GGCTT‑CCGTTTCCGACGGGGAAAAAAAACAACCCAAAAATCCCCCTATTTTTgcc                                                                                                                   >  2:187035/1‑145 (MQ=255)                                              gATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGGGAAGACAAACCACTGAAATATCCGCATATGTTTtccgccgcccc                                                                                                           >  1:1094727/1‑147 (MQ=255)                                                 aaCGGTATTCTGTTTATCGAAAACGTTGGCCACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGa                                                                                                      <  1:1164778/151‑1 (MQ=255)                                                 aaCGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAAGGGAAGACAAACCACTGAAAAATCCCCATATGTTTTCCGCCGCCCCGCTg                                                                                                       >  1:730679/1‑150 (MQ=255)                                                  aCGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAGGGGGAGACAAAACAATGAAATATCCGCATATTTTTTCCGCCGCCTCGCTGa                                                                                                      >  2:1150901/1‑150 (MQ=255)                                                   ctgTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCCAAAGGGAAGAAAAACCACTGAAAAATCCGCCTATTTTTTCCGCCCCCTCCCTGa                                                                                                      >  1:1347187/3‑149 (MQ=255)                                                    ggTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGAtgc                                                                                                   <  2:18008/151‑1 (MQ=255)                                                       aTTCCGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGAtgct                                                                                                  <  1:984635/149‑1 (MQ=255)                                                          cTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTTCCGAAGGGGAAGACAAACCACTGAAATATCCGCGTATGTTTTGCCCCGCCTCGCCGATGCTGCTc                                                                                              >  1:3795/1‑150 (MQ=255)                                                              ttATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTTCCGAAGGGGAAGACAAAACACTGAAATATACGCATATGTTTGTCGCCGCCTCGCTGATGCTGCTcaaca                                                                                          >  1:80144/1‑150 (MQ=255)                                                               tATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGG‑GCTTTACCGTT‑CCCAAGGAGAAAACAAAACACTGAAAAAACCGCATATGTTTGCCGCCCCCCCCCCGGTGCTGCTccacaa                                                                                         >  1:978467/1‑150 (MQ=255)                                                                  cGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGGAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                                                                                     <  2:1003441/151‑1 (MQ=255)                                                                  cGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGt                                                                                      <  2:173218/150‑1 (MQ=255)                                                                   gAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAAGGGAAAAAAAAACACCGAAATATTCGCATATGTTTGCCGCCCCCCCCCTGATGCTGCCCAACCAAGtt                                                                                     >  2:410990/1‑150 (MQ=255)                                                                     aaaCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GGTTT‑CCGTTACCGAAGGGGGAGACAAAACACCGAAATATTCGCATATTTTTGCCCCCGCCTCCCTGGTGCTGGTCAACaaaagt                                                                                     >  1:117427/1‑145 (MQ=255)                                                                      aaCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAAGGGAAGACAAACCACTGAAATATTCGCATATGTTTGTCGCCGCCTCCCTGGTGCTGCTCCAAAAAAATGAcc                                                                                 >  2:803954/1‑151 (MQ=255)                                                                       aCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCt                                                                                <  1:289092/151‑1 (MQ=255)                                                                         gTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GGTTT‑CCCTTTACCCAGGGGGAGAAAAAACCCTGAAAAATCTGCCTTTATTTGCCGGCGCCCCCCCGATGCCGCCCaaaaaaagttacc                                                                                 >  2:668389/1‑141 (MQ=255)                                                                         gTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCCTTTCCGAAGGGGAAAACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCCCAACCAAAATGACCTGt                                                                              >  2:580293/1‑151 (MQ=255)                                                                            ggCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTg                                                                            <  2:627043/150‑1 (MQ=255)                                                                            ggCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGc                                                                           <  1:617622/151‑1 (MQ=255)                                                                             gCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGcc                                                                          <  2:1213527/151‑1 (MQ=255)                                                                               aaCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGGTTTT‑CCTTTCCCAAAGGGGAA‑ACAAACCACCGAAAAATCCGCATATGTTTTCCGCCCCCCCGCTGATTCTGTTCAACAAAAATGACCCTTTGcc                                                                          >  2:148455/1‑149 (MQ=255)                                                                               aaCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGAAGGGGGAGAAAAAACACTGAAATATTCGCATATGTTTGCCGCCGCCCCCCTGATGCTGCTCCACAAAAATGACCTGTTTcc                                                                          >  2:552974/1‑149 (MQ=255)                                                                                aCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTa                                                                       <  1:976042/151‑1 (MQ=255)                                                                                     gTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTAtc                                                                     <  1:289941/148‑1 (MQ=255)                                                                                     gTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCa                                                                  <  1:614485/151‑1 (MQ=255)                                                                                      tATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCCTT‑CCGCTATCGGAGGGGAAGACAAAACACCGAAATATCCCCCTATTTTTTCCGCCGCCCCCCTGATGGTGCTTAACAAAATTGACCCGTTGCCCGAtctc                                                                   >  1:142934/1‑149 (MQ=255)                                                                                            cGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTACCGAAGGGGAAGACAAAACACTGAAATATCCGCATATGTTTTCCGCCGCCCCGCTGATGGTGCTCAACAAAATTGACCTGTTGCCGGATCTCAACttt                                                             >  1:637494/1‑149 (MQ=255)                                                                                                cAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGtt                                                       <  2:1347114/151‑1 (MQ=255)                                                                                                 aGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGGGAAGACAAACCACTGAAAAATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGtt                                                       >  1:516151/1‑150 (MQ=255)                                                                                                       gATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAAGGGGAGACAAAACAACGAAAAATACGCAAATGTTTTCCCCCCCCCCGCCGATGCTGCCCCACAAAATTTTACTGTTGCCGGATCTCAACTTTTACGTTGaaaa                                                  >  2:1291921/1‑149 (MQ=255)                                                                                                          ctcGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCCTTAACGAAGGGGGAGACAAACCACTGAAAAAACCGCATATGTTTGCCGCCGCCTCGCCGGTGCTGCTCAACAAAATTGACCTGTTGCCGGATCTCCACTTTGACGTTGaaaa                                                  >  2:416901/1‑146 (MQ=255)                                                                                                          ctcGGTGAAAAACACAAAGTGGCGGT‑CCTTT‑TCTTTAACCGAGGTGGAGACAAACCACTGAAATATTCGCATACGCTTTCCGCCGCCCCGCCCATGCTTTTCCAAAAAATTGACCAGTTTCCGTATGTCGAAttttatgattcaaag                                                 >  2:1225074/1‑135 (MQ=255)                                                                                                          ctcGGTGAAAAACACAAAGTGGCGGA‑GCTTT‑CCGTTACCCAAGGGGAAAACAAAACACTGAAAAATTCGCATATTTTTTCCCCCGCCCCGCTGGTGCTGCCCAAAAAAGTTGACCTGTTGCCGTATCTCCACTTTTACGTTGaaaa                                                  >  1:811153/1‑146 (MQ=255)                                                                                                           tcGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTAACCGAGGGGGAGACAAAACACCGAAAAATCCGCATATGTTTTGCGCCGCCTCCCCGATGCTGCTCAACAAAATTGAACCCTTGCCCTATCTCCACCTTTTCGTTGtaaa                                                  >  2:1047745/1‑145 (MQ=255)                                                                                                            cGGTGAAAAACACAAAGTGGCGGT‑GGTTT‑TCGGTACCGAAGGTGGAGACAAACAACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCCGATTCTGCTCAACAAAGTTGAACTGTTGCCGGATCTCCACTTTGACGTTGAAAAgg                                                >  1:773647/1‑145 (MQ=255)                                                                                                             ggTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCg                                          <  1:1307151/151‑1 (MQ=255)                                                                                                             ggTGAAAAACACAAAGTGGCGGC‑GCTTT‑CCGTTACCGAAGGTGAAGAACAAACACAGAAATATTCGCATATGTTTTTCGCCGCCCCCCTGATTCTGCTCAACAAAACTTACCCTTTGTCGTATCTCCACTTTTGTGATGAAAAGAGCAt                                            >  1:1174784/1‑149 (MQ=255)                                                                                                               tGAAAAACACAAAGTGGCGGT‑GCTTT‑TCGTTTCCGAAAGGGAAGAAAAAACACTGAAATATACGCATATGTTTTCCGCCGCCCCGCTGATGCTGCCCCACAAAATTGACCTGTTTCCGTATCTCCACTTTTACGTTGAAAAAGGCATCGcc                                        >  2:609699/1‑151 (MQ=255)                                                                                                                gAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCAGCg                                          <  2:258869/148‑1 (MQ=255)                                                                                                                      cacaAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGccc                                 <  1:1142177/151‑1 (MQ=255)                                                                                                                       acaAAGTGGCGGT‑GCTTT‑CCGTTTCCGAAGGGGAAGACAAAACACTGAAATATCCGCATATGTTTTCCGCCCCCCCGCTGGTGCTGCTCAACAAAGGTGGCCTGTTGCCGTATCTCAAATTTTACGTTGGAAAAAGCATAGCCCGCGccc                                 >  2:139099/1‑150 (MQ=255)                                                                                                                         aaaGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCcgcg                              <  2:435493/151‑1 (MQ=255)                                                                                                                                cGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCa                        <  2:280026/150‑1 (MQ=255)                                                                                                                                 ggT‑GCTTT‑CCGTTACCGAAGGTGAAGACAACCCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAAt                      <  1:1208417/151‑1 (MQ=255)                                                                                                                                 ggT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAAt                      <  1:731623/151‑1 (MQ=255)                                                                                                                                  gT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATc                     <  1:1334102/151‑1 (MQ=255)                                                                                                                                     gCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCa                   <  2:509812/151‑1 (MQ=255)                                                                                                                                      cTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATcc                    >  1:906147/1‑149 (MQ=255)                                                                                                                                         t‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGaa                >  2:256364/1‑150 (MQ=255)                                                                                                                                            cGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGaa                <  2:947836/148‑1 (MQ=255)                                                                                                                                             gTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTg            >  2:157086/1‑151 (MQ=255)                                                                                                                                               tACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCTGAAATTGa           <  1:158592/150‑1 (MQ=255)                                                                                                                                                        tGAAGACAATAAAATGAAATATCCGCATATGTTTGGGGATGCCTCGCTGATGCTGCTCAACAAAGGTGAAGTGTTGCGGCATCTAAACTTTGACGGTGAAATATGTATCGCCTGGGCCCGCGAAGTCAAGCCAGAAATTGAAATCATCCt  >  2:1086181/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGT‑GCTTT‑CCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCT  >  NC_000913/2851505‑2851803

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