BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,055	Δ1 bp	coding (509/1458 nt)	ascF →	PTS system enzyme IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,055	0	G	.	99.6%	1036.5 / ‑2.5	225	coding (509/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base G (0/1); major base . (106/117); minor base C (0/1); total (106/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 228 total aligned reads displayed.

CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGC  >  NC_000913/2837968‑2838142

 aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                                   <  1:1045407/93‑1 (MQ=255)

  aaaGGATGGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                                  <  2:2616380/93‑1 (MQ=255)

  aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                                  <  1:5571990/93‑1 (MQ=255)

  aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                                  <  1:5362960/93‑1 (MQ=255)

   aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                                 <  2:2584770/93‑1 (MQ=255)

    aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                                <  2:4864468/93‑1 (MQ=255)

    aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                                >  2:3058929/1‑93 (MQ=255)

    aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                                <  1:1312396/93‑1 (MQ=255)

      gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTa                                                                              <  2:1424155/93‑1 (MQ=255)

       aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                                             <  2:3328294/93‑1 (MQ=255)

        tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                                            <  1:3825554/93‑1 (MQ=255)

        tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                                            <  1:2766439/93‑1 (MQ=255)

          gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                                          >  1:883987/1‑93 (MQ=255)

           ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAcc                                                                         >  1:1353645/1‑93 (MQ=255)

            cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                                        >  1:5186403/1‑93 (MQ=255)

            cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                                        >  1:2885142/1‑93 (MQ=255)

             gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                                       <  1:1155288/93‑1 (MQ=255)

              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                                      <  2:727244/93‑1 (MQ=255)

              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                                      <  1:2211618/93‑1 (MQ=255)

              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                                      <  1:1208212/93‑1 (MQ=255)

              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                                      <  2:5119067/93‑1 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     >  2:2342652/1‑93 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     >  2:3644068/1‑93 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     >  2:1241390/1‑93 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   <  2:1878550/93‑1 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   >  2:4783306/1‑93 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   >  2:4658699/1‑93 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   >  1:4578633/1‑93 (MQ=255)

                  tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAATCAATATGt                                                                  >  2:892425/1‑93 (MQ=255)

                   aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                                 <  1:288923/93‑1 (MQ=255)

                   aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                                 <  2:5232944/93‑1 (MQ=255)

                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCTGC‑CATCAAATTTAAAACCAATATGTCGCt                                                              <  2:1208314/93‑1 (MQ=255)

                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                              <  2:4449716/93‑1 (MQ=255)

                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                             <  2:4081254/93‑1 (MQ=255)

                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                             <  2:4241879/93‑1 (MQ=255)

                         tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGc                                                           <  1:2145436/93‑1 (MQ=255)

                           tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                                         >  1:1037356/1‑93 (MQ=255)

                            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                                        <  2:2494952/93‑1 (MQ=255)

                              aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                                      >  2:3388727/1‑93 (MQ=255)

                              aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                                      <  1:6039463/93‑1 (MQ=255)

                               tgtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGc                                                     <  2:2839792/93‑1 (MQ=255)

                               tgtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGc                                                     <  2:402966/93‑1 (MQ=255)

                                  gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                                  <  1:1077787/93‑1 (MQ=255)

                                   aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGTCGATTGCCATc                                                 >  1:3259327/1‑93 (MQ=255)

                                    ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                                >  2:5896242/1‑93 (MQ=255)

                                       gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                                             <  2:4505174/93‑1 (MQ=255)

                                       gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                                             <  1:1763614/93‑1 (MQ=255)

                                         gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                                           <  1:1526481/93‑1 (MQ=255)

                                          gcGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                                          <  2:3338869/92‑1 (MQ=255)

                                            ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGggtg                                        >  1:4809342/1‑93 (MQ=255)

                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGTCATCGCGGGTGTGc                                       >  2:4973586/1‑93 (MQ=255)

                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                                       >  2:4558829/1‑93 (MQ=255)

                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                                       <  1:2632578/93‑1 (MQ=255)

                                              tGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                                      <  2:5294945/93‑1 (MQ=255)

                                               gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                                     >  1:858070/1‑93 (MQ=255)

                                               gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                                     <  2:1451778/93‑1 (MQ=255)

                                                cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                                    >  2:1766596/1‑93 (MQ=255)

                                                cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                                    <  1:1696772/93‑1 (MQ=255)

                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                                   <  1:3385260/93‑1 (MQ=255)

                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                                   >  1:4071262/1‑93 (MQ=255)

                                                  ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                                  >  2:2739515/1‑93 (MQ=255)

                                                    cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                                <  2:2002690/93‑1 (MQ=255)

                                                     ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                               >  1:6115894/1‑93 (MQ=255)

                                                      tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                              >  1:2063724/1‑93 (MQ=255)

                                                       cttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                             >  2:2517782/1‑93 (MQ=255)

                                                          ccTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg                          <  2:2319709/93‑1 (MQ=255)

                                                          ccTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg                          <  1:7372/93‑1 (MQ=255)

                                                          ccTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg                          >  2:426836/1‑93 (MQ=255)

                                                            tGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt                        <  2:6052194/93‑1 (MQ=255)

                                                            tGCCGCTGAAGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt                        >  1:4185611/1‑93 (MQ=255)

                                                             gCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt                       >  1:1829194/1‑93 (MQ=255)

                                                             gCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt                       >  2:3292683/1‑93 (MQ=255)

                                                               cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa                     <  1:5861805/93‑1 (MQ=255)

                                                               cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa                     <  2:2416149/93‑1 (MQ=255)

                                                               cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa                     >  2:5928681/1‑93 (MQ=255)

                                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt                    >  2:37059/1‑93 (MQ=255)

                                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt                    >  2:3899067/1‑93 (MQ=255)

                                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt                    >  2:1002203/1‑93 (MQ=255)

                                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGCGTGCTGGTACATCCGAGTTTTAt                    >  2:5768378/1‑93 (MQ=255)

                                                                 cTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt                   >  2:2799128/1‑93 (MQ=255)

                                                                 cTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt                   <  1:312282/93‑1 (MQ=255)

                                                                 cTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt                   >  2:1503606/1‑93 (MQ=255)

                                                                   gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa                 <  1:5046500/93‑1 (MQ=255)

                                                                   gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa                 >  1:2524560/1‑93 (MQ=255)

                                                                     tGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAc               <  2:3047040/93‑1 (MQ=255)

                                                                     tGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAc               <  1:4888640/93‑1 (MQ=255)

                                                                       gTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTg             >  2:1585587/1‑93 (MQ=255)

                                                                         cgcggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGAt           <  2:2434556/89‑1 (MQ=255)

                                                                          gcggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATg          <  1:4637006/90‑1 (MQ=255)

                                                                           cggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATgg         <  1:5651854/91‑1 (MQ=255)

                                                                            ggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGc        >  2:3540139/2‑93 (MQ=255)

                                                                              cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGa      <  1:2934478/93‑1 (MQ=255)

                                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaaa    >  1:3717867/1‑93 (MQ=255)

                                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaaa    >  1:3749144/1‑93 (MQ=255)

                                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaaa    >  1:4316052/1‑93 (MQ=255)

                                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaaa    >  2:4200021/1‑93 (MQ=255)

                                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATCGCGaaa    >  1:2375877/1‑93 (MQ=255)

                                                                                 cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAg   <  1:3425874/93‑1 (MQ=255)

                                                                                  tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGc  <  2:6138289/93‑1 (MQ=255)

                                                                                  tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGc  <  1:5880220/93‑1 (MQ=255)

CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGC  >  NC_000913/2837968‑2838142

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