Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	53.7%	80.0 / 108.4	67	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/36);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 7.66e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-12

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGAT  >  NC_000913/2921028‑2921326                                                                                                                                                        |                                                                                                                                                       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                                                           >  2:115093/1‑150 (MQ=255)                                                                                             tttACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATa                                                           >  1:1227762/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:891903/1‑150 (MQ=255)                                                                                                            cAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGt                                            >  2:726065/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)                                                                                                                                                     cccccccaCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc      <  2:993779/149‑1 (MQ=255)                                                                                                                                                     cccccccaCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGAAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGa       <  2:442659/148‑1 (MQ=255)                                                                                                                                                     cccacccaGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCg    >  2:852619/1‑151 (MQ=255)                                                                                                                                                       cacccaGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGat  >  1:92564/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGAT  >  NC_000913/2921028‑2921326

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