Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,173	0	T	C	50.7%	62.2 / 104.4	69	I310V (ATT→GTT)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (31/3);  new base C (0/35);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-11

AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC  >  NC_000913/2921023‑2921322                                                                                                                                                        |                                                                                                                                                        agGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAaataat                                                                                                                                                         >  2:1137552/1‑151 (MQ=255)    ccAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑ccc                                                                                                                                                     >  2:597206/1‑151 (MQ=255)     cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑cc                                                                                                                                                      >  1:703385/1‑149 (MQ=255)      aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑cccac                                                                                                                                                   >  1:725654/1‑151 (MQ=255)       gcgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑cccacc                                                                                                                                                  >  2:961687/1‑151 (MQ=255)           cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGc                                                                                                                                              >  1:724819/1‑151 (MQ=255)               aaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCacac                                                                                                                                          <  1:496138/151‑1 (MQ=255)                 gCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACAccg                                                                                                                                        >  1:442653/1‑151 (MQ=255)                     atgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACAccgccg                                                                                                                                     >  2:809033/1‑150 (MQ=255)                              cccTCAACGTTGAGGAGGTTGCAGAACCCCATCCCCACCGACACCACCATCCCCATTACAACGGGTTTTTTCCGCCAATTGTTCAG‑CCATCCCGTGCGGGGCATCCGCCAGTGGGAAATAAT‑CCCACCCCGCACACCGCCGATAaaaccgc                                                                                                                           <  1:1015040/149‑1 (MQ=255)                                   accGTTGACGTAGTTGCAAAACCCCATCCCCACCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATTTTCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAAC‑CCCCCCCAGCACACCGCCGATAaaaccgcaaac                                                                                                                       <  2:274882/148‑1 (MQ=255)                                   aaCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAaaaccgcaaacc                                                                                                                      >  1:460044/1‑151 (MQ=255)                                   aaCGTTGACGTAGTTGCAAAACCCCACCACCACCGACAGCAACATGCCCATTACGAGCGGTGTTTTACGCGAAATGTTCAG‑CGCCCCCGTGCGGCGCATCAGCCAGTCGAAAATAAC‑CCCCCCCAGCACACCGCCGATAaaaccgcaaacc                                                                                                                      <  1:543370/151‑1 (MQ=255)                                    aCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                                     >  1:1148286/1‑151 (MQ=255)                                       ttGACGTAGTTGCAGAACCCCACCCCCATCGACAGCAACATGCCCTTTACGATCGGTTTTTTACGCGAAATGTCCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAAccgcc                                                                                                                   <  1:983732/150‑1 (MQ=255)                                          acggtttggCAGAACCCCACCCCCTGCGAAAAAAACAGGCCAATTCCGAGCGGGTTTTTACGCAAAATGTCCAG‑CGACCCCGGGCGGCACATCGGCAAGTGGAAAATAAC‑CCCCCCAAGCAACCCGCCGATAAAACCGCAAACCGCcggaa                                                                                                               <  1:839551/143‑1 (MQ=255)                                           cGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaac                                                                                                              >  1:1256032/1‑151 (MQ=255)                                              aGTTGCAAAACACCACCCCCATCGACAGCAACATGCCCTTTACGTCCGGTGTTTTCCGCGCAATTTCCAG‑CGACCCCGGGCGGCGCATCACCCAGCCGAAAAAAAC‑CCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                             <  1:740555/149‑1 (MQ=255)                                                    aGAACCCCACCCCCACCGACAGCAACATGCCCATTACGATCGGTTTTTTCCGCGAAATGTTCAG‑CGACCCCGTGCGGCGCATCACCCAGTCGAAAATAAC‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc                                                                                                       <  2:215449/149‑1 (MQ=255)                                                      acccccatccccatcGACAGCACCATCCCCATTACGATCGGTGTTTTACGCGCAATTTTCAG‑CGCCCCCGTGCGGCCCATCCGCCAGTGGAAAATAAC‑CCCCCCCAGAACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGc                                                                                                        <  1:1271432/142‑1 (MQ=255)                                                      aaCCCCACCCCCATCGACAGCAACAGCCCCATTACGATGGGTGTTTTCCGCGAAATGTCCAC‑CGCCCCCGTGCGGCGCACCGCCCAGTCGAAAAAAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa                                                                                                   <  1:1156178/151‑1 (MQ=255)                                                      aaCACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca                                                                                                    >  1:540714/1‑150 (MQ=255)                                                       acacCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa                                                                                                   >  1:1041723/1‑150 (MQ=255)                                                          ccatccccaccGACAGCAACATCCCCATTACGATCGGTGTTTTACGCGCAATGTCCAG‑CGACCCCGTGCGGCGCACCAGCCAGTCGAAAAAAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca                                                                                                    <  2:1079798/146‑1 (MQ=255)                                                              ccccaccGACAGAAACACCCCCATTACGTTCGGTGTTTTACGCGAAATGTTCAG‑CGTCCCCGGGCGGCGCACCGGCCATTCGAAAATAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                                                                                                <  1:1137540/144‑1 (MQ=255)                                                                  accGACAGAAACATCCCCATTACGAGCGGTGTTTTCCGCGCATTGTTCAG‑CGACCCGGGGCGGCGCACCAGCCAGTGGAAAATAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa                                                                                       <  2:101082/149‑1 (MQ=255)                                                                      aCAGAAACAGCCCCATAACGAGGGGTGTTTTACGCGAAATGTCCAG‑CGACCCCGGGCGGCGCATCAGCCAGTGGAAAATAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                                     <  2:957299/149‑1 (MQ=255)                                                                       cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  >  1:252857/1‑151 (MQ=255)                                                                       cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  >  1:711552/1‑151 (MQ=255)                                                                           accATGCCCTTTAGGACCGTTGTTTAACGGGCAATGTCCGC‑CGACCCGGGGCGCCGCATCAGCCAGCGGAAAAAAAC‑CCCCCCCAACACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACa                                                                              <  1:974069/149‑1 (MQ=255)                                                                           aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                               >  1:329602/1‑150 (MQ=255)                                                                              gtccccTTTGCGATCGGTTTTTTCGGCGCAATTTAAAG‑GCCCCCCGTGCGGCCAATCAGCCGGCAGAAAAAAAC‑CCCCCCAAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACAt                                                                             <  2:85820/146‑1 (MQ=255)                                                                                 cccATTACGATCGGTGTTTTACGCGAAATGTTCAG‑CGACCCCGGGCGGCGCATCAGCCAGTCGAAAATAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                         <  2:737984/150‑1 (MQ=255)                                                                                  ccATTCCGATCGGTTTTTTACGCGCATTTTCCAG‑CGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                       <  2:169194/151‑1 (MQ=255)                                                                                  ccATTACGATCGTTTTTTTACGCGAAATGTTCAG‑CGCCCCCGTGCGGCCCACCACCCAGGCGAAAAAAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                         <  1:137650/149‑1 (MQ=255)                                                                                     ttgcggtCGGTGTTTTCCGCACATTGTGCCG‑CGCCCCCGGGCGGCAAATCGGCCAGTCGAAAAAATC‑CCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTg                                                                     <  1:223336/144‑1 (MQ=255)                                                                                     ttccgttgggTTTTTTACGCCCATTGTCCAG‑CGACCCGGTGCGGCCCACCACCCAGTGGAAAAAACC‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTg                                                                     <  1:1030288/142‑1 (MQ=255)                                                                                          aTCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑c                                                               >  1:498109/1‑150 (MQ=255)                                                                                            gggTTTTTGACCCGAAATTTCCAG‑CGACCCGGTGCGGC‑CAACAGGCCGGCGGAAAAAATCCCCACCCAGCCCTCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTg                                                                 <  2:725664/146‑1 (MQ=255)                                                                                             gttgttttaccgaaaaTGTCCGG‑CGACCCCGGGCGGCGCACCACCCATGCGAAAATAAC‑CCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑Ca                                                              <  1:863495/134‑1 (MQ=255)                                                                                             ggTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CAcc                                                            >  2:936908/1‑150 (MQ=255)                                                                                               tgtTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CAcc                                                            >  1:669460/1‑148 (MQ=255)                                                                                                gtTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CAcc                                                            >  2:920455/1‑147 (MQ=255)                                                                                                 ttttACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGAt                                                        >  1:161817/1‑150 (MQ=255)                                                                                                  tttACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATa                                                       >  1:1227762/1‑150 (MQ=255)                                                                                                  tttACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATa                                                       >  2:115093/1‑150 (MQ=255)                                                                                                   ttACGCGAAATTTCCGG‑CGCCCCCGTGCGGCGCATCAGCCAGTGGAAAAAAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATaa                                                      <  1:1323838/150‑1 (MQ=255)                                                                                                     aCGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAAcc                                                   >  2:735734/1‑151 (MQ=255)                                                                                                      cgcgAAATGTCCGG‑CGCCCCCGGGCGGCCCACCAGCCGGCCGAAAATAAC‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAAc                                                    <  1:9100/149‑1 (MQ=255)                                                                                                         cgaaagTTCAGCCGACCCCGGGCGGCCCACCACCCAGTGGAAAAAAAC‑CCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAAcc                                                   <  1:316948/143‑1 (MQ=255)                                                                                                            aatttccagctcccccGGGCGGCGCATCAGCAAGTGGGAAATAAC‑CCCCCCAACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAAcc                                                   <  2:892855/132‑1 (MQ=255)                                                                                                             tGTCCAG‑CGACCCCGTGCGGCGCACCAGCAAGTCGAAAATAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAAcc                                            <  2:1049357/150‑1 (MQ=255)                                                                                                                acAG‑CGATCCGGGGGGGCGCATCACCCAGGCGAAAATAAC‑CCCCCCAACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCa                                           <  1:653800/147‑1 (MQ=255)                                                                                                                 cAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGt                                        >  2:726065/1‑150 (MQ=255)                                                                                                                 cAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGt                                        >  2:891903/1‑150 (MQ=255)                                                                                                                  aG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGta                                       >  2:89658/1‑150 (MQ=255)                                                                                                                   g‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGtaat                                     >  2:539498/1‑151 (MQ=255)                                                                                                                    gcccccccGGGCGGCGCATCACCCAGCGGAAAATAAC‑CCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGtaa                                      <  2:703394/145‑1 (MQ=255)                                                                                                                     cGACCCCGTGCGGCGCATCACCCGGCGGAAAAAAAC‑CCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCa                                           <  2:1085788/144‑1 (MQ=255)                                                                                                                     cGACCCCGGGCGGCCCAACACCCAGTCGAAAAAAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGtaat                                     <  1:1292567/150‑1 (MQ=255)                                                                                                                      gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGtaata                                    >  1:887254/1‑150 (MQ=255)                                                                                                                      gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACAGCGGGAACGGAGGCACCAAAGCCGCCTTTAAGAACCGACATCCCG‑GTGACTGCCACCAAATAAACGGGGACCCAGGtaat                                     >  2:432580/1‑149 (MQ=255)                                                                                                                              gcggcgCATCAGCCAGTCGAAAATAAC‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAaagaag                              <  2:678526/148‑1 (MQ=255)                                                                                                                                ggcgCACCAGCCAGTCGAAAATAAC‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGaagt                             <  2:872702/147‑1 (MQ=255)                                                                                                                                ggccCACCACCCAGTCGAAAATAAC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGaagtaa                           <  1:961677/149‑1 (MQ=255)                                                                                                                                  cgcATCACCCAGCGGGAAAAACC‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAaagaa                               <  1:276527/143‑1 (MQ=255)                                                                                                                                       cAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg                  >  1:357207/1‑151 (MQ=255)                                                                                                                                              tCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc           >  2:126507/1‑151 (MQ=255)                                                                                                                                                 aaaaTAAC‑CCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg         <  2:530212/149‑1 (MQ=255)                                                                                                                                                  aaaaaaC‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt             <  1:6392/145‑1 (MQ=255)                                                                                                                                                   aataac‑CCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTAc      <  2:134663/151‑1 (MQ=255)                                                                                                                                                       at‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc  >  2:301774/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC  >  NC_000913/2921023‑2921322

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