Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,173	0	T	C	41.1%	101.9 / 32.7	73	I310V (ATT→GTT)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (35/8);  new base C (0/30);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTA  >  NC_000913/2921024‑2921317                                                                                                                                                       |                                                                                                                                                  ggCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑c                                                                                                                                                 <  2:656209/151‑1 (MQ=255) ggCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑c                                                                                                                                                 <  2:8321/151‑1 (MQ=255)    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)                                                                                                                                                       |                                                                                                                                                  GGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑T‑CCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTA  >  NC_000913/2921024‑2921317

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