BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,987,280	A→G	R78R (CGT→CGC)	yecH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,280	0	A	G	97.6%	526.7 / 0.3	169	R78R (CGT→CGC)	yecH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (69/96); total (70/99)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 179 total aligned reads displayed.

GTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987191‑1987369

  aTTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCa                                                                                      >  2:1441221/1‑93 (MQ=255)

    tttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCaaa                                                                                    >  2:1096967/1‑93 (MQ=255)

    tttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCaaa                                                                                    >  1:1883485/1‑93 (MQ=255)

      tATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATAAGTGGCGGCAAATc                                                                                  >  2:2011940/1‑93 (MQ=255)

        tCGACAAAGACCTCTCTACACTCGACACACAAAACTAATGCTCTCCGCCGACAGGTTTACTTCCGCTGCTTCAAATCAGTGGAGGCCAACCtt                                                                                >  2:2745392/1‑93 (MQ=255)

         cGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTa                                                                               >  2:4682769/1‑93 (MQ=255)

         cGAAAAACACCTCTGTGCCCTCGACACACAAAACTAATGCTCTGCGGCGACAGGTTTACTGTCGCCGGTTCAAATCAGTGGGGGAAAATTTTa                                                                               <  2:1582431/93‑1 (MQ=255)

          gACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTAc                                                                              <  1:1238799/93‑1 (MQ=255)

              aaGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTg                                                                          <  1:2866678/93‑1 (MQ=255)

              aaGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTg                                                                          >  1:841554/1‑93 (MQ=255)

              aaGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTg                                                                          <  1:4570040/93‑1 (MQ=255)

               aGACCTCTCTCCCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGa                                                                         <  2:4302942/93‑1 (MQ=255)

               aGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGa                                                                         <  2:2282276/93‑1 (MQ=255)

               aGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGa                                                                         >  2:3775521/1‑93 (MQ=255)

                gAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGAt                                                                        <  2:3217070/93‑1 (MQ=255)

                gAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGAt                                                                        <  2:2964621/93‑1 (MQ=255)

                 aaCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATc                                                                       >  1:3870298/1‑93 (MQ=255)

                  aCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCg                                                                      >  2:2853499/1‑93 (MQ=255)

                  aCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCg                                                                      >  1:3679916/1‑93 (MQ=255)

                    tctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    <  2:2878040/93‑1 (MQ=255)

                    cctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    <  2:1552783/92‑1 (MQ=255)

                     ctctACCCTCGGCACACAAAACTAATGCTCTCCGGCGGCAAGTTTACTGCCGCCGTTTCGAATCAGTGGCGGCGAATCTTACTCTGCTCGGTc                                                                   >  2:1235808/1‑93 (MQ=255)

                      tctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCg                                                                  >  2:3396417/1‑93 (MQ=255)

                       ctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGa                                                                 <  1:3960787/93‑1 (MQ=255)

                         aCCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaa                                                               >  1:721774/1‑93 (MQ=255)

                          cccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaa                                                              <  2:4056150/93‑1 (MQ=255)

                          cccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaa                                                              >  1:3115387/1‑93 (MQ=255)

                           ccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaaa                                                             <  2:4471380/93‑1 (MQ=255)

                           ccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGCGGCGGCAAATCTTACTCTGATCGGTCGaaaaa                                                             <  1:3242199/93‑1 (MQ=255)

                            cTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAc                                                            <  1:3558533/93‑1 (MQ=255)

                            cTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAc                                                            <  1:4616094/93‑1 (MQ=255)

                            cTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAc                                                            <  1:2430068/93‑1 (MQ=255)

                             tCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAcc                                                           >  2:3990987/1‑93 (MQ=255)

                             tCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAcc                                                           <  1:1564913/93‑1 (MQ=255)

                              cGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCg                                                          >  1:3386963/1‑93 (MQ=255)

                               gACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGt                                                         <  1:2831867/93‑1 (MQ=255)

                               gACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGt                                                         >  1:881387/1‑93 (MQ=255)

                               aacacaaAAACCGAATGCCCCCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGt                                                         <  1:1771866/92‑1 (MQ=255)

                                acacacAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTc                                                        >  2:32650/1‑93 (MQ=255)

                                acacacAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTc                                                        <  1:1291709/93‑1 (MQ=255)

                                 cacacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCt                                                       >  2:4586784/1‑93 (MQ=255)

                                 cacacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCt                                                       >  1:3204315/1‑93 (MQ=255)

                                  acacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCtt                                                      >  1:3887509/1‑93 (MQ=255)

                                  acacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCtt                                                      <  1:1216709/93‑1 (MQ=255)

                                   cacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCttt                                                     >  1:3314848/1‑93 (MQ=255)

                                    acaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCTAAAAACCGTCTTTc                                                    <  2:3192567/93‑1 (MQ=255)

                                    acaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTc                                                    <  1:3525579/93‑1 (MQ=255)

                                       aaaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGaa                                                 <  1:984523/93‑1 (MQ=255)

                                        aaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAg                                                <  2:2335001/93‑1 (MQ=255)

                                        aaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAg                                                <  2:4008899/93‑1 (MQ=255)

                                         aCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAgg                                               <  2:2288459/93‑1 (MQ=255)

                                          cTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGt                                              <  2:4670078/93‑1 (MQ=255)

                                          cTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGt                                              <  2:3516668/93‑1 (MQ=255)

                                           tAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGta                                             >  2:523623/1‑93 (MQ=255)

                                            agtGCTCTCGGGCGGCAGGTTTCCTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtat                                            <  2:17237/91‑1 (MQ=255)

                                            aaTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtat                                            <  1:4713433/93‑1 (MQ=255)

                                             aTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtata                                           >  2:3087779/1‑93 (MQ=255)

                                                 cctcCGGCGGCGGGTTTACTGCCGCCGGTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAtt                                       <  1:4106702/92‑1 (MQ=255)

                                                  ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTTGAAGGTATAAAttt                                      >  2:3291181/1‑93 (MQ=255)

                                                  ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAttt                                      <  2:2643606/93‑1 (MQ=255)

                                                     cggcggCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGcc                                   <  2:1404186/93‑1 (MQ=255)

                                                      ggcggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCt                                  <  1:2853765/93‑1 (MQ=255)

                                                       gcgtcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCtt                                 <  2:656006/93‑1 (MQ=255)

                                                       gcggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCtt                                 <  2:1643205/93‑1 (MQ=255)

                                                          gcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtg                              <  2:3019402/93‑1 (MQ=255)

                                                          gcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtg                              >  2:3186462/1‑93 (MQ=255)

                                                           cAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgc                             <  1:3044414/93‑1 (MQ=255)

                                                           cAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgc                             <  1:420719/93‑1 (MQ=255)

                                                            aGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgct                            <  2:601428/93‑1 (MQ=255)

                                                            aGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgct                            <  2:1317479/93‑1 (MQ=255)

                                                              gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          >  1:2551951/1‑93 (MQ=255)

                                                               tttACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGcc                         <  1:1089403/93‑1 (MQ=255)

                                                                ttACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCa                        >  2:2480233/1‑93 (MQ=255)

                                                                    tGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    <  2:2141963/93‑1 (MQ=255)

                                                                    tGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    <  1:3654417/93‑1 (MQ=255)

                                                                    cgccgccGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    <  2:1409256/92‑1 (MQ=255)

                                                                      ccgccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACg                  <  1:666749/93‑1 (MQ=255)

                                                                      ccgccgTGTCAAATCAGTTGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACg                  <  1:2655505/93‑1 (MQ=255)

                                                                       cgccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGc                 >  1:3775105/1‑93 (MQ=255)

                                                                        gccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCt                <  2:4387251/93‑1 (MQ=255)

                                                                         ccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTa               >  2:990822/1‑93 (MQ=255)

                                                                         ccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTa               <  1:1497768/93‑1 (MQ=255)

                                                                         ccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTa               <  2:1999464/93‑1 (MQ=255)

                                                                          cgTTTCAAATCAGTGGCGTCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAc              <  1:1848123/93‑1 (MQ=255)

                                                                          cgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAc              <  1:2158643/93‑1 (MQ=255)

                                                                           gTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAcc             >  1:2358624/1‑93 (MQ=255)

                                                                             ttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAg           >  1:3399468/1‑93 (MQ=255)

                                                                              tCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGc          <  2:3350371/93‑1 (MQ=255)

                                                                              tCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGc          <  2:950441/93‑1 (MQ=255)

                                                                               cAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGct         <  2:3891507/93‑1 (MQ=255)

                                                                                aaaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctc        >  1:2294580/1‑93 (MQ=255)

                                                                                aaaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctc        <  2:1386540/93‑1 (MQ=255)

                                                                                aaaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctc        <  1:1927619/93‑1 (MQ=255)

                                                                                 aaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctct       >  2:1053987/1‑93 (MQ=255)

                                                                                 aaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctct       >  1:1015196/1‑93 (MQ=255)

                                                                                 aaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctct       >  2:693017/1‑93 (MQ=255)

                                                                                 aaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctct       >  1:1015044/1‑93 (MQ=255)

                                                                                  aTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCtctc      >  1:584502/1‑93 (MQ=255)

                                                                                    cAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTccc    >  2:3917986/1‑93 (MQ=255)

                                                                                      gTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGc  <  2:3632425/93‑1 (MQ=255)

GTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987191‑1987369

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: `‑`