BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,636,979	A→T	20.3%	D144E (GAT→GAA)	yhiN ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,636,979	0	A	T	20.3%	378.3 / 96.8	183	D144E (GAT→GAA)	yhiN	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (64/81); new base T (10/27); total (74/109)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 185 total aligned reads displayed.

GCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTC  >  NC_000913/3636881‑3637077

gCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGAt                                                                                                   <  2:4062519/100‑1 (MQ=255)

gCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGAt                                                                                                   <  1:1655954/100‑1 (MQ=255)

  aaTCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCa                                                                                                 <  1:2978216/100‑1 (MQ=255)

      ttATAACCAAACGGCGACGCGCCCAGCCCCTCCATTGACGGCCCACCCGTCGCGATTACCAGCTTTTCGTAACCGACAGTCATGCCGTTCAGATCAAGCg                                                                                             <  2:1408955/100‑1 (MQ=255)

      ttATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCg                                                                                             >  1:261484/1‑100 (MQ=255)

         taACCAACCGGCGACGCGCCCAGCCCCGGCATTGGCAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTg                                                                                           <  2:4111202/99‑1 (MQ=255)

         taACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGa                                                                                          <  2:3035370/100‑1 (MQ=255)

           aCCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAg                                                                                        <  1:286061/100‑1 (MQ=255)

           aCCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAg                                                                                        <  1:3834293/100‑1 (MQ=255)

           aCCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAg                                                                                        >  1:2221274/1‑100 (MQ=255)

            ccAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGc                                                                                       <  1:1175670/100‑1 (MQ=255)

            ccAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGc                                                                                       <  1:4007398/100‑1 (MQ=255)

            ccAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGc                                                                                       <  1:4007637/100‑1 (MQ=255)

             cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGTCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGcc                                                                                      >  1:1996972/1‑100 (MQ=255)

             cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCTGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGcc                                                                                      <  2:3980328/100‑1 (MQ=255)

             cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGcc                                                                                      <  1:2655920/100‑1 (MQ=255)

             cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGcc                                                                                      >  1:1973368/1‑100 (MQ=255)

                   gCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTc                                                                                >  2:2056562/1‑100 (MQ=255)

                     gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCAt                                                                              <  1:461278/100‑1 (MQ=255)

                     gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCAt                                                                              <  2:1570837/100‑1 (MQ=255)

                     gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCAt                                                                              >  1:3377868/1‑100 (MQ=255)

                     gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCAt                                                                              <  1:1478207/100‑1 (MQ=255)

                       cgcgCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATcc                                                                            >  2:1015452/1‑100 (MQ=255)

                        gcgcCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCt                                                                           >  1:365599/1‑100 (MQ=255)

                         cgcCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCtt                                                                          >  2:975460/1‑100 (MQ=255)

                          gcCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTc                                                                         <  1:2167025/100‑1 (MQ=255)

                          gcCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTc                                                                         >  1:1426775/1‑100 (MQ=255)

                           cccAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCg                                                                        <  1:3121711/100‑1 (MQ=255)

                             cAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGcc                                                                      <  2:1512081/100‑1 (MQ=255)

                             cAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGcc                                                                      <  2:1942103/100‑1 (MQ=255)

                              aGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCca                                                                     <  1:4034686/100‑1 (MQ=255)

                               gCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCca                                                                     >  2:2341732/1‑99 (MQ=255)

                                ccccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCcaca                                                                   >  2:4173326/1‑100 (MQ=255)

                                 cccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCacac                                                                  <  2:2110095/100‑1 (MQ=255)

                                   cGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTc                                                                <  2:1780124/100‑1 (MQ=255)

                                   cGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTc                                                                <  1:273136/100‑1 (MQ=255)

                                     gCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAg                                                              <  1:2377995/100‑1 (MQ=255)

                                      cATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGc                                                             <  1:945972/100‑1 (MQ=255)

                                      cATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGc                                                             <  1:2695468/100‑1 (MQ=255)

                                        ttGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCAc                                                           >  2:2758626/1‑100 (MQ=255)

                                          gACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACtt                                                         >  2:3027588/1‑100 (MQ=255)

                                           aCAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTc                                                        <  2:2725576/100‑1 (MQ=255)

                                            cAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCg                                                       >  1:632476/1‑100 (MQ=255)

                                             aGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGc                                                      <  1:1287386/100‑1 (MQ=255)

                                               cccACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTa                                                    <  1:847948/100‑1 (MQ=255)

                                                ccaccaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTAc                                                   <  2:60414/100‑1 (MQ=255)

                                                ccaccaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTAc                                                   <  1:3264797/100‑1 (MQ=255)

                                                caacgaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTAc                                                   <  2:1715650/95‑1 (MQ=255)

                                                 caccaGTCGCGATAACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACg                                                  <  1:3455739/100‑1 (MQ=255)

                                                  accaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGc                                                 <  2:1526764/100‑1 (MQ=255)

                                                   ccaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCa                                                >  2:1110632/1‑100 (MQ=255)

                                                      gTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATc                                             <  2:4100766/100‑1 (MQ=255)

                                                      gTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATc                                             <  1:26318/100‑1 (MQ=255)

                                                      gTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATc                                             >  1:576079/1‑100 (MQ=255)

                                                      gTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCGGTGTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATc                                             >  2:195719/1‑100 (MQ=255)

                                                       tCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCt                                            <  1:1336245/100‑1 (MQ=255)

                                                        cgcgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTg                                           >  1:1478696/1‑100 (MQ=255)

                                                            aTGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAgg                                       <  1:2388339/100‑1 (MQ=255)

                                                            aTGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAgg                                       >  2:2648350/1‑100 (MQ=255)

                                                            aTGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAgg                                       <  2:674660/100‑1 (MQ=255)

                                                             tGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGt                                      <  1:20832/100‑1 (MQ=255)

                                                              gACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTc                                     <  2:2239715/100‑1 (MQ=255)

                                                                  aGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACAt                                 >  2:3236556/1‑100 (MQ=255)

                                                                   gCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACAtt                                >  1:2154232/1‑100 (MQ=255)

                                                                   gCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACAtt                                >  2:1866976/1‑100 (MQ=255)

                                                                      tttCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGcc                             >  2:2459011/1‑100 (MQ=255)

                                                                          gCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtc                         >  1:3569979/1‑100 (MQ=255)

                                                                          gCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTAGCTACGCAATATGAAGGTCAAATTGCCCCtc                         >  1:766885/1‑100 (MQ=255)

                                                                            aaCCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtctc                       <  1:3733507/100‑1 (MQ=255)

                                                                            aaCCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtctc                       <  1:2931730/100‑1 (MQ=255)

                                                                            aaCCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtctc                       <  2:2952105/100‑1 (MQ=255)

                                                                            aaCCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtctc                       <  1:2307887/100‑1 (MQ=255)

                                                                             aCCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCg                      <  1:3311616/100‑1 (MQ=255)

                                                                              ccGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGc                     <  1:2343821/100‑1 (MQ=255)

                                                                               cGACAGTCATGCGGGTCAGCTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGc                     <  2:1185455/99‑1 (MQ=255)

                                                                                gACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCAc                   <  1:3245885/100‑1 (MQ=255)

                                                                                gACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCAc                   <  1:3789931/100‑1 (MQ=255)

                                                                                gACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCAc                   <  2:1277389/100‑1 (MQ=255)

                                                                                 aCAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACt                  <  1:385146/100‑1 (MQ=255)

                                                                                    gTCATTCCGTTCAGATCAAACGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCAt               >  2:1252180/1‑100 (MQ=255)

                                                                                    gTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCAt               >  1:1606495/1‑100 (MQ=255)

                                                                                    gTCATGCCATTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCAt               >  1:3521386/1‑100 (MQ=255)

                                                                                     tCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATc              <  1:2501865/100‑1 (MQ=255)

                                                                                        tGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccac           >  1:799742/1‑100 (MQ=255)

                                                                                         gCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacc          <  1:1829313/100‑1 (MQ=255)

                                                                                         gCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacc          <  2:3590015/100‑1 (MQ=255)

                                                                                          ccGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacca         <  1:126089/100‑1 (MQ=255)

                                                                                           cGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAg        <  1:2454702/100‑1 (MQ=255)

                                                                                           cGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAg        >  1:3652160/1‑100 (MQ=255)

                                                                                            gTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGc       <  2:1387112/100‑1 (MQ=255)

                                                                                             ttCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCa      <  2:2019662/100‑1 (MQ=255)

                                                                                             ttCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCa      <  2:2159074/100‑1 (MQ=255)

                                                                                              tCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCAt     <  1:4123902/100‑1 (MQ=255)

                                                                                               cAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATg    <  2:703170/100‑1 (MQ=255)

                                                                                                 taTCTAGCTTCAAGCCTTGTACATCCTTCGCCACACACATCGCTTCGCTCCGCAATCTGAAACTCACATTCCCCTACGCACACCCTACCACCAGCa      >  2:2266756/2‑96 (MQ=255)

                                                                                                 naTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTCCGCAATCTGACCGTCCCATTCCCCTTCTCGCACTCATCCAGCAGCATGTc  >  2:2422623/2‑100 (MQ=255)

                                                                                                 gATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTc  >  2:3819973/1‑100 (MQ=255)

                                                                                                 gATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTc  >  1:1461658/1‑100 (MQ=255)

                                                                                                 gATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTc  <  1:784439/100‑1 (MQ=255)

                                                                                                 gATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTc  <  1:2199177/100‑1 (MQ=255)

GCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTC  >  NC_000913/3636881‑3637077

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