Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,166	0	G	A	44.1%	73.3 / 33.7	68	S312F (TCC→TTC)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (32/6);  new base A (0/30);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.19e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAC  >  NC_000913/2921016‑2921313                                                                                                                                                        |                                                                                                                                                      ccGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCgg                                                                                                                                                       <  1:610521/151‑1 (MQ=255)        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)                                                                                                                                                        |                                                                                                                                                      CCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAG‑CGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAAT‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTG‑CACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAC  >  NC_000913/2921016‑2921313

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: ‑