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%	848.7 / ‑4.6	252	coding (497/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (120/131); total (120/132)
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.99e-01

Only 100 of 269 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAC  >  NC_000913/2837951‑2838129

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                         >  1:750587/1‑93 (MQ=255)

 aTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑g                                                                                       >  2:3444510/1‑93 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑gg                                                                                      <  1:6766498/93‑2 (MQ=255)

   gAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggc                                                                                     >  1:3729574/1‑91 (MQ=255)

    aGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggca                                                                                    >  1:1414390/1‑90 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggcat                                                                                   <  2:5003526/93‑5 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggcat                                                                                   >  2:7350647/1‑89 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggcat                                                                                   >  2:3809807/1‑89 (MQ=255)

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                                  <  1:6218125/93‑1 (MQ=255)

        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTg                                                                                 >  1:1317698/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGc                                                                                <  2:588173/93‑1 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGcc                                                                               >  1:2133102/1‑93 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGcc                                                                               <  1:4553990/93‑1 (MQ=255)

              tGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcca                                                                           <  2:2470011/93‑3 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccat                                                                          >  1:4483214/1‑90 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                       <  2:6287378/93‑1 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                       <  2:5163451/93‑1 (MQ=255)

                  aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                      >  1:6349327/1‑93 (MQ=255)

                  aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                      >  2:2021163/1‑93 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                     >  1:2243660/1‑93 (MQ=255)

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

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                    >  1:6437246/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                    >  2:6691963/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                    >  2:6494141/1‑93 (MQ=255)

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

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAttt                                                                  >  1:5056392/1‑93 (MQ=255)

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAttt                                                                  >  2:526545/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTa                                                                 >  2:6929337/1‑93 (MQ=255)

                        aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                                >  2:2802410/1‑93 (MQ=255)

                           gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                             <  1:3894952/93‑1 (MQ=255)

                           gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                             <  1:4001795/93‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                           <  1:5755380/93‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                           <  1:845335/93‑1 (MQ=255)

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

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                         >  2:1370345/1‑93 (MQ=255)

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

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                        <  1:4395821/93‑1 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                        <  1:2853627/93‑1 (MQ=255)

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

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                        <  1:2684429/93‑1 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                       <  1:309453/93‑1 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                       >  1:3608845/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                  <  2:2665726/93‑1 (MQ=255)

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

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                               >  2:3671407/1‑93 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                               >  2:4576454/1‑93 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                            <  1:6676516/93‑1 (MQ=255)

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

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                          >  1:6653675/1‑93 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                          >  1:1126705/1‑93 (MQ=255)

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

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

                                                 gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                                       >  2:4649309/1‑93 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                      <  1:3166517/93‑1 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                    <  1:1467254/93‑1 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                    >  1:732174/1‑93 (MQ=255)

                                                      tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                                  <  2:570306/93‑1 (MQ=255)

                                                        gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                                >  1:4739969/1‑93 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                               >  2:4686696/1‑93 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                               <  1:263828/93‑1 (MQ=255)

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

                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                              >  1:6966924/1‑93 (MQ=255)

                                                           aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                             <  1:1562818/93‑1 (MQ=255)

                                                           aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                             <  1:3982921/93‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                            >  1:5605909/1‑93 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                            <  1:6573141/93‑1 (MQ=255)

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

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

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

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCTGGTGTGCTgg                       >  1:4263114/1‑93 (MQ=255)

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

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

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

                                                                     cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                   >  1:796091/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGTTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                  >  1:4797963/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                 <  1:4983170/93‑1 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg               >  2:4510287/1‑93 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg               <  2:3040020/93‑1 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg               >  1:5521203/1‑93 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa              <  2:546746/93‑1 (MQ=255)

                                                                            cTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt            <  2:6063402/93‑1 (MQ=255)

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

                                                                             tGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt           >  1:5477020/1‑93 (MQ=255)

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

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

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

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

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

                                                                                   tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg     >  2:4218678/1‑93 (MQ=255)

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

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

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

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAC  >  NC_000913/2837951‑2838129

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