Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,173	0	T	C	43.3%	89.8 / 14.1	60	I310V (ATT→GTT)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (28/6);  new base C (0/26);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.45e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACC  >  NC_000913/2921023‑2921323                                                                                                                                                         |                                                                                                                                                       agGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAaataat                                                                                                                                                        <  2:143421/151‑1 (MQ=255)   gCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcc                                                                                                                                                      <  1:21569/151‑1 (MQ=255)    ccAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATccc                                                                                                                                                     >  2:872117/1‑151 (MQ=255)         gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                                 >  1:672942/1‑150 (MQ=255)          ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccca                                                                                                                                                <  1:435010/150‑1 (MQ=255)             tAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                            >  2:175523/1‑151 (MQ=255)              aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcaca                                                                                                                                           <  2:383517/151‑1 (MQ=255)                aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                         >  2:127342/1‑151 (MQ=255)                  ccGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc                                                                                                                                       >  1:1007175/1‑151 (MQ=255)                    gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgccg                                                                                                                                     >  2:528072/1‑151 (MQ=255)                             ccACTCAACGTTGACGTAGTTGCAGAACCCCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAaaaccg                                                                                                                            <  2:390964/151‑1 (MQ=255)                             ccACTCAACGTTGACGTAGTTGCAGAACCCCACCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAATCCCCCCCAGCACACCGCCGATAaaac                                                                                                                              <  2:795898/149‑1 (MQ=255)                              cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                           >  2:744385/1‑151 (MQ=255)                                 tCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaa                                                                                                                         >  2:863860/1‑150 (MQ=255)                                     cGTTGACGTAGTTGCAGACCACCACCACCATCGACAGCAACATGCCCATTACGATCGTTTTTTT‑ACGCGCAATGTTCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAATCCCACCCAGCCCACCGCCGATAaaaccgcaaaccgc                                                                                                                    <  1:111156/151‑1 (MQ=255)                                                    aGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                                     >  1:768892/1‑151 (MQ=255)                                                      aaCCCCACCCCCATCGACAGAAACATCCCCATTCCGATGGGTTTTTT‑CCGCGAAATGTCCAG‑CGACCCCGTGCGGCGCACCAGCCGGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca                                                                                                    <  1:168502/150‑1 (MQ=255)                                                      aaCCCCACCCCCACCGCCAGAACCATGCCCATTACGACCGGTGTTTT‑CCGCGAAATGTTAAG‑CGCCCCCGGGGGGGGCATCACCCGGCGGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                                     <  1:967661/149‑1 (MQ=255)                                                        caccatcaccatCGACAGCAACATCCCCATTACGAGCGGTGTTTT‑ACGCGCAATGTTCAG‑CGACCCCGTGCGGCGCATCAGCCAGCCGAAAATAACCCCCCCAAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaaa                                                                                                  <  2:624170/150‑1 (MQ=255)                                                            accaccatcGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTCCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                              <  2:48586/148‑1 (MQ=255)                                                                ccatcGACAGCAACATGCCCTTTACGATCGGTTTTTT‑CCGCGAAATGTTCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                              <  1:489298/146‑1 (MQ=255)                                                                 catcGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                        >  1:828707/1‑151 (MQ=255)                                                                    cGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                        >  2:600888/1‑148 (MQ=255)                                                                     gACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                                     >  2:286331/1‑150 (MQ=255)                                                                     gACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                                     >  2:929867/1‑150 (MQ=255)                                                                      aCAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTCCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                                     <  2:236854/149‑1 (MQ=255)                                                                           aaCATGCCCATGACGATCGGTGTTTT‑CCGCGAAATGTCCAG‑CACCCCCGGGCGGCGCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAAt                                                                                   <  1:850133/146‑1 (MQ=255)                                                                              aTGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATcc                                                                           >  1:330133/1‑151 (MQ=255)                                                                                ccccttgcgtaTCGGTTTTTC‑ACGCGAACTGTTCAG‑CGCCCCTGCGCGGCGCCTCAGCCGCTCAAAAACATCCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                          <  1:903706/140‑1 (MQ=255)                                                                                 cccATTACGATGGGTTTTTT‑CCGCGCATTGTCCAG‑CGACCCCGTGCGGCGCATCAGCCGGCCGAAAATAACCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                         <  1:857905/150‑1 (MQ=255)                                                                                 cccATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                        >  1:109457/1‑151 (MQ=255)                                                                                 cccATTACGATCGGTGTTTT‑ACGCGAAATTTCCAG‑CGACCCCGTGCGGCCCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                          <  2:352032/149‑1 (MQ=255)                                                                                   cATTTCGGTCGGTTTTTTAACGCG‑AATTTTCCGCCGACCCCGGGCGCCGCATCAGCCATTCGGAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                       <  1:854787/151‑1 (MQ=255)                                                                                        cGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTg                                                                 >  1:923066/1‑151 (MQ=255)                                                                                          aTCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                                >  1:198955/1‑150 (MQ=255)                                                                                                  ttt‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                        >  1:788010/1‑150 (MQ=255)                                                                                                  ttt‑ACGCGCAATGTCCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                        <  1:601355/150‑1 (MQ=255)                                                                                                        gcgcAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                                  >  2:824832/1‑151 (MQ=255)                                                                                                        gcgcAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                                  >  1:397592/1‑151 (MQ=255)                                                                                                          gcAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                 >  2:825275/1‑150 (MQ=255)                                                                                                          gcAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                 >  2:164697/1‑150 (MQ=255)                                                                                                            aaTGTCCAG‑CGACCCCGTGCGGCCCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                    <  1:793/145‑1 (MQ=255)                                                                                                            aaTGTCCAG‑CGACCCCGCCCGCCGCACCAGCCAGCCGGAAATAACCCCACCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                                  <  1:106406/147‑1 (MQ=255)                                                                                                                  cAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                         >  1:594617/1‑150 (MQ=255)                                                                                                                      cGACCCCGTGCGGCGCACCAGCCAGTCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                            <  1:929688/144‑1 (MQ=255)                                                                                                                       gaccccGTGCGGCGCATCACCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGta                                        <  2:492825/144‑1 (MQ=255)                                                                                                                       gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                                     >  2:14215/1‑150 (MQ=255)                                                                                                                          cccggcgggcgCATCACCCAGTGGAAAATAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATaaa                                   <  1:425823/142‑1 (MQ=255)                                                                                                                               gcggccCATCAGCCAGTCGGAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                              <  1:947260/149‑1 (MQ=255)                                                                                                                                cggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                             >  1:3928/1‑149 (MQ=255)                                                                                                                                cggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                             >  2:709362/1‑149 (MQ=255)                                                                                                                                 ggcgCATCAGCCAGTCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaag                           <  2:873427/150‑1 (MQ=255)                                                                                                                                 gccccATCACCCAGTCGAAAAAACCCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                             <  1:863780/144‑1 (MQ=255)                                                                                                                                     cATCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaa                            <  1:986645/145‑1 (MQ=255)                                                                                                                                           caagcggaaaaaaaCCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTg                <  1:711153/143‑1 (MQ=255)                                                                                                                                            cAGGCGAAAATACCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                      <  1:970431/144‑1 (MQ=255)                                                                                                                                             aGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt              >  1:379304/1‑151 (MQ=255)                                                                                                                                              gccgaaaaaaaCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATa             <  1:872034/146‑1 (MQ=255)                                                                                                                                                ggaaaaTAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa           <  2:132112/148‑1 (MQ=255)                                                                                                                                                cGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc            >  1:376305/1‑150 (MQ=255)                                                                                                                                                  aaaaTAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt         <  2:462711/150‑1 (MQ=255)                                                                                                                                                        atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc   >  2:769953/1‑151 (MQ=255)                                                                                                                                                        atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc   >  2:548030/1‑151 (MQ=255)                                                                                                                                                        atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc   >  1:125283/1‑151 (MQ=255)                                                                                                                                                         tCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAcc  >  2:607631/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTT‑ACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACC  >  NC_000913/2921023‑2921323

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