BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,703,593	A→G	18.6%	I239T (ATC→ACC)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,593	0	A	G	18.6%	193.0 / 38.4	98	I239T (ATC→ACC)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (51/28); new base G (13/5); total (64/34)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 101 total aligned reads displayed.

GCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTG  >  NC_000913/1703496‑1703690

gCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGata                                                                                                 >  1:319536/1‑100 (MQ=255)

 ccccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGatat                                                                                                >  2:493097/1‑100 (MQ=255)

 ccccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGatat                                                                                                >  2:1940781/1‑100 (MQ=255)

  cccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACAAtata                                                                                               >  2:1635048/1‑100 (MQ=255)

   ccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGtatat                                                                                              <  1:759227/100‑2 (MQ=255)

     aGTCTTCCTGCGGTAGACGCGCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCg                                                                                            <  1:1955677/100‑1 (MQ=255)

        cTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGt                                                                                           >  1:4225819/1‑98 (MQ=255)

         ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGc                                                                                         <  1:3235560/99‑1 (MQ=255)

         ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCt                                                                                        <  1:3587005/100‑1 (MQ=255)

           ccTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                      <  1:3379149/100‑1 (MQ=255)

            cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatgg                                                                                     <  1:2365351/100‑5 (MQ=255)

            cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                     <  1:1400297/100‑1 (MQ=255)

            cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                     <  1:1420863/100‑1 (MQ=255)

            cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                     >  2:3847394/1‑100 (MQ=255)

             tGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatg                                                                                      <  1:1744711/98‑4 (MQ=255)

               cGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCa                                                                                  >  2:2780049/1‑100 (MQ=255)

                ggTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGACAGAGGCCATATTTGAAATAACTGCCTGGGGATGCATGCACGAAGTACCGTGCTGACt                                                                                    >  2:3947668/1‑97 (MQ=255)

                 gTAGACGCTCGCCATTTGTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCTCATAACTGCCTCGGTATCCATGCACGGTATATCGTGcatgggtccc                                                                                <  1:373254/100‑10 (MQ=255)

                  tAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCg                                                                               >  1:4048564/1‑100 (MQ=255)

                   aGACGCTCGCCATTTTTCGGACGTTCTGCCGGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                              <  1:4304371/100‑1 (MQ=255)

                   aGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                              <  1:2046210/100‑1 (MQ=255)

                   aGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                              <  2:319537/100‑1 (MQ=255)

                      cccTCGCCATTTTTCAGACGTTCTTCCTGGGGCTCGAGGCCATATTTCACATACCTCCCTCGGGATCCATGCACGTTATACCGTGCTGACTCAGCGGCTg                                                                           <  1:112577/98‑1 (MQ=255)

                        cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                         >  1:4019660/1‑100 (MQ=255)

                        cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                         <  1:3519677/100‑1 (MQ=255)

                        cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                         <  1:2513150/100‑1 (MQ=255)

                         tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCa                                                                        <  1:4399128/100‑1 (MQ=255)

                            ccATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCa                                                                     >  1:1553421/1‑100 (MQ=255)

                              aTTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATg                                                                   >  2:1750170/1‑100 (MQ=255)

                                  ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                               >  1:541544/1‑100 (MQ=255)

                                  ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                               <  1:3506927/100‑1 (MQ=255)

                                  ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                               <  1:3386453/100‑1 (MQ=255)

                                         gTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCCGCATGGTACCGTGTaa                                                        >  2:3499881/1‑100 (MQ=255)

                                          ttctccCTGTGGATCAAGGCCATATTCCACATACCTGCCACGAGCCCCATGCACGATATACCGTGCTCACTCCTCGGCTGCCAGCATGGTACCGTGTaa                                                        <  1:3013243/99‑1 (MQ=255)

                                            cttcCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAAt                                                     >  2:3946627/1‑100 (MQ=255)

                                             ttcCTGTGGATCGAGGCCATATGTCACATGACTGCCGCGGGATCCGTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCATGTaaaa                                                      <  1:3766004/98‑1 (MQ=255)

                                                cTGTGGATCGAGGCCATATTTCACATAACTGCCGCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTATAAAATGAc                                                  >  1:2889351/1‑99 (MQ=255)

                                                 tgtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                >  2:3763934/1‑100 (MQ=255)

                                                 tgtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                >  2:1717918/1‑100 (MQ=255)

                                                  gtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTCCCGTGTAAAATGACTcg                                               >  1:4053412/1‑100 (MQ=255)

                                                   tgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgc                                              >  1:1923059/1‑100 (MQ=255)

                                                     gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATCCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                            >  1:1456411/1‑100 (MQ=255)

                                                     gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                            >  1:4281303/1‑100 (MQ=255)

                                                     gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                            >  1:945105/1‑100 (MQ=255)

                                                      ntCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGCAAAATGACTcgccac                                           >  1:2669253/2‑98 (MQ=255)

                                                      aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccgc                                           >  2:3349271/1‑100 (MQ=255)

                                                         gAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtgt                                        >  1:673473/1‑100 (MQ=255)

                                                         gAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtgt                                        <  2:1923066/100‑1 (MQ=255)

                                                           ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTgg                                      >  1:112990/1‑100 (MQ=255)

                                                           ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTgg                                      >  2:35516/1‑100 (MQ=255)

                                                            gCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTCACTCAGCGGCTGCCAGCATGGTACCCTGCAAAATGACTCGCCGCt                                          >  1:2119017/1‑95 (MQ=255)

                                                              cATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGGATGGTACCGGGTAAAATGATTCGCCGCTGTgg                                      >  1:2598724/1‑97 (MQ=255)

                                                                tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAgg                                 <  1:2142065/100‑1 (MQ=255)

                                                                 atTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGa                                <  1:833123/100‑1 (MQ=255)

                                                                  tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAATATGACTCGCCGCTGTGGATAGGAc                               >  1:3263091/1‑100 (MQ=255)

                                                                  tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGAc                               >  2:3626561/1‑100 (MQ=255)

                                                                  tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGAc                               >  2:3246962/1‑100 (MQ=255)

                                                                     cacaTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAg                            <  1:1890581/100‑1 (MQ=255)

                                                                        aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                           <  1:2874551/98‑1 (MQ=255)

                                                                         tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATgg                        >  1:3596555/1‑100 (MQ=255)

                                                                         tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATgg                        >  2:1700036/1‑100 (MQ=255)

                                                                         tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                         >  1:4068178/1‑99 (MQ=255)

                                                                         tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                         >  1:3718174/1‑99 (MQ=255)

                                                                         tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGGGTAAAATGACTCGCCGCCGTGGATAGGACa                              >  1:3713639/1‑94 (MQ=255)

                                                                             tGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTg                    >  1:3836509/1‑100 (MQ=255)

                                                                              gCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTgg                   >  1:1808172/1‑100 (MQ=255)

                                                                              gCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTgg                   >  2:3979506/1‑100 (MQ=255)

                                                                               ccTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGa                  <  2:4019690/100‑1 (MQ=255)

                                                                                cTCGGGATCCATGCACGATATACCGTGCTGACTCAGTGGCTGCCAGCATGGTACCCTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTg                    <  1:2413126/97‑1 (MQ=255)

                                                                                 tCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGaaa                >  1:1590765/1‑100 (MQ=255)

                                                                                 tCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGa                  <  1:2738944/98‑1 (MQ=255)

                                                                                   gggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata              >  2:918302/1‑100 (MQ=255)

                                                                                   gctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata              >  2:3604573/10‑100 (MQ=255)

                                                                                   gctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACGCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata              >  1:2462578/10‑100 (MQ=255)

                                                                                    ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGAATGGTACCGTGTAAAATAACTCGCCAGTGTGGATCGGACAAGGTGGCGTAGAaata              >  2:2569618/1‑99 (MQ=255)

                                                                                    ctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTCAAATGACTCGCCGCTGTGGATAGGACAAAATGGCGTGGAaataa             >  1:4320748/9‑100 (MQ=255)

                                                                                    ctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaataa             >  2:336121/9‑100 (MQ=255)

                                                                                     gATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaataat            >  2:1409568/1‑100 (MQ=255)

                                                                                      gagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATc           >  2:3460941/7‑100 (MQ=255)

                                                                                      aTCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATc           <  2:1553427/100‑1 (MQ=255)

                                                                                       tCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg          >  1:1023631/1‑100 (MQ=255)

                                                                                       tCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg          <  1:3786790/100‑1 (MQ=255)

                                                                                       agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg          >  2:1572273/6‑100 (MQ=255)

                                                                                       agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg          >  2:495916/6‑100 (MQ=255)

                                                                                       agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg          >  2:494647/6‑100 (MQ=255)

                                                                                       agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATAGCGTGGAaat               >  1:2738054/6‑95 (MQ=255)

                                                                                        gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATTGTACCGTGTAAAATGACTCGCCGCTGGGGATAGGACAAGATGGCGTGGAAATAATCg          >  1:2444502/5‑99 (MQ=255)

                                                                                        gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg         >  1:960091/5‑100 (MQ=255)

                                                                                        gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg         >  2:258772/5‑100 (MQ=255)

                                                                                        gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGGGGAAATAATCgg         >  2:2776740/5‑100 (MQ=255)

                                                                                        ccATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAACGTGTAAAAGGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg         <  1:1050865/100‑1 (MQ=255)

                                                                                          aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTc       <  1:1087290/100‑1 (MQ=255)

                                                                                           ngCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCg      >  2:2198970/2‑100 (MQ=255)

                                                                                            gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGa     >  2:765128/1‑100 (MQ=255)

                                                                                            gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGa     >  1:1291993/1‑100 (MQ=255)

                                                                                             cACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGAc    >  2:3377849/1‑100 (MQ=255)

                                                                                             cACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGAc    >  2:2104904/1‑100 (MQ=255)

                                                                                              aCGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACt   >  2:2223628/1‑100 (MQ=255)

                                                                                              aCGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACt   >  2:2223712/1‑100 (MQ=255)

                                                                                               cGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTg  <  1:2484916/100‑1 (MQ=255)

GCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTG  >  NC_000913/1703496‑1703690

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