Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	45.5%	99.4 / 19.3	66	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/30);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 8.50e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACT  >  NC_000913/2921027‑2921319                                                                                                                                                         |                                                                                                                                               cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGTCAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑ccca                                                                                                                                                <  1:936085/151‑1 (MQ=255)  aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑cccac                                                                                                                                               >  2:335388/1‑151 (MQ=255)    cgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑cccaccc                                                                                                                                             <  2:526518/151‑1 (MQ=255)     gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑cccaccc                                                                                                                                             >  2:102895/1‑150 (MQ=255)      ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑cccaccc                                                                                                                                             >  2:202481/1‑149 (MQ=255)        aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGc                                                                                                                                          >  2:283239/1‑150 (MQ=255)             gCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACAccg                                                                                                                                    <  1:1107416/151‑1 (MQ=255)             gCCGATGATCAGCCACTCAACGTTGACGTAGTTGCAGAACCCCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACAccg                                                                                                                                    <  2:1021685/151‑1 (MQ=255)               cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACAccgcc                                                                                                                                  >  1:974196/1‑151 (MQ=255)                   gatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGAt                                                                                                                               >  2:2520/1‑150 (MQ=255)                      catcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATaaa                                                                                                                            >  2:1013023/1‑150 (MQ=255)                          cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                       >  1:515976/1‑151 (MQ=255)                            cTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAaaaccgca                                                                                                                      >  1:152560/1‑150 (MQ=255)                               aaCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAaaaccgcaaacc                                                                                                                  >  2:338614/1‑151 (MQ=255)                               aaCGTGGACGTAGTTGCAGAACCCCACCCCCATCAACAGCAACATCCCCATTACGAGCGGTTTTTTACGCGCAATTTCCAGCGATCCCGTGGGGCGCATCAGCCAGTCGAAAATAA‑C‑CCCCCCCAGCACACCGCCGATAaaaccgcaaa                                                                                                                    <  1:896039/149‑1 (MQ=255)                                 cgtgcgatatgTTGCAAAACACCAGCACCGCCGACACCAACCCGCCCATTCCGATCGTTGTTTGCCGCCCAATGTTCAGCCGCCCCGTGCGGCGCATCGGCCGTCCGAAAATAA‑C‑CCCCCCCAGCCCGCCGCCGATAaaaccgcaaaccg                                                                                                                 <  1:1075914/140‑1 (MQ=255)                                     gccGTAGTTGCAGAACCCCACCCCCATCGACACCAACATGCCCATTACGACCGGTGTTTTACGCGCAATGTCCAGCGACCCGGGGCGGCGCATCAGCCAGTCGGAAATAA‑C‑CCCCCCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                                 <  1:370391/144‑1 (MQ=255)                                       cGTAGCTGCAGAACACCATCACCATCGCCAGCAACATGTCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGTCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGCTAAAACCGCAAACCGCcggaac                                                                                                          >  1:908115/1‑151 (MQ=255)                                        gttgTTGCAAAACTCCACCACCCTCAACAGAAACAGGCCAATTACGATCGGTGTTTAACGCGCAGTGTTAAGCGTTCCCGGGCGGCGCATCGGCCAGTGGAAAAAAC‑C‑CCCCCCCAAACCCCCGCCGATAAAACCGCAAACCGCcggaac                                                                                                          <  2:1047701/147‑1 (MQ=255)                                        gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                         >  2:360663/1‑151 (MQ=255)                                         tggtggaaaaaCCCCATCCCCATCAACACAAAAACGCCCATACCGAGGGGTGTTTTCCCCCCAATTTCCAGGGACCCGGTGCGCCAAACCGGCCGGTCGAAAAAAA‑C‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCcggaac                                                                                                          <  1:154375/140‑1 (MQ=255)                                           gTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga                                                                                                       >  1:165387/1‑150 (MQ=255)                                           gTTGCAAACCCCCACCCCCGTCGACGTCAACCTGCCAATTACGAGCTGTGTTTTCCGGGAAATATCCAGCGAGCCGGTCGGGCGCGTCACCAAGTGGGAAATAA‑C‑CCCCCCCGTCACGCCGCCGATAAAACCGCAAACCGCcgga                                                                                                            <  2:12553/145‑1 (MQ=255)                                            ttccaaaaCCCCACCCCCATCGACAGCAACATCCCCATTCCGATCGGTTTTTCCCGCCAAATGTTCAGCGCTCCCGGGCGGCGCATCACCCGGTCGAAAAAAC‑C‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                      <  2:49419/144‑1 (MQ=255)                                            tgacaaaaCCCCCTCCCCATCGACAGCACCAGCCCCATTACGACCGGTTTTTTCCGCGCATGTTTCAGCTACCCCGTGGGGCCCACCACCCAGTGGAAAAAAA‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCcggaac                                                                                                          <  2:1045565/140‑1 (MQ=255)                                                      ccatccccaccGACACCACCAGCCCAATTCCGATCGTTGTTTCCCGCGAATTGTCCAGCGTCCCCGTGCGGCGCACCAGCCAGTCGAAAATAA‑T‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                                                                                            <  2:789692/150‑1 (MQ=255)                                                      ccacccccaccGACAGAAACAGCCCCTTTACGTTGGGTGTTTTCCGCAAAATTTCCAGGGACCCGGGGCGGCCCTTCACCCAGCCGAAAAAAA‑C‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg                                                                                                     <  1:922300/135‑1 (MQ=255)                                                             catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                    >  1:898868/1‑151 (MQ=255)                                                             catcGACACCAACATGCCCATTACGATCGGTTTTTTCCGCGCATTGTCCAGCGACCCGGTGCGGCGCATCAGCCAGTCGAAAATAA‑C‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                    <  2:1167550/151‑1 (MQ=255)                                                                    aGCAACATCCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGGGGGGCGCATCAGCCAGTGGAAAAAAA‑C‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                              <  2:145620/150‑1 (MQ=255)                                                                       aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                           >  2:1104189/1‑150 (MQ=255)                                                                       aaCATCCCCATAACGATCGTTTTTTCACGCGAAATGTTCAGCGACCCCGGGCGGCGCATCAGCAATTCGAAAAAAC‑C‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                           <  2:440010/150‑1 (MQ=255)                                                                            gCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                      >  1:329734/1‑150 (MQ=255)                                                                             cccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                    >  2:640741/1‑151 (MQ=255)                                                                              ccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                   >  1:1074022/1‑151 (MQ=255)                                                                                aTTACGATCGGCGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGt                                                                  >  1:96741/1‑150 (MQ=255)                                                                                  tACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGcc                                                               >  2:945273/1‑151 (MQ=255)                                                                                     gATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                            >  2:554335/1‑151 (MQ=255)                                                                                        cGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                         >  2:860164/1‑151 (MQ=255)                                                                                           tgtTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                         >  1:872510/1‑148 (MQ=255)                                                                                             ttttCCGCGAATGGTCCAGCGCCCCCGTGCGGCCCACCCCCCAG‑CGGAAAAAACC‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAg                                                       <  1:202479/148‑1 (MQ=255)                                                                                                 aCGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                >  1:632273/1‑151 (MQ=255)                                                                                                  cgcgCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                               >  1:842906/1‑151 (MQ=255)                                                                                                  cgcgCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                               >  2:1167040/1‑151 (MQ=255)                                                                                                   gcgcAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                              >  2:778143/1‑151 (MQ=255)                                                                                                   gcgaattttccAGCGTTCCGGTGGGCCGAACCACCCAGACGAAAAAAA‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                                 <  2:142010/138‑1 (MQ=255)                                                                                                    cacattGTCCAGCGCCCCCGGGCGGCGCACCAGCCAGTCGAAAATAA‑C‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                               <  1:132473/144‑1 (MQ=255)                                                                                                       attttccggctaTCCCGGGCGCCGGATCCCCCAGCCGGAAAAAA‑C‑CCCCCCTACCAACCCGCCGATAAACCCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                            <  2:238486/138‑1 (MQ=255)                                                                                                       aatttccAGCGACCCCGTGCGGCCCATCAGCCAG‑CGGAAAAAACC‑CCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                          <  2:1231471/145‑1 (MQ=255)                                                                                                       aatttccAGCGACCCCGGGCGGCCCTCCAGCCAGCGGAAAAAAA‑C‑CCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                               <  2:1059424/140‑1 (MQ=255)                                                                                                       aaTGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                          >  2:470055/1‑151 (MQ=255)                                                                                                         tGTCCACCGACCCCGGGCGGCGCACCAGCCAGTCGAAAATAT‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                               <  2:1292167/144‑1 (MQ=255)                                                                                                           ttCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                      >  2:1177111/1‑151 (MQ=255)                                                                                                            tCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                     >  1:121941/1‑151 (MQ=255)                                                                                                                cgcccccGGGCGGCTCCACAGCCAGT‑GGAAAAAA‑TCCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                     <  1:779890/143‑1 (MQ=255)                                                                                                                cGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCCTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                  >  1:651384/1‑150 (MQ=255)                                                                                                                 gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                                 >  2:923044/1‑150 (MQ=255)                                                                                                                  aTCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                                >  2:895555/1‑150 (MQ=255)                                                                                                                    cccgcggggccacaaaGCCAGTCGGAAAAAA‑C‑CCCACCCGGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                             <  2:335061/136‑1 (MQ=255)                                                                                                                    cccGGGCGGCGCACCAGCCGGTCGAAAATAA‑C‑CCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                      <  2:999675/142‑1 (MQ=255)                                                                                                                     cggggcggcgCATCACCCGATCGAAAAAAA‑C‑CCCACCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa                            <  1:365558/147‑1 (MQ=255)                                                                                                                        ggcggccCATCAGCCAGCCGAAAAAAC‑C‑CCCCCCCAACCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa                            <  2:192599/147‑1 (MQ=255)                                                                                                                         gcggcgCAACCCCCGGTCAAAAAAAA‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                          <  1:1288306/149‑1 (MQ=255)                                                                                                                                  cacccAGTGGAAAATAC‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                 <  1:605041/146‑1 (MQ=255)                                                                                                                                  cAGCCAGAAGAAAATAC‑C‑CCCCCCAAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAg                   <  1:1150477/147‑1 (MQ=255)                                                                                                                                   aGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGt              >  2:1236772/1‑151 (MQ=255)                                                                                                                                     ccAGTGGAAAATAA‑C‑CCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGt              <  2:74790/149‑1 (MQ=255)                                                                                                                                     ccAGTCGAAAAAAA‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAATCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt             <  1:43997/150‑1 (MQ=255)                                                                                                                                         tCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATa         >  2:1001156/1‑150 (MQ=255)                                                                                                                                             aaaaaa‑C‑CCCCCCCACCACCCCGCCGATAAACCCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTa    <  1:794926/151‑1 (MQ=255)                                                                                                                                             aaaTAC‑C‑CCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTa    <  1:633668/151‑1 (MQ=255)                                                                                                                                                aaaacccccccccaCCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt     <  1:731442/143‑1 (MQ=255)                                                                                                                                                 aaaaccccccccaGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt  <  2:152562/146‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               CAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACT  >  NC_000913/2921027‑2921319

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