Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	41.3%	92.1 / 35.0	64	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 (34/3);  new base A (0/26);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.02e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACG  >  NC_000913/1435579‑1435865                                                                                                                                                     |                                                                                                                                           gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                          <  2:157307/151‑1 (MQ=255)  aaTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                          >  2:793198/1‑150 (MQ=255)         aaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCg                                                                                                                                  <  1:488878/151‑1 (MQ=255)                aTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATaa                                                                                                                            >  2:762122/1‑150 (MQ=255)                   gggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGc                                                                                                                        >  2:125999/1‑151 (MQ=255)                     gACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCgg                                                                                                                      >  1:1042871/1‑151 (MQ=255)                        aaGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                    >  2:412831/1‑150 (MQ=255)                         aaaaTCTTTCTTTTAATAAACCCCCTTTTTTTAGTACCCGAAAAAAGTGAAACCCAAGATGACCTTTTTGGTGTAATGAGGAAGTACACTTTTTTAATTAAGATAAAAGGGGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCag                                                                                                                  <  1:550998/149‑1 (MQ=255)                            atgtttCTTTTAATAAACCCTCTTTTTTTAGTATCCGCAAAAAGTGAAACGCCAGAGAACATTTTTTGTGAAATGAGGGATTTCACATTTTTAATTAAGATAAAGGGAGAAAAAATAAGCAAAAAAACCCCGCCATAACAGCGGGCaggag                                                                                                               <  2:497989/147‑1 (MQ=255)                                 tgtGTTTAGAAAGCCTCGTGTTTTACTTACCCCATAAGGTGAAAAACCAGATGACGCTTTTGTCGAAATGCCGCGGTTCACATTTTAATTTTATACAAAAGGGGAAAGATTAACGAAAAAAAGCCCGCCATAACAGCGGGCaggag                                                                                                               <  1:137948/144‑1 (MQ=255)                                  cTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                          >  2:491682/1‑150 (MQ=255)                                     ttaaaaaacCCTCTTTTTTTAGTATCCGAAAAAAGTGTACCGCCAGATGACTTTTTTTGTGTAATGGGGGAGTTCACATTTTTAATTAAGATCAAAGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                        <  1:338909/144‑1 (MQ=255)                                         tAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                   >  1:84558/1‑150 (MQ=255)                                         tAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                   >  2:440090/1‑150 (MQ=255)                                         aaaacCCTCTTTTTGTAGTACCCAAATAAAGTATAACCCCAGATGACTTTTTTTGTGAAAGGCGGGATTACACTTTTTAAATTTAGAAAAAAGGGAAAAAAAAAAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                  <  2:191201/150‑1 (MQ=255)                                          aaaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                 >  1:242775/1‑151 (MQ=255)                                           aaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                 >  2:891702/1‑150 (MQ=255)                                              cctcctGTTTTAATTTCCGAAAAAAGTGCAACACCAGACGACTTTTTTTGTTTAATAAGGGAGTTAACTTTTTAAATTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                                                              <  1:58729/150‑1 (MQ=255)                                                tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                            >  1:904488/1‑150 (MQ=255)                                                      tttAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                      >  1:1013437/1‑150 (MQ=255)                                                           tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                >  1:709236/1‑151 (MQ=255)                                                                  aTAAAGTGTAACGCCAGATGACTCTTTTTGTGTAATGAGGGAGTTCACATTTTAAATTTAGAACAAAGGAGGAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                           <  1:346576/149‑1 (MQ=255)                                                                  aTAAAGTGTAACGCCAGATGACATTTTTTGTATAATGACGGATTTCTCATTTTTAATTTAGAACAAGGGAGGAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                           <  2:773949/149‑1 (MQ=255)                                                                  aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                         >  2:732155/1‑151 (MQ=255)                                                                  aTAAAGTAAACCGCAAGATGACTTTTTTTGTTTAAGGACGGAGTTCATATTTTTAATTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAata                                                                          <  2:98015/150‑1 (MQ=255)                                                                     aaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcat                                                                      >  1:639628/1‑151 (MQ=255)                                                                       gtgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                    >  2:344647/1‑151 (MQ=255)                                                                        tgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg                                                                   >  1:312664/1‑151 (MQ=255)                                                                             cccaaGATGACTTTTTTTGTGAAAGGAGGGATTTCACTTTTTAAATTTAGAAAAAAGGAGAAAAAATAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTaga                                                               <  2:632284/146‑1 (MQ=255)                                                                             cGCCAAATGACTTTTTTTGTATAAGGACGGAGTTCACATTTTTAATTAAAAAAAAGGGAGGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                                                              <  2:898690/151‑1 (MQ=255)                                                                                   aTGCCATTTTTTGTGTAATGACGGATTTCACATTTTTAATAGAGATCAAGGAAGAAAAAAAAAGCACAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc                                                         <  1:614366/150‑1 (MQ=255)                                                                                    taccatttttttTGTAATGACGGAGTTCTTATTTTTATTGTAGAACAAAGGGGAAAGAATAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                           <  1:611001/141‑1 (MQ=255)                                                                                    tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAt                                                       >  1:602974/1‑151 (MQ=255)                                                                                      acatTTTTTGTGTAATGAGGGAGTTCACTTTTTTAATTAAGAAAAAGGGGAAAAAAATAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGAt                                                            <  2:367702/144‑1 (MQ=255)                                                                                         cTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                   >  2:188211/1‑150 (MQ=255)                                                                                          tttttGTGAAATGGGGGATTTCACTTTTTTAATTTAAAAAAAGGGAGAAAAAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                   <  1:240384/149‑1 (MQ=255)                                                                                            tttGTTAAATGACGGATTTCACTTTTTTATTTTAGAAAAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                <  1:209970/150‑1 (MQ=255)                                                                                            tttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                >  1:531107/1‑150 (MQ=255)                                                                                            tttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCg                                               >  2:969913/1‑151 (MQ=255)                                                                                             tggtgtAATGAGGGAGTCCATATTTTTATTTAAAATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCg                                               <  1:628229/148‑1 (MQ=255)                                                                                                 gtAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                           >  2:766924/1‑150 (MQ=255)                                                                                                   aaTGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc                                        >  1:473746/1‑151 (MQ=255)                                                                                                        cGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTc                                    >  2:988466/1‑150 (MQ=255)                                                                                                         ggATTTCATTTTTTAAATTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCa                                   <  1:125766/150‑1 (MQ=255)                                                                                                           aGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                >  1:1027991/1‑151 (MQ=255)                                                                                                             ttcaattttttATTTTAAAAAAAAAAAGAAAAAAAAAACAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt                               <  2:541427/144‑1 (MQ=255)                                                                                                             ttCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                              >  2:735577/1‑151 (MQ=255)                                                                                                                acaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATc                           >  2:417110/1‑151 (MQ=255)                                                                                                                  aTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATcc                          >  1:602909/1‑150 (MQ=255)                                                                                                                  aTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATc                           >  1:89573/1‑149 (MQ=255)                                                                                                                   tttttAATTAAGATAAAAGGAGGAAAAATAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                              <  1:283690/145‑1 (MQ=255)                                                                                                                     tttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                       >  2:1121700/1‑150 (MQ=255)                                                                                                                      ttaattaaaaaaaaaGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACa                      <  2:1042745/138‑1 (MQ=255)                                                                                                                        aaTTTAGAAAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAg                     <  1:893758/149‑1 (MQ=255)                                                                                                                          tttAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                 >  2:1012891/1‑151 (MQ=255)                                                                                                                          tttAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                  <  1:500789/150‑1 (MQ=255)                                                                                                                            aaaaaaaagggaggaAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                 <  1:1055253/140‑1 (MQ=255)                                                                                                                             agacaaaaGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGa                <  2:904374/144‑1 (MQ=255)                                                                                                                             aaataaaaGGAGAAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt              <  1:26366/146‑1 (MQ=255)                                                                                                                              aataaagggagaaAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt             <  1:493025/144‑1 (MQ=255)                                                                                                                                      gaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc      >  1:312780/1‑150 (MQ=255)                                                                                                                                       aggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCa     >  1:694920/1‑150 (MQ=255)                                                                                                                                       agaaagaaTACGAAGAAAAATCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCa     <  1:417879/147‑1 (MQ=255)                                                                                                                                         gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTACAGTCAACg  >  2:255126/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                           GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACG  >  NC_000913/1435579‑1435865

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