Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,039	0	T	G	42.6%	100.5 / 15.3	61	N108T (AAC→ACC)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (26/0);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184                                                                                                                                                    |                                                                                                                                                   aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGtttt                                                                                                                                                 <  1:503876/151‑1 (MQ=255)   ctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGgttt                                                                                                                                                 >  2:615253/1‑149 (MQ=255)     gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGc                                                                                                                                             <  2:771096/151‑1 (MQ=255)     gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGc                                                                                                                                             <  1:254978/151‑1 (MQ=255)     gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGGGGgttt                                                                                                                                                 >  1:1121359/1‑147 (MQ=255)       ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCGAGGCAGAGGGTTTTCGcc                                                                                                                                            >  1:322626/1‑150 (MQ=255)        tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCa                                                                                                                                           <  2:300493/150‑1 (MQ=255)         ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGAGGGTTTTCGCCAc                                                                                                                                          >  2:756549/1‑150 (MQ=255)                gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGTCCGAGGGAGAGGGTTTCTTCCCCTGTc                                                                                                                                      >  2:577283/1‑147 (MQ=255)                  gCTTAGCGTTGAAACAAGTTTTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                <  2:735664/151‑1 (MQ=255)                  gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                <  2:288148/151‑1 (MQ=255)                          ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                         <  2:1020437/150‑1 (MQ=255)                           tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCCGAGGGTTTCCGCCACTGTCAGCCCTTccccc                                                                                                                          >  2:932607/1‑148 (MQ=255)                           tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAg                                                                                                                        <  2:999719/150‑1 (MQ=255)                            gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGa                                                                                                                       <  1:615488/150‑1 (MQ=255)                            gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTCACCCAGGGAGAGGGTTTCCCCCACCCTCCGCCCCTCCCCCCGa                                                                                                                       >  2:5447/1‑150 (MQ=255)                             aaccAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCAGCGCCTTTTGCAGGTAATCACCCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAc                                                                                                                      <  1:299211/147‑1 (MQ=255)                                 aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCAACCAGGCGGAGGGGTTCCGCCCCCATCAGCCCCTCCACCCGACTTAAc                                                                                                                 >  2:859383/1‑151 (MQ=255)                                   gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTGCCCCCCCCCCCGGGAGGGGGTTTTTCCCCCTGTCAGTCCTTCCCCCAGaagcaa                                                                                                                  >  1:732918/1‑143 (MQ=255)                                      gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGCAGAGGGTTTTCGCCCCCGTCCGCCCCTCCCCCCGACTTAACTcc                                                                                                              >  2:190418/1‑149 (MQ=255)                                      gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGGAGAGGGTTTTCGCCACCGTCAGCCCCTCCCCCCCCCCTAACTcc                                                                                                              >  1:902073/1‑149 (MQ=255)                                       ttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCCGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACCTAACTcc                                                                                                              >  1:1223034/1‑148 (MQ=255)                                        tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                           <  1:655198/150‑1 (MQ=255)                                         cAGGACTAACGCCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCGGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTcc                                                                                                              <  1:831653/146‑1 (MQ=255)                                           ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                       <  1:1029466/151‑1 (MQ=255)                                            gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCCACTTAAATCCTGATaaa                                                                                                       >  2:213740/1‑150 (MQ=255)                                              cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAACCCCTCCCCCCGACCTAACTCCTGATaaa                                                                                                       >  2:1246047/1‑148 (MQ=255)                                               tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGCAGAGGGTTTTTCCCCCCGTTAGGCCCTCCACCCGACTTAACACCTCATaaaac                                                                                                     >  1:1087210/1‑147 (MQ=255)                                                aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACTTAAATCCTGATaaaaccc                                                                                                   >  2:706446/1‑146 (MQ=255)                                                aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                  <  2:1041228/151‑1 (MQ=255)                                                    tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCCGGGCAGAGGGTTTCTGCCCCCGTTAGCCCCTCCACCCGACATTAATCCTCGTAAAAcgcgcg                                                                                                >  1:248606/1‑149 (MQ=255)                                                      aaCGCCCAGCGCTTGTAAGCAACGTTGAGCATCCGGCGCCATTGGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCgcgcgc                                                                                               <  2:675548/148‑1 (MQ=255)                                                          cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCCCCCCCCTCAGCCCTTCCACCCGACATAAATCCTGATAAAACGCGCGCcccccc                                                                                         >  2:775370/1‑150 (MQ=255)                                                            cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCGGCGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTTCTGATAAAACACGCcccccccc                                                                                         >  2:153801/1‑148 (MQ=255)                                                                gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCCCCCCCGTCAGCCCTTCCACCCGACTTAACTCCTGATAAACCGCGCGCCCCCCCAGCTc                                                                                    >  2:457331/1‑149 (MQ=255)                                                                 ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAg                                                                                  <  2:672590/150‑1 (MQ=255)                                                                    gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCCCTGTCCGCCCTTCCCCCAGACCTAACTCCTGATAAAACGCGCGCCCCCCCCCTTTAGCgg                                                                               >  1:960939/1‑150 (MQ=255)                                                                      aggCAACGTTGGGCATCCTGCGCCATTTGCAGGGAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGAGAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                            <  1:203265/149‑1 (MQ=255)                                                                      aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTCACCGAGGCAGAGGGTTTCCCCCACCGTTAGCCCTTCCACCCCACTTAACACCTGATAAAACGCGCGCCCCCCCCGCTCAGCgg                                                                               >  2:256800/1‑148 (MQ=255)                                                                               tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAAACCCCCGCCCCCCCCGTTCAGCGGCAccg                                                                          >  2:955463/1‑142 (MQ=255)                                                                                gggCATCCTGCGCCATTTGCAGTTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                   <  1:437010/150‑1 (MQ=255)                                                                                gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCCACCGAGGGGGAGGGTTTTCGCCCCCGTCATCCCCTCCCCCCGACTTAACTCCTGATAAAACGCGCGCACCCCCCCCTCAGCGGCGCGGCGAgcc                                                                    >  1:278365/1‑147 (MQ=255)                                                                                gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                   <  2:900345/150‑1 (MQ=255)                                                                                    aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGAGGGAGAGGGTTTTCCCCCCTCTCCGACCCTCCCCCCCACCTAAACCCCCATAAAACCCGCGCCCCCCCCCTTCAGCGGCCCGGCggggcc                                                                   >  1:761970/1‑141 (MQ=255)                                                                                         gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGAACGAGGGCGGGGGTTTTCGCCCCCGTCCGCCCTTCCACCCGACTTAACTCCTGATAAAACGCCCCCCCCCCCCCGTCAGCGGCACCGCGGGGCAGCGGAGcct                                                          >  1:98353/1‑148 (MQ=255)                                                                                          cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACCGAGGCCGAGGGTTTTCGCCCCCCTCAGGCCCTCCACCACACTTCACTCCTGTTAAAACGCGCGCCCCCCCCGTTTATCCGCACGGCGAGGCAGCGGAGcgt                                                          >  1:1124676/1‑147 (MQ=255)                                                                                            ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑c                                                     <  2:1138761/151‑1 (MQ=255)                                                                                             cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑Ca                                                    <  2:906032/151‑1 (MQ=255)                                                                                                        aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTc                                          <  2:1134096/150‑1 (MQ=255)                                                                                                         atcaatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTcgc                                        <  2:108451/151‑1 (MQ=255)                                                                                                            attcaTGCCGTTGCGGCGGGGCCACTGCCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCGCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTcgcg                                       <  1:970231/147‑1 (MQ=255)                                                                                                            aatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGtc                                     <  2:6548/151‑1 (MQ=255)                                                                                                              tcaTGCCGTTGCGGCGGGGCCGCTGCCCCAGGGAGAGGTTTTTCGCCCCCGTCAGGCCTTCCCCCCCACTTCTCCCCTGCTTAAACCCGCGCACCCCCCGTTCTGCGCCACGGCGCGGCCGCGTAGCTTTGTCCACCTTccctcccgtgtc                                    >  1:579857/1‑141 (MQ=255)                                                                                                                  gCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGcc                                <  2:902121/150‑1 (MQ=255)                                                                                                                   ccGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCaa                              <  1:924621/151‑1 (MQ=255)                                                                                                                       tGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATcc                           <  2:79008/150‑1 (MQ=255)                                                                                                                        gcggcgGGGCCACTGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGGCTCGCCAATc                            >  2:246447/1‑148 (MQ=255)                                                                                                                                gCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                 >  1:1229726/1‑151 (MQ=255)                                                                                                                                   aCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc               >  1:948570/1‑150 (MQ=255)                                                                                                                                        ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg         >  1:692207/1‑151 (MQ=255)                                                                                                                                          gAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg       <  2:547907/151‑1 (MQ=255)                                                                                                                                           aGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc      >  2:496487/1‑151 (MQ=255)                                                                                                                                            ggCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc      <  1:538993/150‑1 (MQ=255)                                                                                                                                               agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc  >  1:403610/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGT‑CACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184

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