BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,101,039	(C)_7→8	coding (585/720 nt)	yggN ←	putative EcfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,101,032	1	.	C	99.2%	1095.8 / ‑0.5	269	coding (592/720 nt)	yggN	putative EcfF
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (149/117); total (151/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 302 total aligned reads displayed.

CGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATT  >  NC_000913/3100932‑3101126

cGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑g                                                                                                 >  2:5060716/1‑101 (MQ=255)

  gCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cc                                                                                              >  2:5060721/1‑101 (MQ=255)

  gCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cc                                                                                              >  1:5060720/1‑101 (MQ=255)

    tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cccc                                                                                            <  2:5060239/101‑1 (MQ=255)

      tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cccccc                                                                                          >  2:5060726/1‑101 (MQ=255)

        ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cccccccc                                                                                        <  1:5060077/101‑2 (MQ=255)

        ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cccccccc                                                                                        >  1:5060727/1‑100 (MQ=255)

         tCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑cccccccca                                                                                       <  2:5060151/101‑3 (MQ=255)

           cAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑ccccccccagg                                                                                     >  2:5060730/1‑97 (MQ=255)

            agagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑ccccccccaggc                                                                                    <  2:5060414/101‑6 (MQ=255)

            agagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑ccccccccaggc                                                                                    <  1:5060294/101‑6 (MQ=255)

             gagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCt                                                                                    >  2:5060734/1‑101 (MQ=255)

              agTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCtt                                                                                   <  1:5060272/101‑1 (MQ=255)

                tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTcc                                                                                 >  1:5060738/1‑101 (MQ=255)

                  acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCa                                                                               <  1:5060635/101‑1 (MQ=255)

                   caACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTcc                                                                                 >  2:5060743/1‑98 (MQ=255)

                   caACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAg                                                                              >  1:5060742/1‑101 (MQ=255)

                       gcgGATACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGcaca                                                                          <  2:5060214/101‑1 (MQ=255)

                         ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACAtt                                                                        <  1:5059788/101‑1 (MQ=255)

                         ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACAtt                                                                        >  2:5060749/1‑101 (MQ=255)

                         ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACAtt                                                                        <  1:5060365/101‑1 (MQ=255)

                           cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCt                                                                      >  1:5060750/1‑101 (MQ=255)

                             aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCGCTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGt                                                                    <  1:5060107/101‑1 (MQ=255)

                             aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGt                                                                    >  1:5060754/1‑101 (MQ=255)

                              cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTa                                                                   >  2:5060756/1‑101 (MQ=255)

                               aaaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                  >  2:5060761/1‑101 (MQ=255)

                                   aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                              <  2:5060261/101‑1 (MQ=255)

                                   aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                              >  2:5060765/1‑101 (MQ=255)

                                   aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                              >  1:5060766/1‑101 (MQ=255)

                                    tCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                             <  1:5060304/101‑1 (MQ=255)

                                     cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGtt                                                            <  2:5060336/101‑1 (MQ=255)

                                      tttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTa                                                           >  1:5060771/1‑101 (MQ=255)

                                       ttGTCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                          >  2:5060773/1‑101 (MQ=255)

                                       ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                          <  2:5060491/101‑1 (MQ=255)

                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAcc                                                         <  2:5060278/101‑1 (MQ=255)

                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                          >  1:5060776/1‑100 (MQ=255)

                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                       >  2:5060781/1‑101 (MQ=255)

                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                       >  2:5060784/1‑101 (MQ=255)

                                           cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                      <  1:5060458/101‑1 (MQ=255)

                                           cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                      >  1:5060786/1‑101 (MQ=255)

                                            aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCTCTCGGTTTGGATTGAGGATTGCA‑g                                                                                                 <  1:5059968/57‑1 (MQ=255)

                                            aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                     >  1:5060790/1‑101 (MQ=255)

                                            aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                     >  2:5060791/1‑101 (MQ=255)

                                            aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                     <  2:5060129/101‑1 (MQ=255)

                                                tgctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCttt                                                 >  1:5060796/1‑101 (MQ=255)

                                                tgctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCttt                                                 <  2:5060316/101‑1 (MQ=255)

                                                 gctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtttt                                                >  1:5060799/1‑101 (MQ=255)

                                                 gctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtttt                                                <  1:5060265/101‑1 (MQ=255)

                                                  ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                               >  2:5060802/1‑101 (MQ=255)

                                                   tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                              >  1:5060806/1‑101 (MQ=255)

                                                   tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GACCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                              >  2:5060805/1‑101 (MQ=255)

                                                    ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                             <  2:5059662/101‑1 (MQ=255)

                                                       aaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                          >  2:5060811/1‑101 (MQ=255)

                                                       aaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                          <  2:5060177/101‑1 (MQ=255)

                                                       aaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                          >  1:5060812/1‑101 (MQ=255)

                                                         tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                        >  2:5060816/1‑101 (MQ=255)

                                                           tttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCt                                      >  2:5060820/1‑101 (MQ=255)

                                                           tttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCt                                      <  1:5060160/101‑1 (MQ=255)

                                                            ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCtt                                     >  1:5060821/1‑101 (MQ=255)

                                                             tttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                    <  1:5060468/101‑1 (MQ=255)

                                                              ttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                   <  1:5060222/101‑1 (MQ=255)

                                                               tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                  >  2:5060827/1‑101 (MQ=255)

                                                               tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                  >  1:5060826/1‑101 (MQ=255)

                                                                 cTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCATC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                <  2:5060105/101‑1 (MQ=255)

                                                                 cTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                >  1:5060833/1‑101 (MQ=255)

                                                                 cTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                >  1:5060830/1‑101 (MQ=255)

                                                                   gCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                              >  2:5060837/1‑101 (MQ=255)

                                                                      ttttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                            <  1:5060285/84‑1 (MQ=255)

                                                                       tttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAtt                          <  1:5060499/101‑1 (MQ=255)

                                                                        ttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                         >  1:5060845/1‑101 (MQ=255)

                                                                        ttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                         <  2:5060207/101‑1 (MQ=255)

                                                                        ttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                         <  1:5060204/101‑1 (MQ=255)

                                                                        ttCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                         >  2:5060842/1‑101 (MQ=255)

                                                                          ccACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCa                       >  2:5060850/1‑101 (MQ=255)

                                                                           cACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                      <  2:5060452/101‑1 (MQ=255)

                                                                           cACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                      >  2:5060854/1‑101 (MQ=255)

                                                                           cACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                      >  1:5060851/1‑101 (MQ=255)

                                                                            aCTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAtt                     >  2:5060862/1‑101 (MQ=255)

                                                                            aCTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAtt                     >  1:5060858/1‑101 (MQ=255)

                                                                            aCTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAtt                     <  1:5060388/101‑1 (MQ=255)

                                                                            aCTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAtt                     <  2:5060298/101‑1 (MQ=255)

                                                                             cTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTa                    <  1:5060568/101‑1 (MQ=255)

                                                                             cTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTa                    >  2:5060865/1‑101 (MQ=255)

                                                                               cGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                  >  1:5060866/1‑101 (MQ=255)

                                                                               cGGTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                  <  1:5060033/101‑1 (MQ=255)

                                                                                ggTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                 <  2:5060847/101‑1 (MQ=255)

                                                                                ggTTTGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                 >  2:5060869/1‑101 (MQ=255)

                                                                                  tttGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt               >  2:5060874/1‑101 (MQ=255)

                                                                                   ttGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg              <  1:5060416/101‑1 (MQ=255)

                                                                                   ttGGATTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg              >  1:5060875/1‑101 (MQ=255)

                                                                                       aTTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCt          <  1:5060442/101‑1 (MQ=255)

                                                                                       aTTGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCt          <  2:5060223/101‑1 (MQ=255)

                                                                                        ttGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg         <  2:5060597/101‑1 (MQ=255)

                                                                                         tGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGTGCCCATTTCATTAATGCTGTCCTGt        <  1:5060188/101‑1 (MQ=255)

                                                                                         tGAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGt        <  2:5059987/101‑1 (MQ=255)

                                                                                          gAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa       <  1:5060147/101‑1 (MQ=255)

                                                                                          gAGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa       <  2:5060058/101‑1 (MQ=255)

                                                                                           aGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa      <  1:5060341/101‑1 (MQ=255)

                                                                                           aGGATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa      >  2:5060880/1‑101 (MQ=255)

                                                                                             gATTGCA‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaaaa    >  1:5060883/1‑101 (MQ=255)

CGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATT  >  NC_000913/3100932‑3101126

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