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	99.5%	677.2 / NA	214	R78R (CGT→CGC)	yecH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (114/99); minor base T (1/0); total (115/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

TGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987190‑1987369

  tATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCCGCGGCAGGGTTACTGCCGCCGGTTCAAATCCGTggcggc                                                                                       >  2:4752885/1‑93 (MQ=255)

    ttttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCaa                                                                                     >  2:1490278/1‑93 (MQ=255)

    ttttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCaa                                                                                     >  1:3372219/1‑93 (MQ=255)

       tATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAAtt                                                                                  >  1:961808/1‑92 (MQ=255)

         tCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCtt                                                                                >  2:2465692/1‑93 (MQ=255)

          caaCAAAGACCTCCCTCCCCGCGACACACAAACCTAAGGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTa                                                                               <  2:3069938/91‑1 (MQ=255)

          cGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTa                                                                               <  2:3322766/93‑1 (MQ=255)

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

              aaaGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTActct                                                                           <  2:2761710/93‑1 (MQ=255)

              aaaGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTActct                                                                           <  1:7875765/93‑1 (MQ=255)

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

                 gAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGAt                                                                        >  2:2559889/1‑93 (MQ=255)

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

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

                  aaCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATc                                                                       <  1:7268297/93‑1 (MQ=255)

                   aCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTgctag                                                                      >  2:5862857/1‑89 (MQ=255)

                    ctctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCgg                                                                     <  2:2558259/93‑1 (MQ=255)

                     tctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    >  2:7203357/1‑93 (MQ=255)

                        ctACCCTCGACACACAAAACTAATGCTCCCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGa                                                                 <  1:3582892/93‑1 (MQ=255)

                         tACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaa                                                                >  1:4247324/1‑93 (MQ=255)

                         tACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaa                                                                >  1:633482/1‑93 (MQ=255)

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

                            ccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaaa                                                             >  1:3599928/1‑93 (MQ=255)

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

                             cTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAc                                                            >  1:6155321/1‑93 (MQ=255)

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

                               cGACACACAAAACTAAGCCTCCCCGGCGGCAGGTTTCCTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCg                                                          <  2:1654538/93‑1 (MQ=255)

                               cGACACAAAAAACAAAGGCTCTCCGGCGGCGGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCg                                                          <  1:1436747/93‑1 (MQ=255)

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

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

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

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

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

                                 acacacAAAACTAATGCTCGCCGACGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGGTCGGTCGAAAAACCGTc                                                        <  2:5405150/93‑1 (MQ=255)

                                 acacacAAAACTAATGCTATCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTc                                                        <  2:824178/93‑1 (MQ=255)

                                  cacacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCt                                                       <  2:7711129/93‑1 (MQ=255)

                                  cacacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCt                                                       <  1:3714790/93‑1 (MQ=255)

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

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

                                     acaAAACTAGTGCTCTCGGGCGGCAGTTTTACGGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTc                                                    <  2:2880077/93‑1 (MQ=255)

                                     acaAAACTAATGCTCTCCGGCGGCATGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTc                                                    <  2:7368466/93‑1 (MQ=255)

                                     acaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTc                                                    >  1:6454261/1‑93 (MQ=255)

                                      caAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCg                                                   <  2:6602440/93‑1 (MQ=255)

                                      caAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAACGTCTTTCg                                                   >  2:4066354/1‑93 (MQ=255)

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

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

                                        aaaCTAATGCTCTCCGGAGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGaa                                                 >  1:3203367/1‑93 (MQ=255)

                                         aaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAg                                                >  2:7839586/1‑93 (MQ=255)

                                           cTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGt                                              >  1:2844793/1‑93 (MQ=255)

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

                                            tAATGCTCGCCGGCGGCAGGTTTACTGCTGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGta                                             <  1:1870594/93‑1 (MQ=255)

                                             aaTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtat                                            >  2:208370/1‑93 (MQ=255)

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

                                             aaTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtat                                            >  2:2310765/1‑93 (MQ=255)

                                              aTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTCCGAAGGtata                                           <  2:5092381/93‑1 (MQ=255)

                                               tGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATaa                                          >  2:2045809/1‑93 (MQ=255)

                                                gCTCTCCGGCGGCAGGTTTACTGCCGCCGGTTCAAATCAGTGGCGGCAAATCTTACTCTGCTCGGTCGAAAAACCGTCTTTCGGAGGGATaat                                         >  2:3166295/1‑92 (MQ=255)

                                                gCTCTCCGGCGGCAGGTTTACTCCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATaaa                                         <  1:4230214/93‑1 (MQ=255)

                                                 ctctCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAt                                        <  1:233883/93‑1 (MQ=255)

                                                 cccgccGGCGCAAGTTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAt                                        <  2:3233991/89‑1 (MQ=255)

                                                   ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAttt                                      >  2:5450418/1‑93 (MQ=255)

                                                   ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAttt                                      >  2:1927541/1‑93 (MQ=255)

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

                                                    tcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTg                                     >  2:5740119/1‑93 (MQ=255)

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

                                                        gcggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCtt                                 >  1:6949652/1‑93 (MQ=255)

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

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

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

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

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

                                                             aGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgct                            >  1:5235393/1‑93 (MQ=255)

                                                              ggTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATACATTTGCCTTTtgctc                           >  1:4127735/1‑92 (MQ=255)

                                                              ggTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctc                           >  1:4069731/1‑92 (MQ=255)

                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          <  2:2629360/93‑1 (MQ=255)

                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          <  2:7519013/93‑1 (MQ=255)

                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          <  2:7557802/93‑1 (MQ=255)

                                                                    cTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaa                     <  1:2403896/93‑1 (MQ=255)

                                                                     tGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAACACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    >  2:7392004/1‑93 (MQ=255)

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

                                                                     tGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    >  1:4420002/1‑93 (MQ=255)

                                                                     tGCCGCCGTGTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAGGGTATAAATTTGCCTTTTGCTGCCAGaaa                    <  2:7880403/93‑1 (MQ=255)

                                                                      gccgccGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAAc                   >  2:1774569/1‑93 (MQ=255)

                                                                      gccgccGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAAc                   <  2:2231612/93‑1 (MQ=255)

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

                                                                         gccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCCGAAACGCt                >  1:2969182/1‑93 (MQ=255)

                                                                            gTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGTCTTTTGCTGCCAGAAACGCTAcc             >  2:4426602/1‑93 (MQ=255)

                                                                            gTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAcc             <  1:7429252/93‑1 (MQ=255)

                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCa            >  1:5856182/1‑93 (MQ=255)

                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAACGTCTTTCGAAGGGATAAATTTGCCTTTTGCTGCCAGAAACGCTACCa            >  2:6315060/1‑93 (MQ=255)

                                                                              ttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAg           <  2:7152064/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                                    tCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTcc     <  2:213408/93‑1 (MQ=255)

                                                                                    tCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTcc     >  1:1959335/1‑93 (MQ=255)

                                                                                    tCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTcc     <  2:4052480/93‑1 (MQ=255)

TGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987190‑1987369

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