Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	40.9%	96.2 / 41.7	66	intergenic (‑108/+33)	uspF/ompN	stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (37/2);  new base A (0/27);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.66e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAT  >  NC_000913/1435579‑1435876                                                                                                                                                      |                                                                                                                                                      gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                                     <  2:110379/151‑1 (MQ=255)    tCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaaaa                                                                                                                                                   >  1:753683/1‑150 (MQ=255)      ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                                                >  1:12842/1‑151 (MQ=255)        aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                              >  2:356013/1‑151 (MQ=255)        aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                              >  2:796129/1‑151 (MQ=255)             tCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCAt                                                                                                                                         >  2:994512/1‑151 (MQ=255)                 tAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAc                                                                                                                                      >  2:635019/1‑150 (MQ=255)                    ggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCg                                                                                                                                  >  2:101009/1‑151 (MQ=255)                      aaaaaaatgtttttGTTAAAAAACCCTCCTGTTTTAGTATCCGAAAAAAGTGTAACGCCAGAGGACATTTTTTGTGAAATGACGGATTTCACTTTTTTAATTT‑AGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCgg                                                                                                                                 <  2:797993/137‑1 (MQ=255)                           aaTCTTTCTTTTAAAAAACCCCCTTTTTTTATTACCCGAAAAAAGTGAACCGCCAAATGAAATTTTTTGTGAAAGGAGGGAGTCCACATTTTAATTTA‑AAACCAAAGGAGAAAGAAAAAGAAAAAAAGCCCCGCCATAACAGCGGGCagga                                                                                                                           <  2:128049/151‑1 (MQ=255)                            aTCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                                            >  1:713169/1‑149 (MQ=255)                             tCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                                         >  1:1109159/1‑151 (MQ=255)                              cGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                                         >  2:325079/1‑150 (MQ=255)                               gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                       >  1:324832/1‑151 (MQ=255)                               gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                       >  2:1016287/1‑151 (MQ=255)                                ttCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAtt                                                                                                                      >  2:834587/1‑151 (MQ=255)                                   ttttCAAAAACCCTCTTTTTTTATTACCCGAAAAAAGTGTACCGCAAGATGACATTTTTTGTTAAAGGCCGGATTTCACTTTTTTAATTA‑AAAAAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                                         <  2:252071/143‑1 (MQ=255)                                     ttaaaaaacCCTCTTTTTTTAGTACCCAAAAAAAGTGAACCGCAAGATGACTTTTTTTGTGAAAGGACGGATTCCACTTTTTTAATTT‑AGAAAAAAGGAGAAAAAAAAAGCAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                                   <  2:538930/144‑1 (MQ=255)                                         tAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                              >  2:185845/1‑150 (MQ=255)                                          aaacCCTCTTTTTTTAGTACCCGAAAAAAGTGAAACCCCAGATAACTTTTTTTGTGAAATGACGGAGTTCACTTTTTAAATTA‑AGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGaa                                                                                                                 <  2:797072/146‑1 (MQ=255)                                                   tGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                                                                   >  2:545833/1‑151 (MQ=255)                                                     ttttAGTATCCGAAAAAAGTGTAACGCCAGAGGACTTTTTTTGTTAAATGACGGGGTTCACATTTTTATTTA‑AGATAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa                                                                                                          <  2:396678/142‑1 (MQ=255)                                                      tttAGTATCCGCAAAAAGTGTAACGCCAGATGACATTTTTTGTGAAATGACGGAGTTCACATTTTTAATTT‑AGAAAAAAGGAGAAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                 <  1:415607/150‑1 (MQ=255)                                                         aGTATCCCAAAAAAGTGTAACGCCAAATGACTCTTTTTGTGTAATGACGGAGTTCAATTTTTTAATTT‑AAAACAAAGGAGGAAAAAAAAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa                                                                                               <  1:973979/149‑1 (MQ=255)                                                          gttccctaaaaaaGTGAAACCCCAGATGACTTTTTTTGTGAAAAGACGGTGTTCACTTTTTTATTTT‑AAATCAAAGGAGAAAGAATAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                                                                              <  1:736590/139‑1 (MQ=255)                                                              ccGTATAATGTGCAACGCGAAATGTCATTTTTGGTTGATAGAGGTGGTTCATATTTTTAATTT‑AGAAAAAAGTAAAAAAAATAAGAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                                                                        <  1:716185/151‑1 (MQ=255)                                                               cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                                                                        >  2:232144/1‑150 (MQ=255)                                                               cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                                                                        >  1:34179/1‑150 (MQ=255)                                                                 cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                      >  1:268216/1‑150 (MQ=255)                                                                  aTAAAGTGTACCGCCAAATGACTTTTTTTGTGAAATGACGGATTTCACTTTTTTAATTT‑AGATAAAGGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAata                                                                                     <  1:786051/150‑1 (MQ=255)                                                                  aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                                    >  1:136806/1‑151 (MQ=255)                                                                    aaagttaaccGCCAGATGCCTTTTTTTGTTTAAGGACGGAGTTCACTTTTTTAATTA‑AGACAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatca                                                                                  <  1:670073/142‑1 (MQ=255)                                                                    aaaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatc                                                                                   >  1:587382/1‑150 (MQ=255)                                                                    aaaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatc                                                                                   >  2:740957/1‑150 (MQ=255)                                                                        tgtAACGCAAGATGACCTTTTTTTTATAATGACGGTGTT‑AAATTTTTAATTTAAAAAAAAAGGAAGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                               <  1:754009/150‑1 (MQ=255)                                                                           aaCCCCAGATGACATTTTTTGTGAAAGGAGGGAGTTCACATTTTTAATTT‑AGAAAAAAGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                           <  2:495325/151‑1 (MQ=255)                                                                              gCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                           >  1:61614/1‑148 (MQ=255)                                                                               ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                           >  2:1075288/1‑147 (MQ=255)                                                                                 agacaacatTTTTTGTGTAATGAGGGGTTTCAATTTTTTAATTA‑AGAAAAAAGGAAAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGAt                                                                       <  2:338725/144‑1 (MQ=255)                                                                                   aTGACTTTTTTTGTATAATGCCGGAGTTCATTTTTTTAATTT‑AGAACAAAGGAGAAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                                      <  1:46168/148‑1 (MQ=255)                                                                                   aTGACATTTTTTGGGTAATGAGGGAGTTCACATTTTTAATTT‑AAAAAAAAGGAGGAAAAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc                                                                    <  1:797911/150‑1 (MQ=255)                                                                                      acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                                                                >  2:1081478/1‑151 (MQ=255)                                                                                       ctttttttGTGTAATGGGGGAGTTCACATTTTTAATTA‑AGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                                                 <  2:947403/146‑1 (MQ=255)                                                                                          tttttGTGTAAAGCCGGAGTTCATTTTTTTAATGG‑AGAAAAAAGGGGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCa                                                             <  2:114638/150‑1 (MQ=255)                                                                                           ttttGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                           >  1:36895/1‑151 (MQ=255)                                                                                             ttGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                         >  2:873204/1‑151 (MQ=255)                                                                                                         ggAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCAt                                             >  2:97395/1‑151 (MQ=255)                                                                                                          gagttccttttttaAATTT‑AGAAAAAGGGAGAAAGAAAAAGAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTc                                               <  2:285056/139‑1 (MQ=255)                                                                                                           aGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                           >  1:1095839/1‑151 (MQ=255)                                                                                                             ttCACATTTTTAATTT‑AGACAAAAGGAGGAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                           <  1:120338/149‑1 (MQ=255)                                                                                                               cactTTTTAAATTT‑AGAAAAAAGGAGGAAGAATAAGGAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                       <  1:487824/151‑1 (MQ=255)                                                                                                                   tttttAATTA‑AGAAAAAGGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt                                          <  2:144194/144‑1 (MQ=255)                                                                                                                       tAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                               >  2:602734/1‑151 (MQ=255)                                                                                                                        aaTTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                              >  2:658869/1‑151 (MQ=255)                                                                                                                        aaTTC‑AGATAAAGGGAAAAAAAAGAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACa                                 <  2:307415/148‑1 (MQ=255)                                                                                                                         aTTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                             >  1:277800/1‑151 (MQ=255)                                                                                                                              agaaaaagggaggaAGAAAAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGa                           <  1:398751/141‑1 (MQ=255)                                                                                                                               gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt                        >  2:661914/1‑151 (MQ=255)                                                                                                                               gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt                        >  1:150935/1‑151 (MQ=255)                                                                                                                                 tCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                       >  1:721446/1‑150 (MQ=255)                                                                                                                                   aaaGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                     <  2:843881/150‑1 (MQ=255)                                                                                                                                     aggaggAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt                  >  2:724678/1‑151 (MQ=255)                                                                                                                                      ggaggaAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt                  <  1:785618/150‑1 (MQ=255)                                                                                                                                           aaaaaTAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAAc              <  1:101033/146‑1 (MQ=255)                                                                                                                                                     aGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAt  >  1:917209/1‑151 (MQ=255)                                                                                                                                                     aGAAAAAGCCCGCCATAACAGCGGGCAGGAGGACTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAt  >  1:543309/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAT  >  NC_000913/1435579‑1435876

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