BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,907,998	A→T	15.5%	intergenic (‑57/+31)	eno ← / ← pyrG	enolase/CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,907,998	0	A	T	15.5%	398.1 / 57.4	160	intergenic (‑57/+31)	eno/pyrG	enolase/CTP synthetase
Reads supporting (aligned to +/- strand): ref base A (64/71); new base T (21/4); total (85/75)
Fisher's exact test for biased strand distribution p-value = 8.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 160 total aligned reads displayed.

GTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCAT  >  NC_000913/2907900‑2908096

 tCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTGACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTgg                                                                                                  >  2:2757371/1‑100 (MQ=255)

  cGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAggg                                                                                                 >  2:3115650/1‑100 (MQ=255)

  cGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTACAGACAAACAACGCGTACCCAggg                                                                                                 >  1:856678/1‑100 (MQ=255)

  cGAGGATTTCACGACCGAAGATTAATACGATTTTGGACATTAGGTTTTCGTCAAGTCACTAGTTAAACTGAAAATCCAGACAAACAACGCGTACCCAggg                                                                                                 >  2:3555752/1‑100 (MQ=255)

   gatgatTTCACGACCGATGATTTTTACGATTTTGTACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAggg                                                                                                 <  1:818957/99‑1 (MQ=255)

    atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                               <  1:1969550/100‑1 (MQ=255)

      gatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTAcg                                                                                             >  2:74480/1‑100 (MQ=255)

       atTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcg                                                                                             <  1:2468608/99‑1 (MQ=255)

       atTTCACCACCGATGATTTTTACGATTTTGTACAGTAGGTTTTCCTCAAGACACTAGTTAAGCTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcg                                                                                             <  1:1841246/99‑1 (MQ=255)

        tttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                           <  2:1530640/100‑1 (MQ=255)

        tttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                           <  2:1876851/100‑1 (MQ=255)

          tCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGtt                                                                                         <  1:3471701/100‑1 (MQ=255)

          tCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGtt                                                                                         <  1:874525/100‑1 (MQ=255)

           cACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTg                                                                                        <  2:561247/100‑1 (MQ=255)

           cACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTg                                                                                        >  2:2253373/1‑100 (MQ=255)

            aCGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgt                                                                                       <  2:2097354/100‑2 (MQ=255)

             cGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGcc                                                                                      >  2:362117/1‑100 (MQ=255)

               aCCGATGATTTTTAAGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGAATACCCAGGGTACGCGTTGcc                                                                                      >  2:745682/1‑98 (MQ=255)

                 cGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGctc                                                                                  <  2:210789/100‑1 (MQ=255)

                  gatgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGctct                                                                                 >  1:790987/1‑100 (MQ=255)

                   atgatTTTTACGATTTTGTACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGctct                                                                                 <  1:1394409/99‑1 (MQ=255)

                   atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                <  2:3287770/100‑1 (MQ=255)

                      atTTTTACGATTTTTGACATTAGGTTTTCCTCAAGTCATTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAAc                                                                              <  1:3324625/99‑1 (MQ=255)

                       tttttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGAGTACCCAGGGTACGCGTTGCCGCTCTAACtt                                                                            >  2:939064/1‑100 (MQ=255)

                           tACGATTTTGTACATTACGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGTTACGCGTTGCCGCTCTAACtttttt                                                                        <  2:2811918/100‑1 (MQ=255)

                            aCGATTTTGGACATTATTTTTTCTTTTTGTTATTATTTATACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                       <  1:3040898/100‑1 (MQ=255)

                             cGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTTATCTGAATCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                      <  1:1138386/100‑1 (MQ=255)

                             cGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                      <  1:2606656/100‑1 (MQ=255)

                             cGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAAAAATGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                      >  2:2747454/1‑100 (MQ=255)

                              gATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTTAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttact                                                                     <  1:2747478/100‑1 (MQ=255)

                               agtttGGACATTAGGTTTTCCTCAAGTCACTAGTTGCAGTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTTCCGCTCTAACTTTTttactt                                                                    <  2:2291100/98‑1 (MQ=255)

                               aTTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactt                                                                    <  1:2196836/100‑1 (MQ=255)

                                ttttGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                   <  1:2850549/100‑1 (MQ=255)

                                     gACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTacttactt                                                                >  1:2119408/1‑98 (MQ=255)

                                       cATTAGGTTTTCCTCAAGTCACTAGTTAAATTGAAATTTCATTTTAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                            <  1:259907/100‑1 (MQ=255)

                                       cATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                            <  1:362123/100‑1 (MQ=255)

                                       cATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                            <  1:2993122/100‑1 (MQ=255)

                                        aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCt                                                           <  1:953719/100‑1 (MQ=255)

                                        aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCt                                                           <  2:856668/100‑1 (MQ=255)

                                         ttAGGTTTTCCTCAAGTCACTAGTTATATTGAAATTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTg                                                          <  1:1474222/100‑1 (MQ=255)

                                         ttAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTg                                                          <  1:1826139/100‑1 (MQ=255)

                                          tAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGa                                                         <  2:437422/100‑1 (MQ=255)

                                          tAGGTTTTCCTCAAGTCAATAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGa                                                         <  1:1791054/100‑1 (MQ=255)

                                           aGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                        >  2:2069556/1‑100 (MQ=255)

                                            ggTTTTCATCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGATACACATTACCGCTCTAACTTTTTTACTTACTTCGCCTGACg                                                       >  1:2041956/1‑100 (MQ=255)

                                                  ccTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCACTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTg                                                 >  1:449619/1‑100 (MQ=255)

                                                   cTCAAGGCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCACTAACTTAAATACTTAATTAACATAACGTATCTg                                                 >  1:263149/1‑99 (MQ=255)

                                                       aGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACt                                            <  1:2804477/100‑1 (MQ=255)

                                                       aGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACt                                            >  2:2980112/1‑100 (MQ=255)

                                                          cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGc                                         >  2:1412001/1‑100 (MQ=255)

                                                          cACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTGTACTTACTTCGCCTGACGTTTCTGGAACTCga                                         <  2:2762814/100‑2 (MQ=255)

                                                           aCTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCt                                        >  2:2958467/1‑100 (MQ=255)

                                                           aCTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCt                                        >  1:2910774/1‑100 (MQ=255)

                                                            cTAGTTAAACTGAAACACCAGACAAACAACGCGTACCCATGGTATGCGTTGCCGCTCTAACTTTTTTACTTACTTAGCCTGACGTTTCTGTTATTCGCt                                        >  2:65286/1‑99 (MQ=255)

                                                             tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTgg                                      >  2:3760639/1‑100 (MQ=255)

                                                               gTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGAAACTCGCTggcg                                    >  1:3743371/1‑100 (MQ=255)

                                                                 tAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCAGCTCTAACGTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcggc                                  >  2:3855851/1‑100 (MQ=255)

                                                                     cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                              <  1:2302431/100‑1 (MQ=255)

                                                                     cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                              <  1:357819/100‑1 (MQ=255)

                                                                     cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                              <  2:2567576/100‑1 (MQ=255)

                                                                     cTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                              >  2:1840796/1‑100 (MQ=255)

                                                                      tGAAACTCCAGACAAGCAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                             <  2:2832663/100‑1 (MQ=255)

                                                                      ctaaaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                             <  2:3156706/98‑1 (MQ=255)

                                                                       gAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcac                            <  2:1701333/100‑1 (MQ=255)

                                                                        aaaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcaca                           <  2:1175060/100‑1 (MQ=255)

                                                                        aaaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcaag                           >  1:3156736/1‑98 (MQ=255)

                                                                         aaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaa                          <  1:2790002/100‑1 (MQ=255)

                                                                              cAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCt                     <  2:3055019/100‑1 (MQ=255)

                                                                              cAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCt                     >  2:2942896/1‑100 (MQ=255)

                                                                               aGACAAATAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTg                    <  2:2001259/100‑1 (MQ=255)

                                                                               aGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTg                    >  1:1576178/1‑100 (MQ=255)

                                                                               aGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCGTTCATAAAGCCTg                    >  1:3526800/1‑100 (MQ=255)

                                                                                  gcggcaacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaa                 >  2:231476/5‑100 (MQ=255)

                                                                                  gcggcaacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaa                 >  1:2702113/5‑100 (MQ=255)

                                                                                  caaacaaCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaa                 >  1:1203002/1‑100 (MQ=255)

                                                                                   cggcaccgcgTCCCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACGTACTTCCCCTGACGTTTCTGGAACTCGCTAGCGGCTTTCACAAACCCTGCaaa                >  2:3325248/7‑100 (MQ=255)

                                                                                   aaacaaCGCGTACCCAGGGTACGTGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                <  1:840866/100‑1 (MQ=255)

                                                                                    aacaacGCGTACCCAGTGTACGCGTTGCCGCTCTACCTTTTTTACTTACTGCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAAc               <  2:1784026/100‑1 (MQ=255)

                                                                                    aacaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACAATACAGGAACTCGCTGGCGACTTTCACAAAGCCTGCAAAc               >  1:3327137/1‑100 (MQ=255)

                                                                                    aacaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAAGTTTCTGGAACTCACTGGCGGCTTTCACAAAGCCTGCaaa                >  1:327323/1‑99 (MQ=255)

                                                                                     gcaacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACa              >  1:1783337/2‑100 (MQ=255)

                                                                                     gcaacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACa              >  2:3660138/2‑100 (MQ=255)

                                                                                     gcaacGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACa              >  2:1168051/2‑100 (MQ=255)

                                                                                        acGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCg           >  2:2697910/1‑100 (MQ=255)

                                                                                        acGCGTACCCTGAGTACGCGTTGCCGACCTAACTTTTTTACTTACTTCGCCTGACGTTTGTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCg           >  2:1043194/1‑100 (MQ=255)

                                                                                        acGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAg             >  1:3798942/1‑98 (MQ=255)

                                                                                         cgcgTACCCTGGGTACGCGTTGCCGCGCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGAGGCTTTCACAAAGCCTGCAAACAGCgg          >  2:1904081/1‑100 (MQ=255)

                                                                                          gcgTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCg           >  1:2862012/1‑98 (MQ=255)

                                                                                           cgTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCATTCAAAAAGCCTGCAAACAGCagg         >  1:3626303/1‑99 (MQ=255)

                                                                                            gTACCCAGTGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCggg         <  2:2556838/98‑1 (MQ=255)

                                                                                            gTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAGCTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTg       >  2:2748342/1‑100 (MQ=255)

                                                                                             tACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACGAAGCCTGCAAACAGCGGGTGa      >  2:685639/1‑100 (MQ=255)

                                                                                             tACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGa      >  2:3798315/1‑100 (MQ=255)

                                                                                              aCCCAGGGTTCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc     >  2:152784/1‑100 (MQ=255)

                                                                                              aCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc     >  2:3267539/1‑100 (MQ=255)

                                                                                              aCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc     <  2:1202989/100‑1 (MQ=255)

                                                                                               cccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGTCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc     <  2:784594/99‑1 (MQ=255)

                                                                                                ccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa   <  2:1665278/100‑1 (MQ=255)

                                                                                                ccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa   <  2:1425947/100‑1 (MQ=255)

                                                                                                 cAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCAt  >  1:3138014/1‑100 (MQ=255)

GTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCAT  >  NC_000913/2907900‑2908096

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