BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,731,175	T→C	R1355R (CGT→CGC)	lhr →	Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,731,175	0	T	C	98.2%	333.3 / ‑2.5	108	R1355R (CGT→CGC)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (54/52); total (55/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 115 total aligned reads displayed.

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

gTCGTTATTTCAGGGTGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATTCTCGACCGCTGCGGCATCATCAGTCGCCa                                                                                         <  1:3670508/93‑1 (MQ=255)

gTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCa                                                                                         >  1:151645/1‑93 (MQ=255)

 tCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAg                                                                                        <  1:4487001/93‑1 (MQ=255)

 tCGTTATTGCAGATGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAg                                                                                        <  1:4487002/93‑1 (MQ=255)

 tCGTAATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAg                                                                                        <  2:4193945/93‑1 (MQ=255)

     tATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCgg                                                                                    >  2:2252677/1‑93 (MQ=255)

      aTTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGt                                                                                   <  2:3772962/93‑1 (MQ=255)

       ttGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTg                                                                                  >  2:1216655/1‑93 (MQ=255)

         gCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGAt                                                                                <  1:1397513/93‑1 (MQ=255)

          cAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATa                                                                               >  1:3336491/1‑93 (MQ=255)

           aGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCTCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAg                                                                              >  2:1434931/1‑93 (MQ=255)

            ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGc                                                                             <  2:1342134/93‑1 (MQ=255)

            ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGc                                                                             <  1:1231115/93‑1 (MQ=255)

             gTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGcc                                                                            >  1:508663/1‑93 (MQ=255)

             gTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGcc                                                                            >  1:508662/1‑93 (MQ=255)

             gTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGcc                                                                            >  1:2699302/1‑93 (MQ=255)

             gTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGcc                                                                            <  1:252005/93‑1 (MQ=255)

             gTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGCGATAGcc                                                                            <  2:229249/93‑1 (MQ=255)

              tGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGTCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCg                                                                           <  2:3136198/93‑1 (MQ=255)

              gggAGCCACTAAGGGATGCTGAATGGATGGTGGCGCGAGCGGAAAATGTGCTCGCCCGCTACGGCATCATCACTCGCCAGGCGGTGATAGCCg                                                                           <  2:4027676/92‑1 (MQ=255)

                gAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGaa                                                                         <  1:4610546/93‑1 (MQ=255)

                 aGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGaaa                                                                        >  2:1140332/1‑93 (MQ=255)

                  gCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGaaaa                                                                       >  1:2643624/1‑93 (MQ=255)

                  gCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGCGACAGCCGaaaa                                                                       >  2:1835226/1‑93 (MQ=255)

                  gCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCCGGGATAGCCGaaaa                                                                       >  1:1085475/1‑93 (MQ=255)

                    cACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGCAAata                                                                     >  1:2951218/1‑93 (MQ=255)

                    cACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAata                                                                     >  2:4548995/1‑93 (MQ=255)

                    cACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAAGCGGTGATAGCCGAAAata                                                                     >  2:2467655/1‑93 (MQ=255)

                      gtagaCGATACCGAAAGGATGCTGGCGGTGGCGGAAAATATGCTCGACCGCTACGGCATCATAAGTCGCCAGGCGGTGATAGCCGAAAATATc                                                                   <  2:339781/89‑1 (MQ=255)

                      cTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATc                                                                   <  1:4708464/93‑1 (MQ=255)

                       tAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATcc                                                                  <  1:1496738/93‑1 (MQ=255)

                       tAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATcc                                                                  >  1:3985169/1‑93 (MQ=255)

                       tAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAACATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATcc                                                                  >  2:1077706/1‑93 (MQ=255)

                        aaaCGATACCGAAAGGATGCTTGCGGTGGCGGAAAATATGCTCGACCGCTAGGGCATCAACCGTCGCCAGGCGGTGATAGCCGAAAATAGccc                                                                 <  2:1700047/93‑1 (MQ=255)

                        aaaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATccc                                                                 <  1:4126944/93‑1 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTCCGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                >  1:1454426/1‑93 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                >  1:1454427/1‑93 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                >  1:4279313/1‑93 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                >  1:1354293/1‑93 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCG‑AAATATCCCTg                                                               <  1:2106865/93‑1 (MQ=255)

                          aCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGCCGGTGATAGCCGAAAATATCCCTg                                                               >  2:4039013/1‑93 (MQ=255)

                          aCGATACCCAAAGGATTCTGGCCCTCGCGGAAAATATGCTCGACCCCTACGGCATCATCAGTCGCCAGGCGCTCATAGCCGAAAATATCCCTg                                                               <  1:3632223/93‑1 (MQ=255)

                           cGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTgg                                                              <  1:4407083/93‑1 (MQ=255)

                            gATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGc                                                             >  1:3545734/1‑93 (MQ=255)

                            gATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGc                                                             >  1:1484178/1‑93 (MQ=255)

                            gATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGc                                                             >  2:2660115/1‑93 (MQ=255)

                             aTACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCg                                                            >  2:19563/1‑93 (MQ=255)

                             aTACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCAAAAATATCCCTGGCg                                                            <  1:1247376/93‑1 (MQ=255)

                                ccGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGt                                                         <  2:2779875/93‑1 (MQ=255)

                                  gAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGttt                                                       >  1:1103063/1‑93 (MQ=255)

                                   aaaGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTc                                                      >  1:445040/1‑93 (MQ=255)

                                    aaGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTcc                                                     >  1:4203158/1‑93 (MQ=255)

                                       gATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATc                                                  >  1:4677687/1‑93 (MQ=255)

                                          gctggcgctggcGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGAt                                               <  1:3693507/93‑1 (MQ=255)

                                               cgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCaaa                                          <  1:3471318/93‑1 (MQ=255)

                                               cgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCaaa                                          <  1:3603838/93‑1 (MQ=255)

                                               cgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCGGGCGGGTTTCCATCGATGCaaa                                          <  2:4128492/93‑1 (MQ=255)

                                               cgctggcgGAAAAGATTCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCaaa                                          <  2:112983/93‑1 (MQ=255)

                                                gctggcgGAAAAAATGCTCGACCGCTACGGCATCATCAGTCGCAAGGCGGTGATAGCCGAAAATATCCCTGGCGGATTTCCATCGATGCAAAc                                         <  2:4210317/93‑1 (MQ=255)

                                                 ctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACg                                        >  1:4027150/1‑93 (MQ=255)

                                                    gcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCtt                                     <  1:4740637/93‑1 (MQ=255)

                                                     cgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCttt                                    >  2:1515613/1‑93 (MQ=255)

                                                      ggAAAATATGCTCGACCGCTATGGCTTCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTTCATCGATGCAAACGCTCTg                                   >  2:4220493/1‑93 (MQ=255)

                                                       gAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGt                                  >  1:2792857/1‑93 (MQ=255)

                                                       gAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGt                                  <  1:2282610/93‑1 (MQ=255)

                                                        aaaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTc                                 >  1:3934696/1‑93 (MQ=255)

                                                        aaaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTc                                 <  2:4255011/93‑1 (MQ=255)

                                                        aaaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTc                                 <  1:2864361/93‑1 (MQ=255)

                                                        aaaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGGTTCCATCGATGCAAACGCTTTGTc                                 >  2:158448/1‑93 (MQ=255)

                                                         aaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCg                                >  1:718593/1‑93 (MQ=255)

                                                          aaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGa                               <  2:446747/93‑1 (MQ=255)

                                                           atatGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGaa                              <  1:371856/93‑1 (MQ=255)

                                                            tatGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAg                             >  2:168831/1‑93 (MQ=255)

                                                              tGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTa                           >  2:1319607/1‑93 (MQ=255)

                                                              tGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTa                           >  1:2067170/1‑93 (MQ=255)

                                                                  cGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGa                       <  2:4627050/93‑1 (MQ=255)

                                                                  cGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGa                       >  1:2283807/1‑93 (MQ=255)

                                                                   gACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGaa                      <  1:2633931/93‑1 (MQ=255)

                                                                    aCCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGTGGGTTTCCATCGATGCAAACGTTTTGTCGAAGTATGGaaa                     >  1:4336456/1‑92 (MQ=255)

                                                                    aCCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGCTTCCATCGATGCAAACGCTTTGTCGACGTATGGAAg                     >  2:3906922/1‑93 (MQ=255)

                                                                      cGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGAc                   <  1:215101/93‑1 (MQ=255)

                                                                      cGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGAc                   <  1:1598081/93‑1 (MQ=255)

                                                                        cTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTc                 <  2:4035399/93‑1 (MQ=255)

                                                                         tACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTcc                <  1:330585/93‑1 (MQ=255)

                                                                         tACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTcc                >  1:744947/1‑93 (MQ=255)

                                                                         tACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTcc                <  2:3885363/93‑1 (MQ=255)

                                                                          aCGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCg               >  2:687789/1‑93 (MQ=255)

                                                                          aCGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCg               <  1:4509791/93‑1 (MQ=255)

                                                                          aCGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCg               >  1:3000788/1‑93 (MQ=255)

                                                                           cGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCgg              <  2:4700050/93‑1 (MQ=255)

                                                                           cGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCgg              <  2:75217/93‑1 (MQ=255)

                                                                             gCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGc            >  1:1181440/1‑93 (MQ=255)

                                                                                   tcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTa      <  2:313052/93‑1 (MQ=255)

                                                                                   tcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTa      >  1:3105722/1‑93 (MQ=255)

                                                                                   tcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTa      <  1:1443759/93‑1 (MQ=255)

                                                                                    caGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGGATGGACTAGTCCGGGCGAATTAt     >  2:1053858/1‑93 (MQ=255)

                                                                                     aGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATg    <  1:545421/93‑1 (MQ=255)

                                                                                     aGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATg    <  1:2815661/93‑1 (MQ=255)

                                                                                       tcgccAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCg  >  2:2109995/1‑93 (MQ=255)

                                                                                       tcgccAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCg  >  1:230575/1‑93 (MQ=255)

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

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