BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,798	A→G	H270R (CAC→CGC)	mutL →	MutL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,798	0	A	G	99.7%	1073.2 / ‑5.8	316	H270R (CAC→CGC)	mutL	MutL
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (140/175); total (140/176)
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.98e-01

Only 100 of 350 total aligned reads displayed.

TGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGG  >  NC_000913/4395708‑4395887

 ggCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCgcgc                                                                                        <  1:6644567/93‑1 (MQ=255)

  gCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATcgcgcg                                                                                       >  1:4508387/1‑93 (MQ=255)

  gCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATcgcgcg                                                                                       <  2:5940110/93‑1 (MQ=255)

   ccGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGa                                                                                      <  1:7419078/93‑1 (MQ=255)

    cGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGAt                                                                                     <  1:6417750/93‑1 (MQ=255)

     gATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATc                                                                                    <  2:6125590/93‑1 (MQ=255)

     gATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATc                                                                                    <  1:2706461/93‑1 (MQ=255)

      aTCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCTCCTGATCAATCGCGCGATcc                                                                                   <  1:6188755/93‑1 (MQ=255)

      aTCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATcc                                                                                   >  2:1012971/1‑93 (MQ=255)

       tCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCg                                                                                  >  1:4704624/1‑93 (MQ=255)

       tCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCg                                                                                  >  1:2366410/1‑93 (MQ=255)

        ccAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGc                                                                                 >  2:6776402/1‑93 (MQ=255)

         cAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGcc                                                                                <  1:2478584/93‑1 (MQ=255)

         cAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGcc                                                                                <  2:113156/93‑1 (MQ=255)

         cAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGcc                                                                                <  1:55835/93‑1 (MQ=255)

         cAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGcc                                                                                <  2:2756097/93‑1 (MQ=255)

         cAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGcc                                                                                <  1:7058125/93‑1 (MQ=255)

          aaaTCACACCACGCCCCCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCa                                                                               <  1:7650024/93‑1 (MQ=255)

           aaTCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAg                                                                              <  2:3770128/93‑1 (MQ=255)

           aaTCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAg                                                                              >  1:3957763/1‑93 (MQ=255)

           aaTCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAg                                                                              >  2:152719/1‑93 (MQ=255)

            aTCACACCACGCCCGCACTGCCAGAAATTCAGTATTGCCACGTGAACGGCCGCATGATGCGCGATTGCCTGATAAATCGCGCGAAcaggcagg                                                                             >  2:4929432/1‑93 (MQ=255)

             tCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGc                                                                            <  1:2601722/93‑1 (MQ=255)

                cacaacGCCCGCACCGCCAGAAATTCAGTATCGCTACGTGAACGGTCGCATGATGAGCGACCGCCTGATCAATCGCGCGATCCGCCAGGCCTg                                                                         <  2:1964454/93‑1 (MQ=255)

                  ccacGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCg                                                                       >  1:8629927/1‑93 (MQ=255)

                   cacGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGa                                                                      >  1:8155796/1‑93 (MQ=255)

                    acGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGaa                                                                     >  2:6979845/1‑93 (MQ=255)

                    acGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCAGCCAGGCCTGCGaa                                                                     >  2:5051646/1‑93 (MQ=255)

                    acGCCAGCACTGGCAGAAATTAAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGaa                                                                     >  2:3494354/1‑93 (MQ=255)

                      gCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGa                                                                   <  2:6731012/93‑1 (MQ=255)

                       cccGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGAc                                                                  >  1:2102580/1‑93 (MQ=255)

                        ccGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACa                                                                 >  1:2868070/1‑93 (MQ=255)

                         cGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACaa                                                                >  2:3538442/1‑93 (MQ=255)

                          gCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACaaa                                                               >  1:7695711/1‑93 (MQ=255)

                           cACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCTCCAGGCCTGCGAAGACAAAc                                                              >  2:1356283/1‑93 (MQ=255)

                           cACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAAc                                                              <  1:7821360/93‑1 (MQ=255)

                              tGGCAGAGATTCAGTATTGCTGCGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTgg                                                           <  1:4875874/93‑1 (MQ=255)

                              tGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTgg                                                           <  1:7936581/93‑1 (MQ=255)

                              tGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTgg                                                           <  1:3257843/93‑1 (MQ=255)

                               ggCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTggg                                                          <  2:1765775/93‑1 (MQ=255)

                               ggCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTggg                                                          >  2:6796177/1‑93 (MQ=255)

                                gCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTgggg                                                         <  2:5530196/93‑1 (MQ=255)

                                gCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTgggg                                                         <  1:7735435/93‑1 (MQ=255)

                                 cAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTggggg                                                        <  1:1016749/93‑1 (MQ=255)

                                 cAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTggggg                                                        >  1:5867258/1‑93 (MQ=255)

                                 cAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTggggg                                                        <  1:1888146/93‑1 (MQ=255)

                                   gAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGcc                                                      <  1:1310092/93‑1 (MQ=255)

                                   gAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGcc                                                      >  2:8481128/1‑93 (MQ=255)

                                   gAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGcc                                                      >  1:5286868/1‑93 (MQ=255)

                                      aTTCAGTATTGCTACGGGAACGGTCGCATGATGCGCGATCTCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGAt                                                   <  2:584936/93‑1 (MQ=255)

                                       ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATc                                                  >  2:7452620/1‑93 (MQ=255)

                                       ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATc                                                  >  2:2727054/1‑93 (MQ=255)

                                       ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATc                                                  <  2:6578695/93‑1 (MQ=255)

                                       ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATc                                                  >  1:8415844/1‑93 (MQ=255)

                                        tCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCa                                                 <  1:6894036/93‑1 (MQ=255)

                                          aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGc                                               <  1:1573214/93‑1 (MQ=255)

                                             aTTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAAc                                            <  2:1328284/93‑1 (MQ=255)

                                              ttGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAAcc                                           <  2:4677553/93‑1 (MQ=255)

                                              ttGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAAcc                                           <  2:3188225/93‑1 (MQ=255)

                                               tGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCg                                          >  2:7726141/1‑93 (MQ=255)

                                                gCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCgg                                         >  2:3103859/1‑93 (MQ=255)

                                                gCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCgg                                         >  2:3312737/1‑93 (MQ=255)

                                                gCTACGTGAACGGTCGCATGATGCGAGATCGCCTGATCAATCGCGAGATCCGCCAGGCATGCGAAGACAAACTGGGGGCCGATCAGCAACAgg                                         >  2:736889/1‑93 (MQ=255)

                                                   aCGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCAt                                      <  2:7110526/93‑1 (MQ=255)

                                                    cGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCAtt                                     <  2:3005434/93‑1 (MQ=255)

                                                    cGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCAtt                                     <  1:3843903/93‑1 (MQ=255)

                                                       gAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTtgt                                  <  1:8812473/93‑1 (MQ=255)

                                                       gAACGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTtgt                                  <  2:6405731/93‑1 (MQ=255)

                                                          cGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGtgtt                               <  1:5827472/93‑1 (MQ=255)

                                                          cGGTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGtgtt                               <  1:8280274/93‑1 (MQ=255)

                                                           ggTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGtgttg                              >  2:4197838/1‑93 (MQ=255)

                                                           ggTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGtgttg                              >  1:1139621/1‑93 (MQ=255)

                                                            gTCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGtgttgt                             >  1:5148814/1‑93 (MQ=255)

                                                             tCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTa                            <  2:5415399/93‑1 (MQ=255)

                                                             tCGCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTa                            >  1:6683304/1‑93 (MQ=255)

                                                               gCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATc                          <  2:8340007/93‑1 (MQ=255)

                                                               gCATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATc                          >  1:4953807/1‑93 (MQ=255)

                                                                cATGATGCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGGGTTGTATCt                         >  1:8073084/1‑93 (MQ=255)

                                                                 atgatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTg                        <  2:7508495/93‑1 (MQ=255)

                                                                 atgatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTg                        <  2:5225524/93‑1 (MQ=255)

                                                                  tgatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTgg                       <  1:3558067/93‑1 (MQ=255)

                                                                   gatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTTga                      >  1:444650/1‑93 (MQ=255)

                                                                   gatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGga                      <  2:7999193/93‑1 (MQ=255)

                                                                   gatgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGga                      <  1:3007587/93‑1 (MQ=255)

                                                                     tgCGCGATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGgaga                    >  2:1549352/1‑93 (MQ=255)

                                                                       cgcgATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATc                  <  1:1709832/93‑1 (MQ=255)

                                                                        gcgATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCg                 <  2:8619688/93‑1 (MQ=255)

                                                                         cgATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGa                >  2:2316317/1‑93 (MQ=255)

                                                                          gATCGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAc               >  2:2589573/1‑93 (MQ=255)

                                                                             cGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCca            <  1:4232661/93‑1 (MQ=255)

                                                                             cGCCTGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCca            >  2:7309306/1‑93 (MQ=255)

                                                                                 tGATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATc        <  1:4071712/93‑1 (MQ=255)

                                                                                  gATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCa       >  1:2033070/1‑93 (MQ=255)

                                                                                  gATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCa       >  2:4378315/1‑93 (MQ=255)

                                                                                  gATCAATCGCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCa       <  2:2988213/93‑1 (MQ=255)

                                                                                    tcaatcgCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAgg     <  2:1823502/93‑1 (MQ=255)

                                                                                     caatcgCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGt    <  1:3694524/93‑1 (MQ=255)

                                                                                      aatcgCGCGATCCGCCATGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTg   >  2:8168655/1‑93 (MQ=255)

                                                                                      aatcgCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCAGATCAGGTg   >  1:5009412/1‑93 (MQ=255)

                                                                                      aatcgCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTg   <  2:2174125/93‑1 (MQ=255)

TGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGG  >  NC_000913/4395708‑4395887

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