Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	39.3%	81.0 / 36.5	62	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 (30/7);  new base A (0/24);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 5.42e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATG  >  NC_000913/1435581‑1435877                                                                                                                                                    |                                                                                                                                                       aTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaaaa                                                                                                                                                    <  2:168920/151‑1 (MQ=255) aTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaaaa                                                                                                                                                    <  1:98208/151‑1 (MQ=255)   cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaaaaa                                                                                                                                                   <  1:505770/150‑1 (MQ=255)   cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGaaaaa                                                                                                                                                   <  1:1014340/150‑1 (MQ=255)    ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGCAGAATAAGCAGAAACAGc                                                                                                                                                 >  1:347562/1‑151 (MQ=255)      aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAACCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCATAAAAAGccc                                                                                                                                               <  1:668488/151‑1 (MQ=255)      aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                               >  1:146664/1‑151 (MQ=255)       aaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCg                                                                                                                                              <  1:85438/151‑1 (MQ=255)         tatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                            >  1:151051/1‑151 (MQ=255)                  gaaCAAGAATGGTTCTTTTAAAAAAACCTCCTTTTTTAGTATCCGAATAAGGTTAACCGAAAGATGACATTGTTTGAAAAATGAGGGAGTTCACATTTTTAATTT‑AGAGAAAGGGAGGAAGAATAACCAGAAAAAGCCCGCCATa                                                                                                                                         <  1:598785/143‑1 (MQ=255)                    acacGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTCACGCCAGATGACCCTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAACGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                                                 >  1:333158/1‑151 (MQ=255)                      aaGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                                >  2:677754/1‑150 (MQ=255)                         aaTCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                                             >  1:771226/1‑150 (MQ=255)                               tCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                      >  1:602398/1‑151 (MQ=255)                                      aaaaacCCTCCTTTTTTAGTACCCGAAAAAAGTGTAACGCCAGATGACTTTTTTTGTGAAATGACGGATTTCACTTTTTAAATTA‑AGACAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt                                                                                                                <  2:274687/148‑1 (MQ=255)                                      aaaaacCCTCCTGTTTTATTATCCGAAAAAAGTGTACCGCCAGATGACTCTTTTTGTGTAATGAGGGAGTTCAAATTTTTAATTT‑AGATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                               <  2:614791/149‑1 (MQ=255)                                        aaacCCTCCTGTTTTATTACCCGAATAAGGTTTAACGCCATATGACTTTTTTTGTGTAATGCCGTAGTTCA‑TTTTTTAATTTAAAATCAAAGGGGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa                                                                                                                  <  2:337594/146‑1 (MQ=255)                                        aaaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                              <  2:287135/150‑1 (MQ=255)                                          aaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                                                            >  2:330992/1‑150 (MQ=255)                                                  gTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                   >  1:833371/1‑151 (MQ=255)                                                   ttttATTTTCCGCAAAAGGTGAACCGCAGGATACTACTTTTTTTGTAAGGACGGAGTTCACTTTTTTATTTT‑ATAAAAACGAAGGAAAAATAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAg                                                                                                      <  2:865587/147‑1 (MQ=255)                                                   tttaAGTATCCGCAAAAAGTGAAACGCAAAATAACTTTTTTGGTGTAAGGACGGATTTCACTTTTTTAATTA‑AAACAAAAGGAGGAAAAAAAAAAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                   <  1:720867/150‑1 (MQ=255)                                                      tAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa                                                                                                >  2:391782/1‑150 (MQ=255)                                                          atccaaaaaaaGTTAAACGCAAGATGACTTTTTTTGTGAAATGAGGGATTTCAC‑TTTTTAATTTAAGATAAAAGGAGAAAAAAAAAAAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                                 <  1:158867/137‑1 (MQ=255)                                                          aTCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                            >  2:21110/1‑150 (MQ=255)                                                          aTCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                            >  2:897911/1‑150 (MQ=255)                                                           tCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                                          >  2:280708/1‑151 (MQ=255)                                                               cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                       >  2:435629/1‑150 (MQ=255)                                                                aTAAAGTGAACCGCCAGATGACTTTTTTTGTGTAATGACGGATTTCACATTTTTAATTT‑AGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                       <  2:786984/149‑1 (MQ=255)                                                                  aaggtgtAACGCCAGATAACACTTTTTGTGAAAGGGCGGAGTTCAGATTTTAAATTT‑AGAAAAAAGGAGAAAAAATAAGCAAAAAAGTGCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatca                                                                                   <  1:530362/148‑1 (MQ=255)                                                                         aaCGCCAGATGACTTTTTTTGTGAAATGACGGATTTCATTTTTTTAATTT‑AGAAAAAAGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                                                             <  2:291023/150‑1 (MQ=255)                                                                          acccaaGATAACATTTTTGTTGTAATGACGTAGTTCATTTTTTTAATTA‑AGAAAAAAGAAGAAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                              <  2:159078/143‑1 (MQ=255)                                                                              cAGATGACACTTTTTGCGTAATGACGGAGGTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCAGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                                       >  1:719670/1‑151 (MQ=255)                                                                                    acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                                                                 >  2:45945/1‑151 (MQ=255)                                                                                    acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                                                                 >  1:713149/1‑151 (MQ=255)                                                                                      acTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                                >  1:855456/1‑150 (MQ=255)                                                                                        tttttGTTTAATGACGGATTTCAATTTTTTAATTT‑AGACAAAAGGAGGAAAAATAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAt                                                                   <  1:69382/145‑1 (MQ=255)                                                                                         ttttGTGAAATGAGGGATTTCACATTTTTAATTA‑ATAACAAAGGAGAAAAAAAAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                                                  <  1:346713/145‑1 (MQ=255)                                                                                                 aggcaCGGAGTCCATTTTTTTAATTA‑AGATAAAAGGAAAAAGAAAAACAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc                                                    <  1:391445/147‑1 (MQ=255)                                                                                                 aaTGAGGGAGTTCATTTTTTAAATTA‑AGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                                      <  2:643945/149‑1 (MQ=255)                                                                                                 aaTGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc                                                    >  2:757084/1‑151 (MQ=255)                                                                                                       ggAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCAt                                              >  1:602652/1‑151 (MQ=255)                                                                                                           ttCACATTTTTAATTT‑AGATAAAAGGAGGAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt                                           <  2:707419/150‑1 (MQ=255)                                                                                                                 tttttAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                                     >  2:310469/1‑150 (MQ=255)                                                                                                                 tttttAATTT‑AGAAAAAAGGAGGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                        <  2:56190/147‑1 (MQ=255)                                                                                                                 tttttAATTT‑AAAGAAAAGAAGAAAGAATAAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCa                                         <  2:851591/146‑1 (MQ=255)                                                                                                                  tttaAATTA‑AGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCaa                                    <  2:234983/150‑1 (MQ=255)                                                                                                                   tttAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                   >  2:720265/1‑150 (MQ=255)                                                                                                                   tttAATTA‑AAAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                   <  1:574997/150‑1 (MQ=255)                                                                                                                      taagtaagaaaaagggagaaAAAAAAAGAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAg                                 <  1:466680/136‑1 (MQ=255)                                                                                                                        ttt‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                              >  2:139757/1‑150 (MQ=255)                                                                                                                        ttt‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                             >  2:32040/1‑151 (MQ=255)                                                                                                                            agaaaaagggaggaAAAAAAAGAAAAAAACGACCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                               <  1:719126/138‑1 (MQ=255)                                                                                                                            aaaaaaaaGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                          <  1:972589/146‑1 (MQ=255)                                                                                                                              aTCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                        >  2:895015/1‑151 (MQ=255)                                                                                                                                cAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                      >  2:280511/1‑151 (MQ=255)                                                                                                                                 aaaaaaaaaagaaaaagaaaAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                          <  2:720032/138‑1 (MQ=255)                                                                                                                                 aaaGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                      >  2:504234/1‑150 (MQ=255)                                                                                                                                 aaaGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                      <  1:381957/150‑1 (MQ=255)                                                                                                                                    ggagaaAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt                   <  2:123057/150‑1 (MQ=255)                                                                                                                                        gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg              >  2:849585/1‑151 (MQ=255)                                                                                                                                            aaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAgg          >  1:291546/1‑151 (MQ=255)                                                                                                                                               aaGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg       >  1:20199/1‑151 (MQ=255)                                                                                                                                                    gAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATg  >  2:894693/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       ATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATG  >  NC_000913/1435581‑1435877

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