Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	43.5%	92.4 / 27.4	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 (31/4);  new base A (0/27);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTG  >  NC_000913/1435579‑1435871                                                                                                                                                     |                                                                                                                                                 gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                                <  2:643113/151‑1 (MQ=255)     cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAg                                                                                                                                            <  1:482222/151‑1 (MQ=255)      ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                                           >  2:242050/1‑151 (MQ=255)       gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGcc                                                                                                                                          >  2:337502/1‑151 (MQ=255)        aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                         <  1:450072/151‑1 (MQ=255)            atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                      >  2:134502/1‑150 (MQ=255)            atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                      >  1:510255/1‑150 (MQ=255)                 tAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAc                                                                                                                                 >  1:825160/1‑150 (MQ=255)                    ggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCTATAACAGCg                                                                                                                             <  1:639122/151‑1 (MQ=255)                         aGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                          >  2:586923/1‑149 (MQ=255)                                ttCTGTTAAAAAAACCTCCTTTTTTAGTATCCGAAAAAAGTGTAACGCAAGATGACATTTTTTGTGTAAGGAGGGATTTCACATTTTTAATTTAGATAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCaggag                                                                                                                     <  2:392081/147‑1 (MQ=255)                                     ttaaTAAACCCCCCTTTTTTATTATCCGAAAAAAGTGAAACGCAAGATGACTTTTTTTGTTAAATGAGGGAGTTCACTTTTTTAATTTAGATAAAAGGAGAAAGAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                <  2:793359/144‑1 (MQ=255)                                        aTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                         >  2:743800/1‑151 (MQ=255)                                         tAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                         >  2:338216/1‑150 (MQ=255)                                          aaaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                       >  2:49489/1‑151 (MQ=255)                                            aaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                                                      >  1:55737/1‑150 (MQ=255)                                                 cctTTTTTAGTACCCGCATAAAGTGTACCGCCAGATGACATTTTTTGTGTAAGGACGGAGTTCACATTTTAAATTTAGATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                                                                                 <  2:313306/150‑1 (MQ=255)                                                   ttttttAGTACCCGCAAAAAGTGTACCCCCAAATAACATTTTTTGTGTAATGACGGATTTCACATTTTTAATTTAGATAAAGGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGa                                                                                               <  2:190369/148‑1 (MQ=255)                                                   tGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                                                              >  1:306551/1‑151 (MQ=255)                                                   tGTTTTAGTATCCGAATAAAGTGTAACGCCAGATAACTCTTTTTGTGAAAAGACGGAGTTCACATTTTTAATTTAAATAAAAGGAGGAAAAAAAAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                                                                    <  1:451562/145‑1 (MQ=255)                                                         atttcccGCATAAAGTGAAACGCAAGATGACATTTTTTGTGTAATGACGGAGTTCACTTTTTTAATTTAGATCAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                                                                         <  1:774321/145‑1 (MQ=255)                                                          gTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                       >  2:639038/1‑151 (MQ=255)                                                           tATCCGAATAAAGTGTAACGCCAGAGGACTTTTTTTGTGTAATGAGGGAGTCCACATTTTTATTTTAGATAAAAGGAGGAAGAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                       <  2:744023/150‑1 (MQ=255)                                                           tATCCGAAAAAGGTGTACCGCAAGATGACTTTTTTTGTGTAATGACGGATTTCACTTTTTTAATTTAGATAAAAGGAGAAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                       <  2:475333/150‑1 (MQ=255)                                                            aTCCGCATAAAGTGTAACGACAGATGACTTTTTTTGTGAAAGGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                       <  1:140622/149‑1 (MQ=255)                                                                gaaaaaaGTGTAACGCCAGATGACTTTTTTTGTGTAATGACGGAGTTCACTTTTTTAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                      <  1:55061/142‑1 (MQ=255)                                                                  aTAAAGTGTACCCCCAGATGACTTTTTTTGTGTAATGCCGAAGTTCAAATTTTTAATTAAGAACAAAGGAGAAAAAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                               <  2:679313/151‑1 (MQ=255)                                                                    aaaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatc                                                                              >  1:653634/1‑150 (MQ=255)                                                                       gtgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCATAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                          >  1:514604/1‑151 (MQ=255)                                                                       gtgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                          >  2:501484/1‑151 (MQ=255)                                                                         gtAACGCCAGATGACATTTTTTGTGTAATGACGGAGTTCACATTTTTAATTAAGATAAAGGGAGGAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg                                                                         <  2:926912/150‑1 (MQ=255)                                                                           aaCGCGAGATGACTTTTTTTGTGAAATGACGGATTTCACATTTTAAATTAAGATCAAGGGAGAAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                                                       <  1:683777/150‑1 (MQ=255)                                                                           aaCGCCAGATGACTTTTTTTGTGTAATGACGGATTTCACTTTTTTAATTAAGATAAAGGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                      <  1:244407/151‑1 (MQ=255)                                                                            aCGCCAGATGACATTTTTTGTGTAGTGACGGAGTTCACTTTTTTAATTTAGATCAAAGGGGAAAAAAAAGGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcat                                                                            <  1:686073/144‑1 (MQ=255)                                                                            aCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                        >  2:256070/1‑148 (MQ=255)                                                                                       cacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                          >  1:283069/1‑151 (MQ=255)                                                                                         cTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                         >  2:781937/1‑150 (MQ=255)                                                                                           ttttGTGTAATGACGGATTTCACATTTTTAATTTAGATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAg                                                       <  2:896678/150‑1 (MQ=255)                                                                                             ttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                      >  1:291161/1‑149 (MQ=255)                                                                                             ttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                    >  2:695672/1‑151 (MQ=255)                                                                                               gtatAATGACGGAGTTAACTTTTTTAATTAAAATCAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                    <  1:927851/146‑1 (MQ=255)                                                                                                   aaTGAGGGATTTCAATTTTTTAATTTAAACAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc                                              <  2:121835/151‑1 (MQ=255)                                                                                                   aaTGACGGAGTTCATTTTTTTAATTAAGAACAAAGGAGAAAAAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                                <  1:446583/149‑1 (MQ=255)                                                                                                   aaTGACGGAGTTCACATTTTTAATTTAGATAAAAGGAGAAAGAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg                                               <  2:910172/150‑1 (MQ=255)                                                                                                   aaTGACGGAGTTCACATTTTTAATTTAGACCAAAGGAGGAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                   <  2:420047/146‑1 (MQ=255)                                                                                                   aaTGACGGAGTTCACATTTTTAATTAAGAGAAAGGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                   <  2:447665/146‑1 (MQ=255)                                                                                                    atgaggggtttaaaaTTTTAAATTAAAAAAAAGGGAGGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                   <  2:460852/131‑1 (MQ=255)                                                                                                      gACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTg                                            >  1:42572/1‑150 (MQ=255)                                                                                                      gACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGt                                           >  1:712267/1‑151 (MQ=255)                                                                                                        cGGAGTTCACATTTTTAATTAAGAACAAAGGAGGAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                                <  1:615965/144‑1 (MQ=255)                                                                                                         ggATTTCACATTTTTAATTAAGAAAAAGGGAGGAAAAATAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCt                                             <  1:557904/146‑1 (MQ=255)                                                                                                               cacaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                  >  2:109432/1‑151 (MQ=255)                                                                                                                        aaTTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                         >  2:765079/1‑151 (MQ=255)                                                                                                                          ttaagaaaaaaGGAGAAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                        <  1:311344/142‑1 (MQ=255)                                                                                                                                     ggaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCTACAGGTTGATTTTATAGTc            >  2:603633/1‑151 (MQ=255)                                                                                                                                     ggagaaAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAg              <  1:907964/149‑1 (MQ=255)                                                                                                                                      gaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc            >  1:471496/1‑150 (MQ=255)                                                                                                                                      gaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc            >  1:730370/1‑150 (MQ=255)                                                                                                                                         gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg        >  2:800196/1‑151 (MQ=255)                                                                                                                                           agaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGt       >  1:574362/1‑150 (MQ=255)                                                                                                                                               tAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTg  >  2:555558/1‑151 (MQ=255)                                                                                                                                               tAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTg  >  1:874477/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTG  >  NC_000913/1435579‑1435871

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: ‑