BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,481	(G)_6→7	coding (11/765 nt)	nfeF ←	NADPH‑dependent ferric chelate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,475	1	.	G	100.0%	1193.7 / NA	279	R6P (CGC→CCC)	nfeF	NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (105/174); total (105/174)

Only 100 of 324 total aligned reads displayed.

CCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATG  >  NC_000913/3216376‑3216570

ccgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑g                                                                                                <  1:5081353/101‑1 (MQ=255)

 cgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gg                                                                                               >  2:5081825/1‑101 (MQ=255)

 cgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gg                                                                                               <  1:5081242/101‑1 (MQ=255)

 cgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gg                                                                                               <  1:5081290/101‑1 (MQ=255)

  gccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggg                                                                                              >  1:5081827/1‑101 (MQ=255)

   ccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggg                                                                                             <  2:5081412/101‑1 (MQ=255)

    cgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggg                                                                                            >  1:5081832/1‑101 (MQ=255)

    cgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggg                                                                                            >  2:5081834/1‑101 (MQ=255)

     gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggggg                                                                                           >  2:5081839/1‑101 (MQ=255)

      aGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggggg                                                                                          >  1:5081842/1‑100 (MQ=255)

          cAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggggggtgtt                                                                                      <  2:5081440/101‑6 (MQ=255)

            aTGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTAt                                                                                     <  2:5081252/101‑1 (MQ=255)

              gcgcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTc                                                                                   >  1:5081847/1‑101 (MQ=255)

                gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAt                                                                                 >  1:5081851/1‑101 (MQ=255)

                 cTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAtt                                                                                <  1:5081257/101‑1 (MQ=255)

                 cTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAtt                                                                                >  2:5081852/1‑101 (MQ=255)

                 cTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAtt                                                                                <  1:5081221/101‑1 (MQ=255)

                 cTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAtt                                                                                <  2:5081640/101‑1 (MQ=255)

                   gAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGGGTTATTCATTGt                                                                              >  2:5081855/1‑101 (MQ=255)

                     aaaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTa                                                                            >  2:5081860/1‑101 (MQ=255)

                       aCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATc                                                                          >  2:5081865/1‑101 (MQ=255)

                       aCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggggggcg                                                                                        >  1:5081867/1‑83 (MQ=255)

                         cGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGc                                                                        <  2:5081435/101‑1 (MQ=255)

                         cGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGc                                                                        >  1:5081871/1‑101 (MQ=255)

                          ggCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggggg                                                                                          >  2:5081872/1‑80 (MQ=255)

                           gcgcTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCt                                                                      <  1:5081189/101‑1 (MQ=255)

                             gcTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                    >  2:5081875/1‑101 (MQ=255)

                              cTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCg                                                                   <  1:5081556/101‑1 (MQ=255)

                              cTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTCCg                                                                   <  1:5081317/101‑1 (MQ=255)

                                gATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                 <  1:5081420/101‑1 (MQ=255)

                                gATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                 <  2:5081476/101‑1 (MQ=255)

                                 aTGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGa                                                                >  1:5081882/1‑101 (MQ=255)

                                    cgcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATgg                                                             <  2:5081378/101‑1 (MQ=255)

                                     gcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGt                                                            >  2:5081883/1‑101 (MQ=255)

                                     gcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggggggg                                                                                        >  2:5081886/1‑69 (MQ=255)

                                        cAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                         <  2:5081234/101‑1 (MQ=255)

                                        cAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                         >  2:5081892/1‑101 (MQ=255)

                                         aaCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                         >  2:5081895/1‑100 (MQ=255)

                                          aCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCa                                                       <  2:5081793/101‑1 (MQ=255)

                                          aCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCa                                                       >  1:5081898/1‑101 (MQ=255)

                                            gcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGa                                                     >  1:5081902/1‑101 (MQ=255)

                                             cgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGat                                                    <  1:5081235/101‑1 (MQ=255)

                                             cgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGat                                                    <  1:5081443/101‑1 (MQ=255)

                                               tAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatat                                                  <  1:5080998/101‑1 (MQ=255)

                                                aaCACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatata                                                 >  2:5081906/1‑101 (MQ=255)

                                                 acacAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatatat                                                >  1:5081907/1‑101 (MQ=255)

                                                 acacAGCCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatatat                                                <  1:5081638/101‑1 (MQ=255)

                                                  cacaGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATc                                               <  2:5081576/101‑1 (MQ=255)

                                                   acaGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCt                                              <  1:5081497/101‑1 (MQ=255)

                                                         agtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaataa                                        >  1:5081913/1‑101 (MQ=255)

                                                         agtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaataa                                        <  2:5081678/101‑1 (MQ=255)

                                                         agtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaataa                                        <  2:5081288/101‑1 (MQ=255)

                                                            tCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAAc                                     <  1:5081344/101‑1 (MQ=255)

                                                            tCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAAc                                     >  2:5081915/1‑101 (MQ=255)

                                                              aCGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTc                                   <  1:5081109/101‑1 (MQ=255)

                                                               cGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCg                                  >  2:5081920/1‑101 (MQ=255)

                                                               cGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCg                                  <  1:5081540/101‑1 (MQ=255)

                                                               cGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCg                                  <  2:5081515/101‑1 (MQ=255)

                                                                ggAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc                                 <  2:5081274/101‑1 (MQ=255)

                                                                 gAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa                                <  2:5081073/101‑1 (MQ=255)

                                                                  aaGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaa                               >  2:5081925/1‑101 (MQ=255)

                                                                  aaGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaa                               <  1:5081432/101‑1 (MQ=255)

                                                                  aaGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaa                               >  1:5081924/1‑101 (MQ=255)

                                                                   aGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaaa                              <  1:5081397/101‑1 (MQ=255)

                                                                     cgcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatat                                                  <  1:5081410/79‑1 (MQ=255)

                                                                     cgcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatg                            >  1:5081930/1‑101 (MQ=255)

                                                                     cgcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatg                            <  2:5081628/101‑1 (MQ=255)

                                                                      gcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatga                           >  2:5081932/1‑101 (MQ=255)

                                                                      gcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatga                           <  2:5081467/101‑1 (MQ=255)

                                                                      gcAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatga                           <  1:5081510/101‑1 (MQ=255)

                                                                       cAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTc                                   <  2:5081401/92‑1 (MQ=255)

                                                                       cAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgat                          <  2:5081346/101‑1 (MQ=255)

                                                                        aGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgata                         <  1:5081502/101‑1 (MQ=255)

                                                                         gATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataa                        >  1:5081936/1‑101 (MQ=255)

                                                                          aTCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                       <  2:5081200/101‑1 (MQ=255)

                                                                          aTCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                       <  1:5081452/101‑1 (MQ=255)

                                                                           tCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatg                      >  2:5081940/1‑101 (MQ=255)

                                                                            cATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatga                     <  2:5081610/101‑1 (MQ=255)

                                                                             aTTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatgat                    >  1:5081942/1‑101 (MQ=255)

                                                                              ttGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                   <  2:5081325/101‑1 (MQ=255)

                                                                              ttGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                   >  2:5081945/1‑101 (MQ=255)

                                                                               tGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaataaa                                       <  2:5081561/80‑1 (MQ=255)

                                                                                gCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGt                 <  2:5081488/101‑1 (MQ=255)

                                                                                  gAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTa               >  2:5081948/1‑101 (MQ=255)

                                                                                  gAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTa               <  1:5081460/101‑1 (MQ=255)

                                                                                   aaCGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTaa              <  1:5081384/101‑1 (MQ=255)

                                                                                    aCGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAAt             <  1:5081599/101‑1 (MQ=255)

                                                                                     cgcgCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc            <  1:5081328/101‑1 (MQ=255)

                                                                                     cgcgCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc            <  2:5081506/101‑1 (MQ=255)

                                                                                     cgcgCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc            <  2:5081458/101‑1 (MQ=255)

                                                                                        gcTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg         <  1:5081680/101‑1 (MQ=255)

                                                                                        gcTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg         >  1:5081951/1‑101 (MQ=255)

                                                                                         cTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa        >  1:5081954/1‑101 (MQ=255)

                                                                                          tGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                 <  2:5081425/43‑1 (MQ=37)

                                                                                          tGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAt       <  1:5081573/101‑1 (MQ=255)

                                                                                            cGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa     <  2:5081542/101‑1 (MQ=255)

                                                                                             gggTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaaa    <  1:5081208/101‑1 (MQ=255)

                                                                                             gggTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaaa    <  2:5081334/101‑1 (MQ=255)

                                                                                             gggTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaaa    >  1:5081960/1‑101 (MQ=255)

                                                                                               gTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg  <  1:5081480/101‑1 (MQ=255)

CCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATG  >  NC_000913/3216376‑3216570

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