Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	266,358	0	T	G	41.3%	111.7 / 14.1	64	T66P (ACC→CCC)	yafY	lipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (7/30);  new base G (26/0);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266218‑266499                                                                                                                                               |                                                                                                                                              gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTTGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGGAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAT‑CCTG‑CAAGTTGGGACCgg                                                                                                                                       <  2:352221/149‑1 (MQ=255)    ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAA‑ACCG‑CCAGGGGGGGCCGGATaa                                                                                                                                   >  1:206736/1‑150 (MQ=255)      ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑ACTG‑CAAGGGGGGACCGGATaa                                                                                                                                   >  1:869316/1‑148 (MQ=255)       cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAagtag                                                                                                                               <  2:735888/151‑1 (MQ=255)           ttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCa                                                                                                                            <  1:779292/150‑1 (MQ=255)             gATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑ACTG‑CAAGGGGGGGCCGGGAAAGTAGGCAAcc                                                                                                                         >  1:95868/1‑151 (MQ=255)              aTGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCt                                                                                                                        <  2:366473/151‑1 (MQ=255)               tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑ACCG‑CAAGGGGGGGCCGGGTAAGTAGGCaaac                                                                                                                         >  1:187260/1‑147 (MQ=255)                gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCtgt                                                                                                                      <  1:566349/151‑1 (MQ=255)                gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGGGGGGCCGGATAAGTAGTCAAcc                                                                                                                         >  1:146273/1‑148 (MQ=255)                gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CCACGGGGGCGCCGGGAAGTAGGCAACCtgt                                                                                                                      >  2:557542/1‑151 (MQ=255)                  tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                     <  2:498215/150‑1 (MQ=255)                  tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                     <  2:378181/150‑1 (MQ=255)                  tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                     <  1:787154/150‑1 (MQ=255)                   aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGAc                                                                                                                   <  1:506196/151‑1 (MQ=255)                     cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGGGGGGCCGGATAAGTAGTCAACCTGTGAcc                                                                                                                  >  1:124949/1‑150 (MQ=255)                     cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGGGGGGCCGGAAAAGTAGTCAACCTGTGAcc                                                                                                                  >  1:585908/1‑150 (MQ=255)                     cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTC‑CAAGGGGGGAGCGGGTaaagaggccaac                                                                                                                         >  2:435893/1‑133 (MQ=255)                         aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttct                                                                                                             <  2:852912/151‑1 (MQ=255)                             cacaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CAAGGGGGGGCCGGATAAATAGTCAACCTGTGGCCATTCTTCGc                                                                                                         >  2:166782/1‑151 (MQ=255)                              acaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCt                                                                                                        <  1:803766/151‑1 (MQ=255)                                aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGa                                                                                                      <  1:433897/151‑1 (MQ=255)                                 gcgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAAGCTGTGACCATTCTTCGCTg                                                                                                       <  2:932546/149‑1 (MQ=255)                                     aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CAAGGGGGGGCCGGGTAAGTAGTCAACCTGTGACCATTCTTTGCTGAc                                                                                                     >  1:931363/1‑147 (MQ=255)                                      ttCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑CCTG‑CAAGGGGGGACCGGATAAGTAGTCAACCTGGGACCCTTCTTCGCTGAcatca                                                                                                 >  2:881963/1‑150 (MQ=255)                                         atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑ACCG‑CAAGGGGGGGCCGGATAAAAAGGCAAACTGTGGACATTCTTCGCTGACAtcctcc                                                                                              >  2:358753/1‑150 (MQ=255)                                            aGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CAAGGGGGGACCGGATAAGTAGTCAAACCGTGACCATTCTTCGCTGACATCATcc                                                                                              >  1:396172/1‑147 (MQ=255)                                               ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGCAAAA‑CCTGACAAGCGGGGGCGGGA‑AAGGAGGCAACCTATGTCCATTCTTCCCTGAAAACAACCttgtg                                                                                         >  1:634607/1‑149 (MQ=255)                                                     aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACca                                                                                 <  1:820355/151‑1 (MQ=255)                                                        gTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCCG‑CCAGGGGGGGCCCGGTAAGAAGGCAACCTGTGACCATTCTTTGCTTACATCATCCTGGTGGCCGAcccc                                                                                >  2:418303/1‑147 (MQ=255)                                                          tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAA‑CCTG‑CAAGGTGGGACCGGATAAGGATTAAACCTGTGACCATTCTTCTCTGACATCATCCTTGTGGTTGAACACATc                                                                             >  1:249237/1‑150 (MQ=255)                                                           ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCa                                                                            <  2:69064/150‑1 (MQ=255)                                                            tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAc                                                                           <  1:58637/150‑1 (MQ=255)                                                              tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACCG‑CCAGGGGGGGCCGGGGAAGTAGTCAACCCGTGACCATTCTTTCCCGACATCATCCTGGGGGCTGACcacaac                                                                             >  2:791821/1‑144 (MQ=255)                                                               ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGGGGGGCCCGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGGCCACATCAc                                                                           >  2:420733/1‑147 (MQ=255)                                                                cTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                       >  1:237605/1‑150 (MQ=255)                                                                   cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                                   <  2:537114/151‑1 (MQ=255)                                                                   cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                                   <  1:31510/151‑1 (MQ=255)                                                                    atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                  <  1:589475/151‑1 (MQ=255)                                                                        ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAACCCGG‑CAGGGGGGGGCCGGA‑AAATAAGCAAACTGTGGCCATTCTTCGCCGACATAATCCTGGTGGCTGACCACATCACTGAAACCCGGAgg                                                              >  2:236176/1‑151 (MQ=255)                                                                            gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTg                                                          <  2:143078/151‑1 (MQ=255)                                                                            gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTg                                                          <  1:316226/151‑1 (MQ=255)                                                                                tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACCTAAAA‑CCTG‑CAAGGTGGGGCCGGATAAGTAGTCAACCTGTGCCCATTTTTCGCTGGCATCATTCTGGGGGGCGGCCCCATCACCGAAACGCGGTggggggggg                                                       >  2:144634/1‑145 (MQ=255)                                                                                 agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAA‑ACTG‑CCAGGTGGGGCCGGAGAAGTAGTCAACCTGTGACCCTTTTTTGGTGACAACCTTCCGGTGGCCGGCCACCTAACTGAAAAGCGggggggggggggc                                                     >  1:277523/1‑144 (MQ=255)                                                                                         cTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGAAAAA‑ACTG‑CAAGGAGGGACCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTTGTGGCTGACCACAACACTGAAAAGCGGTGGGGTGAGGGCCCTTAAtt                                             >  1:489844/1‑151 (MQ=21)                                                                                          tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGCCGGAAAA‑CCTG‑CCAGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTCGCTGACATCATCCCGGGGGCTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCCTTAAttt                                            >  2:224839/1‑151 (MQ=18)                                                                                               tCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTT‑CCAGGGGGGACCGGGTAAGGAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAAGCGGTGGGGGGGGGGCCTTTAATTTACt                                         >  2:194063/1‑149 (MQ=18)                                                                                                ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                       <  2:256450/150‑1 (MQ=255)                                                                                                ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                      <  1:915942/151‑1 (MQ=255)                                                                                                     cAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                 <  2:624883/151‑1 (MQ=255)                                                                                                      aGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                <  1:70174/151‑1 (MQ=255)                                                                                                      aGCCATCGGGCGACTGGCGGGGACCTAAAA‑CCTG‑CCAGGGGGGACCCGGTAAGAAGTCAACACGTGACCATTTTTTGCCGAAATCACACTGGGGGGTGGCCCCATCACTGAAACGCGGGGGGGGGAGGTGCTTtaa                                               >  2:443915/1‑136 (MQ=255)                                                                                                       gCCATCGGGCGACTGGCGGGGACGTAAAA‑ACTG‑CCAGGGGGGGCCGGATAAGAAGGCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGAACACATCACTTAAACGCGGTGGGGGGGGGGCCTTAAATTTACTGATgggcct                                >  2:528409/1‑146 (MQ=17)                                                                                                         cATCGGGCGACTGGCGGGTACGTAAAA‑ACTG‑CAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGGGGGGGCCTTTAATTTACTGATGGCCTTc                               >  1:808011/1‑149 (MQ=21)                                                                                                          aTCGGGCGACTGGCGGGGACGTAAAA‑CCTG‑CAAGGGGGGACCGGATAAGTAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGGGGGGGGAGGGCCTTAAATTTACTGATGGCCTTCAgg                            >  1:421853/1‑151 (MQ=18)                                                                                                            cGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGCGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGt                           <  1:386251/150‑1 (MQ=255)                                                                                                             gggCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                         <  2:514353/151‑1 (MQ=255)                                                                                                              ggCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                        <  2:158370/151‑1 (MQ=255)                                                                                                                cGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatca                      <  2:690864/151‑1 (MQ=255)                                                                                                                   cTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGGCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGGGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACtt                   >  2:250805/1‑151 (MQ=255)                                                                                                                            aCGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc           <  1:714522/150‑1 (MQ=255)                                                                                                                               tAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt       >  2:575942/1‑151 (MQ=255)                                                                                                                                     ccTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt  >  1:11426/1‑151 (MQ=255)                                                                                                                                        g‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt  >  2:927929/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                              GTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAA‑CCTG‑CAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266218‑266499

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