BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,044	+G	coding (498/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,043	1	.	G	99.6%	775.5 / ‑4.0	239	coding (497/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (109/129); total (109/130)
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.33e-01

Only 100 of 272 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑C‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c                                                                                           >  2:356498/1‑93 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑gg                                                                                        >  2:2422372/1‑92 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑gg                                                                                        >  2:2731838/1‑92 (MQ=255)

   gAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑ggc                                                                                       <  1:974538/93‑3 (MQ=255)

        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTg                                                                                   >  2:6551447/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  1:6696733/93‑1 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  >  2:6261925/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  2:6848293/93‑1 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGcc                                                                                 <  2:3001974/93‑1 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcg                                                                                <  2:4003969/93‑1 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcg                                                                                >  2:7786884/1‑93 (MQ=255)

            gCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGTTGGTCG‑CGGCATCTGCcgc                                                                               <  2:2710668/93‑1 (MQ=255)

             cTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgcc                                                                              <  1:3546765/93‑2 (MQ=255)

              tGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgcca                                                                             >  1:8770283/1‑91 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            <  1:8745390/93‑4 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            >  1:1871726/1‑90 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            >  2:6059785/1‑90 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            >  2:3889843/1‑90 (MQ=255)

                  aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCaa                                                                        >  1:6342042/1‑93 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCaaa                                                                       <  1:4655642/93‑1 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCaaa                                                                       >  1:6155410/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAt                                                                      >  1:480364/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAt                                                                      >  2:8543736/1‑93 (MQ=255)

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAttt                                                                    >  2:2375680/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTa                                                                   >  1:5197604/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCC‑CGGCATCTGCCG‑ACATCAAATTTa                                                                   >  1:2681579/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTCCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTa                                                                   >  2:155096/1‑93 (MQ=255)

                          cgccgcCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTaaaa                                                                <  2:3698800/93‑1 (MQ=255)

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                           >  2:7236430/1‑93 (MQ=255)

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                           >  2:1896140/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                          >  1:2568145/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                          <  2:9659050/93‑1 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                          <  2:7908041/93‑1 (MQ=255)

                                  ttAACCCTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                        <  1:7723571/93‑1 (MQ=255)

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                        >  2:101607/1‑93 (MQ=255)

                                    aaCCTTTCTGAATGTGATTGGTGACGATGCTTCCTTCTTCCTGCCGCTGATGCTCG‑CGGCATCTGCCGC‑CATCAAGTTTAAACCCAATATGTc                                                      <  2:8986392/93‑1 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    <  2:3346935/93‑1 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    <  1:1672128/93‑1 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    <  1:3120618/93‑1 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   >  1:273111/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  <  1:9969494/93‑1 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                 >  1:4847721/1‑93 (MQ=255)

                                          tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGc                                                <  1:7074893/93‑1 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                              <  1:2222582/93‑1 (MQ=255)

                                             gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                             <  2:714845/93‑1 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                            >  1:1127465/1‑93 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                            >  1:6723119/1‑93 (MQ=255)

                                               aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                           >  1:9559486/1‑93 (MQ=255)

                                               aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                           >  1:247253/1‑93 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGTG‑GCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  1:10038838/93‑1 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        >  2:1404241/1‑93 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  1:7593095/93‑1 (MQ=255)

                                                  ggATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  1:3198437/92‑1 (MQ=255)

                                                   gATTGGTGACGTTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                       <  2:9831584/93‑1 (MQ=255)

                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                       <  2:2767501/93‑1 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                      <  1:5134800/93‑1 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                      >  1:59712/1‑93 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     <  2:5920524/93‑1 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     >  1:3354348/1‑93 (MQ=255)

                                                      tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                                    <  1:398867/93‑1 (MQ=255)

                                                      tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                                    <  2:2552007/93‑1 (MQ=255)

                                                      tGGTGACGGTGCTTTCTCTTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                                    <  2:4150528/93‑1 (MQ=255)

                                                        gTGACGGTGCTTTCTTCTTCCTGCTGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                                  <  2:4913843/93‑1 (MQ=255)

                                                        gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                                  >  1:6521292/1‑93 (MQ=255)

                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                                <  2:5185054/93‑1 (MQ=255)

                                                           aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                               <  2:2150239/93‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATTGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                              <  1:8157917/93‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                              <  1:8448688/93‑1 (MQ=255)

                                                               tGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                           >  1:8029483/1‑93 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                       >  2:7142863/1‑93 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                       <  2:674620/93‑1 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                       >  2:4432739/1‑93 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                       <  2:1429136/93‑1 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                      >  2:9423764/1‑93 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                      >  2:998742/1‑93 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                      <  1:7042540/93‑1 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  1:2878046/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  1:5002761/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  2:6983862/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                   >  2:1820388/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                   >  1:6904256/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                   >  1:5383420/1‑93 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                <  1:1690918/93‑1 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                >  2:8245031/1‑93 (MQ=255)

                                                                           ccTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg               <  1:7256760/93‑1 (MQ=255)

                                                                             tGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             >  2:9208989/1‑93 (MQ=255)

                                                                             tGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             <  2:7538356/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttttt          >  1:7166350/1‑92 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa          <  2:6963284/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa          <  1:440768/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa          <  1:1512040/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa          <  1:2558257/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGGGTGCTGGTACATCCGAGttttt          >  1:575564/1‑92 (MQ=255)

                                                                                 gCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt         >  2:5865582/1‑93 (MQ=255)

                                                                                  cTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt        >  2:7365498/1‑93 (MQ=255)

                                                                                  cTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt        <  1:3864981/93‑1 (MQ=255)

                                                                                  cTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt        >  2:5476730/1‑93 (MQ=255)

                                                                                  cTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt        <  2:1648642/93‑1 (MQ=255)

                                                                                     aTGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGaa     <  2:459472/93‑1 (MQ=255)

                                                                                     aTGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGACGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGCGTTTTATTgta     >  2:4330186/1‑91 (MQ=255)

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑C‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

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