Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	27.1%	99.9 / 17.2	49	intergenic (‑108/+33)	uspF/ompN	universal stress protein F/outer membrane porin N
Reads supporting (aligned to +/- strand):  ref base G (33/2);  new base A (0/13);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 5.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGCGTGCAGGT  >  NC_000913/1435485‑1435977                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            gggCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGc‑c                                                                                                                                                                                                                                                    <  2:29531/251‑1 (MQ=255)   gCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGAAGAAAAAGc‑ccc                                                                                                                                                                                                                                                  >  2:213408/1‑250 (MQ=255)   gCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCGCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAAGTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCg                                                                                                                                                                                                                                                  <  2:532279/251‑1 (MQ=255)        gAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCAt                                                                                                                                                                                                                                              >  2:586056/1‑250 (MQ=255)        gAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCAt                                                                                                                                                                                                                                              >  2:543106/1‑250 (MQ=255)            gTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAAc                                                                                                                                                                                                                                           >  2:629088/1‑249 (MQ=255)              aTTACCGTAAGAAAATGAACCTCTGCATCATCAATTTTTGCCTCTTCCTCAACGTGTCTAATCACGCGTTGAGTTAATTCTGAATCGAAAATATCGAGAGGGACAAGAATCGTTCTTTTCAAAAACCCTCCTGTTTTATTATCCGCATAAAGTGTACCGCAAGAGGACATTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGAAAAAAGGAGAAAAAATAAGAAAAAAAAGC‑CCGCCATAACAGc                                                                                                                                                                                                                                        <  1:461839/250‑1 (MQ=255)                        ggAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCaggagg                                                                                                                                                                                                                              >  1:80390/1‑250 (MQ=255)                        ggAAATGAACATCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAACAGC‑CCGCCATAACAGCGGGCaggagg                                                                                                                                                                                                                              >  2:235602/1‑250 (MQ=255)                                aCCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                                                                                                                                                                     >  2:610018/1‑251 (MQ=255)                                              aaTCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                                                                                                                       >  2:299847/1‑251 (MQ=255)                                               aTCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                                                                                                                      >  1:221035/1‑251 (MQ=255)                                                                    tGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                                                                                                                                  >  1:442797/1‑250 (MQ=255)                                                                              cgcgTTGAGTTAATTCTGAATCGGAAATATCGCTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTACAGCAACAATATCATCGGTAGAGATGCCAt                                                                                                                                                                       >  2:383784/1‑251 (MQ=255)                                                                                   ggatttatttctaaagagaaaaTATCGAGAGGAAAAAGAATCTTTTTTTTCAAAAAACCTCTTTTTTTATTACCCGAAAAAAGTGTAACGCAAGAGAACATTTTTTGTGTAAGGAGGGATTT‑CATTTTTTAAATTTAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                                                                                                                                   <  2:250338/231‑1 (MQ=255)                                                                                    gatttAATTCTAACTCGAAAATACCGATAGGGACAAGAATCTTTCTTTTCATAAAACCTCTTTTTTTAGTCTCCTCATAAAGTGTAACGCGAGAGGACATTTTTTGTATGTGGACG‑AGTTTCATATTTTTAATTTAGAGAAAAGGAGGAAGAAAAAGCAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCa                                                                                                                                                                  <  1:206790/247‑1 (MQ=255)                                                                                     atttatttttGAATCGGAAATATCGAAAGGAAAAAGAATCGTTCTTTTCAAAAACCCTCCTTTTTTAGTATCCGAAAAAAGTGAAACGCCAGATGACTTTTTTTGTTAAATGAGGGAGTT‑CACTTTTTTATTTAAGAAAAAAGGAGGAAAAATAAGCAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                                                                                                                                <  1:537234/242‑1 (MQ=255)                                                                                                aaTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGt                                                                                                                                                     >  1:447352/1‑251 (MQ=255)                                                                                                          atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATc                                                                                                                                           >  1:197585/1‑251 (MQ=255)                                                                                                              aTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                                                       >  1:321782/1‑251 (MQ=255)                                                                                                                    aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                                                                                                                                           >  1:461729/1‑209 (MQ=255)                                                                                                                             gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                                                                                                                          >  1:429098/1‑249 (MQ=255)                                                                                                                                        aaacccccttGTTTTATTATCCGCATAAAGTGTACCCCCAGATGCCTTTTTTGGTGTAATGACGGATTT‑CATATTTTTAATTTAAAAAAAAGGAGGAAAAAAAAGAAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGt                                                                                                             <  1:511569/242‑1 (MQ=255)                                                                                                                                         aacccttcttttttAGTACCCGAATAAAGTGTAACACCAGAGGACTTTTTTTGTGAAATGACGGAGTT‑CACTTTTTTAATTAAGATAAAAGGAAAAAAAATAAGAAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTg                                                                                                            <  2:92031/241‑1 (MQ=255)                                                                                                                                             ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGAc                                                                                                         >  1:406176/1‑250 (MQ=255)                                                                                                                                               ccttttttATTACCCGCATAAAGTGTAACGCCAGATGACTTTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATAAAAGGAGGAAAAATAAGAAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATg                                                                                                      <  1:586056/247‑1 (MQ=255)                                                                                                                                                  gTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTCGTGTAATGACGGAGTT‑CACATAAAGAAAGCAGATCAAAGGAGGAAGAAAAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATgatt                                                                                                   >  1:392831/1‑248 (MQ=255)                                                                                                                                                     taagtacccaaaaaaaGTGAAACGCCAGATGCCATTTTTGGTGTAAGGGCGGAGTT‑CATATTTTTAATTAAGAGAAAAGGAGAAAAAATAAGAAAAAAAACC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGttttt                                                                                                 <  1:611383/237‑1 (MQ=255)                                                                                                                                                      tAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTAt                                                                                               >  1:582711/1‑251 (MQ=255)                                                                                                                                                             tgCATAAAGTGTCACGCCAGCTGACACTTTTTGTGTAATGACGGAGTT‑AACATTTTTAATTTAGATCAAAGAAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCAAACAGGATGATTTTATAGTCAACGTCGGTGGACATGTTTTTATTGAAAt                                                                                         >  2:26392/2‑250 (MQ=255)                                                                                                                                                                  aatgtgtAACGCCAGATAACTTTTTTTTTTAAATGACGGATTT‑CAAATTTTAAATTTAAATAAAAGGAGGAAAAAAAAGAAAAAAACGA‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAg                                                                                   <  2:123710/248‑1 (MQ=255)                                                                                                                                                                               agttGACACTTTTTGTGTAAAGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCGGTCATCTTCATCAATCAGGTTGATTTTATAGTCAACGTAGGGGGACATGTTTTTATTGAAATAGTCATTCGCGCCACTATc                                                                      >  2:87709/4‑251 (MQ=255)                                                                                                                                                                               aGATGACTTTTTTTGTGTAATGACGGATTT‑CACATTTTTAATTGAGAACAAAGGAGGAAAAAAAAGAAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATc                                                                      <  1:634272/251‑1 (MQ=255)                                                                                                                                                                                               gtAATGACGGATTG‑CAAATTTTTAATTTAAAAAAAAGGAGAAAGAAAAAGAAAAAAAACC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAg                                                       <  1:139698/250‑1 (MQ=255)                                                                                                                                                                                                 aaTGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAAcc                                                         >  1:38918/1‑246 (MQ=255)                                                                                                                                                                                                      cGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTat                                               >  1:369638/1‑251 (MQ=255)                                                                                                                                                                                                               acaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAac                                        >  2:279039/1‑250 (MQ=255)                                                                                                                                                                                                                   ttttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGt                                                                                                             >  2:81896/1‑177 (MQ=255)                                                                                                                                                                                                                   tttaaaTTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAAACC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGt                                                                                                             <  1:81896/173‑1 (MQ=255)                                                                                                                                                                                                                                  aaGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGcac                    >  2:395723/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGcac                    >  1:82793/1‑251 (MQ=255)                                                                                                                                                                                                                                     taggaTGAATAAACAGAAAAAGC‑CCGCCATAACAGCGGGCCGGAGGATATAGAAATGATAAAACAGACCTAACGCCACAATATCATCGGTAGAGCGCCCATTGGAAGCGTAGAAGTTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTCTTGATCTATTAAGTCTCGCCACTATGGGACTCTTTAACCAGATCTATATCATCAACACCTGCCGGGGCTTCTGcaccac                 >  2:46493/2‑251 (MQ=255)                                                                                                                                                                                                                                      agaaagaaAAAGCA‑AAAAAGCACCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCaccacc                <  1:618741/248‑1 (MQ=255)                                                                                                                                                                                                                                       ggaAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACAAGATCTTTATCATCAACACCTGCCGGGTTGTATGCACCAccc               >  2:350604/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaaaaaaGCAAAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGc             <  1:629088/245‑1 (MQ=255)                                                                                                                                                                                                                                            aaTAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGc          >  2:143673/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGCg         >  2:386874/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGCg         >  1:289506/1‑251 (MQ=255)                                                                                                                                                                                                                                                aGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGCTGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCAGGTTGTCTGCACCACCCGCAGCGTGc      >  1:566385/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGCGTGCAGGt  >  1:249758/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            GGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTT‑CACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGC‑CCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGCGTGCAGGT  >  NC_000913/1435485‑1435977

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