BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,175,527	A→G	K293R (AAG→AGG)	ycfU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,175,527	0	A	G	97.5%	484.8 / 1.9	157	K293R (AAG→AGG)	ycfU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (66/87); total (68/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

Only 100 of 163 total aligned reads displayed.

AAAAAAATATGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAAGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCGGGA  >  NC_000913/1175437‑1175616

aaaaaaaTATGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGc                                                                                         <  1:3746774/93‑1 (MQ=255)

 aaaaaaTATGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAAGCa                                                                                        >  1:4503376/1‑93 (MQ=255)

     aaTATGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAAGCAgggg                                                                                    >  2:3065979/1‑92 (MQ=255)

       tatGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGa                                                                                  >  2:3599713/1‑93 (MQ=255)

        atgatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGaa                                                                                 <  2:2526173/93‑1 (MQ=255)

        atgatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGaa                                                                                 <  2:252167/93‑1 (MQ=255)

        atgatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGaa                                                                                 >  2:3254321/1‑93 (MQ=255)

          gatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGt                                                                               <  2:1400667/93‑1 (MQ=255)

          gatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGt                                                                               >  1:3422672/1‑93 (MQ=255)

          gatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGt                                                                               <  1:111667/93‑1 (MQ=255)

          gatgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGt                                                                               <  1:1013432/93‑1 (MQ=255)

           atgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTa                                                                              <  1:3365460/93‑1 (MQ=255)

           atgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTa                                                                              <  2:3354409/93‑1 (MQ=255)

           atgGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTa                                                                              <  2:851796/93‑1 (MQ=255)

             gggTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGc                                                                            >  2:618793/1‑93 (MQ=255)

              ggTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCg                                                                           <  1:732069/93‑1 (MQ=255)

                tttACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGAt                                                                         <  2:2057656/93‑1 (MQ=255)

                tttACTGCTGAGCCTGATTGTCGACGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGaa                                                                         <  1:4223328/93‑2 (MQ=255)

                 ttACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATc                                                                        >  1:3119222/1‑93 (MQ=255)

                 ttACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATc                                                                        <  2:4568110/93‑1 (MQ=255)

                    ctgctgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTg                                                                     <  1:3214633/93‑1 (MQ=255)

                    cgtctgAGCCTGATTGTCGCCGTGGCGGCGTTGAGTATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTg                                                                     <  1:3214632/90‑1 (MQ=255)

                     tgctgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGc                                                                    >  2:2519754/1‑93 (MQ=255)

                      tctgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCa                                                                   <  1:3931858/92‑1 (MQ=255)

                      gctgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCa                                                                   >  1:4749302/1‑93 (MQ=255)

                      gctgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCa                                                                   >  1:369404/1‑93 (MQ=255)

                        tgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCaaa                                                                 <  2:2661050/93‑1 (MQ=255)

                        tgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCaaa                                                                 <  2:561614/93‑1 (MQ=255)

                        tgAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCaaa                                                                 >  2:4698338/1‑93 (MQ=255)

                            ccTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCa                                                             >  1:3471843/1‑93 (MQ=255)

                            ccTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCa                                                             >  1:3283407/1‑93 (MQ=255)

                             cTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCaa                                                            <  2:2374367/93‑1 (MQ=255)

                             cTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCaa                                                            <  2:1360587/93‑1 (MQ=255)

                             cTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCaa                                                            <  2:3048511/93‑1 (MQ=255)

                              tGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAg                                                           >  2:4060730/1‑93 (MQ=255)

                              tGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAg                                                           >  1:3974212/1‑93 (MQ=255)

                               gATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAgg                                                          >  2:4160032/1‑93 (MQ=255)

                                  tGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCtt                                                       <  2:3589273/93‑1 (MQ=255)

                                   gTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTa                                                      <  1:4670377/93‑1 (MQ=255)

                                      gCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATAGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACt                                                   <  2:1131754/93‑1 (MQ=255)

                                       ccGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTc                                                  >  1:2724285/1‑93 (MQ=255)

                                       ccGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAAGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTc                                                  <  1:4529619/93‑1 (MQ=255)

                                         gTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcg                                                >  1:1903161/1‑93 (MQ=255)

                                         gTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcg                                                >  2:1739013/1‑93 (MQ=255)

                                          ttGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcgc                                               <  2:1179907/93‑1 (MQ=255)

                                          ttGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcgc                                               <  2:927632/93‑1 (MQ=255)

                                          ttGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcgc                                               <  1:2103712/93‑1 (MQ=255)

                                           tGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCcgcg                                              <  1:1606909/93‑1 (MQ=255)

                                            gcggcgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGa                                             <  1:4181628/93‑1 (MQ=255)

                                            gcggcgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGa                                             >  2:1254689/1‑93 (MQ=255)

                                              ggcgTTTAACATTATTCCCTCACTAGGGCTGATGATAATGGAGAGGCAGGGCGAAGTACCGATCCTGCAAACGCAAGGCTTAACTCCGCGACa                                           <  1:3724481/93‑1 (MQ=255)

                                              ggcgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACa                                           <  2:422765/93‑1 (MQ=255)

                                              ggcgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACa                                           >  1:950821/1‑93 (MQ=255)

                                               gcgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACaa                                          <  1:1937996/93‑1 (MQ=255)

                                                cgTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACaaa                                         <  2:3420329/93‑1 (MQ=255)

                                                 gTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAAt                                        >  2:536611/1‑93 (MQ=255)

                                                 gTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAAt                                        >  2:2677930/1‑93 (MQ=255)

                                                 gTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAAt                                        <  1:4130007/93‑1 (MQ=255)

                                                 gTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAAt                                        >  2:2426578/1‑93 (MQ=255)

                                                  tttAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATc                                       >  1:3259811/1‑93 (MQ=255)

                                                   ttAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCa                                      <  1:1797951/93‑1 (MQ=255)

                                                    tAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCat                                     >  1:1629239/1‑93 (MQ=255)

                                                    tAACATTATTACCTCACTAGGGCTGATGGTAACGGAGAGGCAGGGCGAAGTAGCGATCCCGCAAACGCAAGGCGTAACTCCGCGACAAATCat                                     >  1:4199532/1‑93 (MQ=255)

                                                      aCATTATTACCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCatga                                   >  2:4077362/1‑93 (MQ=255)

                                                         ttattaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATgg                                >  1:1250419/1‑93 (MQ=255)

                                                          tattaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGTTGGt                               >  2:3826344/1‑93 (MQ=255)

                                                          tattaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGt                               >  2:3826345/1‑93 (MQ=255)

                                                           attaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTc                              >  2:4316981/1‑93 (MQ=255)

                                                           attaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTc                              <  1:689679/93‑1 (MQ=255)

                                                            ttaCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCt                             <  1:2936425/93‑1 (MQ=255)

                                                             taCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCtt                            <  2:3149757/93‑1 (MQ=255)

                                                              aCCTCGCTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCttt                           <  1:807128/93‑1 (MQ=255)

                                                              aCCTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCttt                           >  2:2250716/1‑93 (MQ=255)

                                                               ccTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTa                          <  2:3637432/93‑1 (MQ=255)

                                                                cTCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTAt                         >  2:1141153/1‑93 (MQ=255)

                                                                 tCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATg                        >  2:1876697/1‑93 (MQ=255)

                                                                 tCACTAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATg                        >  2:3129160/1‑93 (MQ=255)

                                                                     tAGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGc                    <  1:3437605/93‑1 (MQ=255)

                                                                      aGGGCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCa                   >  1:13260/1‑93 (MQ=255)

                                                                         gCTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAgg                <  2:3402129/93‑1 (MQ=255)

                                                                          cTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAggg               <  1:4301683/93‑1 (MQ=255)

                                                                          cTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAggg               <  1:2876173/93‑1 (MQ=255)

                                                                          cTGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAggg               <  1:2165707/93‑1 (MQ=255)

                                                                           tGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAgggg              <  2:1972525/93‑1 (MQ=255)

                                                                           tGATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAgggg              >  1:4548996/1‑93 (MQ=255)

                                                                            gATGGTATTGTAGAGGCAGTGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGc             <  1:2563006/93‑1 (MQ=255)

                                                                            gATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGc             <  1:2563005/93‑1 (MQ=255)

                                                                            gATGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGc             >  2:4547257/1‑93 (MQ=255)

                                                                             aTGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGcc            <  1:212155/93‑1 (MQ=255)

                                                                             aTGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGcc            <  2:1817791/93‑1 (MQ=255)

                                                                              tGGTAATGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCa           <  1:2269414/93‑1 (MQ=255)

                                                                                 tAACGGAGAGGCGGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCTTGATGGTCTTTATGGTGCAAGGGGCCAgcg        <  1:4246811/93‑1 (MQ=255)

                                                                                  aaTGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAgcgc       >  1:2493131/1‑93 (MQ=255)

                                                                                  aaTGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAgcgc       >  2:1621439/1‑93 (MQ=255)

                                                                                   aTGGAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGcc      <  1:2953661/93‑1 (MQ=255)

                                                                                     ggAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCgg    <  1:272346/93‑1 (MQ=255)

                                                                                     ggAGAGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCgg    <  2:651292/93‑1 (MQ=255)

                                                                                      gagaGGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCggg   >  1:1013414/1‑93 (MQ=255)

                                                                                      gagaAGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCggg   >  2:4255895/1‑93 (MQ=255)

                                                                                       agagGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCGGGa  <  1:2244543/93‑1 (MQ=255)

AAAAAAATATGATGGGTTTACTGCTGAGCCTGATTGTCGCCGTTGCGGCGTTTAACATTATTACCTCACTAGGGCTGATGGTAATGGAGAAGCAGGGCGAAGTAGCGATCCTGCAAACGCAAGGCTTAACTCCGCGACAAATCATGATGGTCTTTATGGTGCAAGGGGCCAGCGCCGGGA  >  NC_000913/1175437‑1175616

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