Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,089,852	0	A	G	63.4%	85.7 / 80.3	82	intergenic (‑139/‑144)	yefM/hisL	antitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base G (29/23);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAG  >  NC_000913/2089714‑2089987                                                                                                                                           |                                                                                                                                        gTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGACCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGa                                                                                                                                         <  2:37692/139‑1 (MQ=255)    aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAg                                                                                                                                      <  2:162251/139‑1 (MQ=255)    aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                      <  2:571087/139‑1 (MQ=255)          ctctTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCctct                                                                                                                                >  2:544840/1‑139 (MQ=255)           tctTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCa                                                                                                                                       <  1:269259/131‑1 (MQ=255)            ctTTTGTACAGTTAATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTa                                                                                                                              <  1:16547/139‑1 (MQ=255)                  tACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                                                        >  1:304110/1‑139 (MQ=255)                   aCAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcct                                                                                                                       <  1:9317/139‑1 (MQ=255)                       ttCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCa                                                                                                                                       <  1:114433/119‑1 (MQ=255)                               aCAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCttt                                                                                                                  >  1:202665/1‑132 (MQ=255)                               aCAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCttt                                                                                                                  <  2:202665/132‑1 (MQ=255)                                 aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtca                                                                                                         >  2:361079/1‑139 (MQ=255)                                     atgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                     >  2:280581/1‑139 (MQ=255)                                      tgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattcat                                                                                                    >  2:358205/1‑139 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGc                                                                                                                            >  2:405842/1‑113 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGc                                                                                                                            <  1:405842/113‑1 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtca                                                                                                         >  1:547604/1‑132 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtca                                                                                                         <  2:547604/132‑1 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTa                                                                                                  >  1:45638/1‑139 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTa                                                                                                  >  1:121550/1‑139 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTa                                                                                                  >  2:69072/1‑139 (MQ=255)                                        aGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTa                                                                                                  >  2:535491/1‑139 (MQ=255)                                          cGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTa                                                                                                                              <  2:242468/109‑1 (MQ=255)                                          cGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTa                                                                                                                              >  1:242468/1‑109 (MQ=255)                                            ttaattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACa                                                                                              <  1:205552/139‑1 (MQ=255)                                             taattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                                                     <  2:358358/99‑1 (MQ=255)                                             taattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                                                     >  1:358358/1‑99 (MQ=255)                                               attaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAt                                                                                           >  1:204583/1‑139 (MQ=255)                                               attaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAt                                                                                           >  2:501315/1‑139 (MQ=255)                                                  aaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCa                                                                                        <  2:474350/139‑1 (MQ=255)                                                  aaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCa                                                                                        <  2:178714/139‑1 (MQ=255)                                                   aCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAt                                                                                       <  2:121550/139‑1 (MQ=255)                                                   aCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAt                                                                                       <  1:358205/139‑1 (MQ=255)                                                      atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGc                                                                                    >  2:41746/1‑139 (MQ=255)                                                       tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTCAGAGGGCAGTCATCTTAGCCCCCTACTTTCCCCGCTCATTCATTAAACAAATCCATTgac                                                                                   >  1:195942/1‑137 (MQ=255)                                                       tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTaaa                                                                                                >  2:405639/1‑126 (MQ=255)                                                       tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTaaa                                                                                                <  1:405639/126‑1 (MQ=255)                                                         tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCAt                                                                                 >  2:351776/1‑139 (MQ=255)                                                          attaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGt                                                                                                                                     <  2:132829/86‑1 (MQ=255)                                                          attaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGt                                                                                                                                     >  1:132829/1‑86 (MQ=255)                                                          attaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                                                                <  2:304110/139‑1 (MQ=255)                                                            taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                                                              <  1:501315/139‑1 (MQ=255)                                                               ttaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                                                     >  2:54348/1‑81 (MQ=255)                                                               ttaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                                                     <  1:54348/81‑1 (MQ=255)                                                                 aGTTTGTAGATCAAGGTCTTGTCATTGAGTCGAAAATCAAGCCTTCGCTACTCTTGGTCGCATCAGCTAAGATGGCTGTCCACTTAGCCCCCTCCTTTCCCCGCTCATTCATTACACAAATCCATTGCCATAAAatata                                                                         >  1:63403/1‑139 (MQ=255)                                                                   tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAAtatata                                                                       >  1:281147/1‑139 (MQ=255)                                                                   tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAt                                                                                       >  1:296908/1‑123 (MQ=255)                                                                   tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAt                                                                                       <  2:296908/123‑1 (MQ=255)                                                                     tGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaa                                                                     >  2:17603/1‑139 (MQ=255)                                                                        aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaaaa                                                                  >  2:250397/1‑139 (MQ=255)                                                                            cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGccc                                                              >  1:608923/1‑139 (MQ=255)                                                                            cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGccc                                                              >  2:574079/1‑139 (MQ=255)                                                                            cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtc                                                                                                          <  2:420712/95‑1 (MQ=255)                                                                            cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtc                                                                                                          >  1:420712/1‑95 (MQ=255)                                                                                 tATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCt                                                         >  1:25786/1‑139 (MQ=255)                                                                                    tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                     >  2:295714/1‑92 (MQ=255)                                                                                    tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                     <  1:295714/92‑1 (MQ=255)                                                                                    tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATCAAATATATAAAAAAGCCCTTGCTTTc                                                      >  2:294166/1‑139 (MQ=255)                                                                                    tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                                                      >  1:152585/1‑139 (MQ=255)                                                                                     gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCt                                                     <  1:280581/139‑1 (MQ=255)                                                                                         gTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg                                                 >  1:214158/1‑139 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGa                                                                                                                                         >  2:201739/1‑50 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGa                                                                                                                                         <  1:201739/50‑1 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTcc                                                                                                                                   <  2:154511/56‑1 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTcc                                                                                                                                   >  1:154511/1‑56 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                                                                <  2:602760/107‑1 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                                                                >  1:602760/1‑107 (MQ=255)                                                                                          tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt                                                >  2:605162/1‑139 (MQ=255)                                                                                           gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCg                                                                                                             <  1:206212/77‑1 (MQ=255)                                                                                           gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCg                                                                                                             >  2:206212/1‑77 (MQ=255)                                                                                           gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                      >  1:73745/1‑52 (MQ=255)                                                                                           gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                      <  2:73745/52‑1 (MQ=255)                                                                                              aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGcccc                                                                                                                         <  2:21091/62‑1 (MQ=255)                                                                                              aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGcccc                                                                                                                         >  1:21091/1‑62 (MQ=255)                                                                                              aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaaa                                            <  2:214158/139‑1 (MQ=255)                                                                                              aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaaa                                            <  2:152585/139‑1 (MQ=255)                                                                                               cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAg                                           >  1:584756/1‑139 (MQ=255)                                                                                                gAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGt                                          >  2:141393/1‑139 (MQ=255)                                                                                                    aTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGtt                                      <  1:351776/139‑1 (MQ=255)                                                                                                         ggCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGt                                 <  1:17603/139‑1 (MQ=255)                                                                                                         ggCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGt                                 <  2:90656/139‑1 (MQ=255)                                                                                                               gCTATTTTTGGTGCCATCAGCTAAGAGGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAg                           >  1:342027/1‑139 (MQ=255)                                                                                                                     tttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCa                     <  1:535491/139‑1 (MQ=255)                                                                                                                                       gagGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTcaca   <  1:428643/139‑1 (MQ=255)                                                                                                                                        agGGCAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACag  <  1:237000/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        GTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAG  >  NC_000913/2089714‑2089987

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