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	99.6%	702.2 / NA	221	R1355R (CGT→CGC)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (97/123); minor base A (1/0); total (98/123)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

gTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCGTCAGTCGCCa                                                                                         <  2:2832837/93‑1 (MQ=255)

gTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCa                                                                                         <  2:590975/93‑1 (MQ=255)

gTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCa                                                                                         <  1:5977939/93‑1 (MQ=255)

gTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCa                                                                                         <  2:3686715/93‑1 (MQ=255)

gTCATTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCa                                                                                         <  2:2096639/93‑1 (MQ=255)

 tCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAg                                                                                        >  1:5267084/1‑93 (MQ=255)

  cGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAgg                                                                                       <  1:71785/93‑1 (MQ=255)

   gTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGc                                                                                      >  1:4950822/1‑93 (MQ=255)

   gTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGc                                                                                      >  2:1457492/1‑93 (MQ=255)

   gTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGc                                                                                      <  2:4738032/93‑1 (MQ=255)

    ttATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCg                                                                                     >  2:4027162/1‑93 (MQ=255)

    ttATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCg                                                                                     >  1:5004110/1‑93 (MQ=255)

        tGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGa                                                                                 <  1:1859699/93‑1 (MQ=255)

        tGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGa                                                                                 <  1:4344436/93‑1 (MQ=255)

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

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

         gCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGACAGGCGGTGAt                                                                                >  1:553954/1‑93 (MQ=255)

          cAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATa                                                                               <  2:3166035/93‑1 (MQ=255)

           aGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAg                                                                              >  2:2244862/1‑93 (MQ=255)

            ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGc                                                                             >  2:2497743/1‑93 (MQ=255)

             gTGGAGCCATGAAACTATACCGAAAGGATGCTGGCGCTGGGGGAAAATATGCTCGACCGCTACGGCATCATCAGGCGCCAGGCGGTGATAGcc                                                                            <  2:760533/93‑1 (MQ=255)

              tGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCg                                                                           <  1:2305796/93‑1 (MQ=255)

               ggAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGa                                                                          >  2:1725101/1‑93 (MQ=255)

               ggAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGa                                                                          >  2:5672139/1‑93 (MQ=255)

               ggAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGa                                                                          <  2:4558595/93‑1 (MQ=255)

                gAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGATGATAGCCGaa                                                                         <  1:1522052/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                <  1:6201548/93‑1 (MQ=255)

                         aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCt                                                                <  1:3828023/93‑1 (MQ=255)

                          aCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTg                                                               <  1:2938069/93‑1 (MQ=255)

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

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

                                ccGAAAGTATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGt                                                         <  1:4846112/93‑1 (MQ=255)

                                ccGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGt                                                         >  1:2696865/1‑93 (MQ=255)

                                 cGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGtt                                                        >  1:1302626/1‑93 (MQ=255)

                                 cGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGtt                                                        <  2:5343955/93‑1 (MQ=255)

                                 cGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGtt                                                        <  1:4096144/93‑1 (MQ=255)

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

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

                                     aGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCa                                                    >  1:3682434/1‑93 (MQ=255)

                                      ggATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCAt                                                   <  2:2933532/93‑1 (MQ=255)

                                       gATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATc                                                  <  1:4205731/93‑1 (MQ=255)

                                       gATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATc                                                  <  2:4352838/93‑1 (MQ=255)

                                         tGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGa                                                <  2:651594/93‑1 (MQ=255)

                                          gctggcgctggcGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGAt                                               >  1:1868279/1‑93 (MQ=255)

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

                                            tggcgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGc                                             >  1:5757015/1‑93 (MQ=255)

                                             ggcgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTACATCGATGCa                                            >  1:5683953/1‑93 (MQ=255)

                                             ggcgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTACATCGATGCa                                            >  2:1915725/1‑93 (MQ=255)

                                                gctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAAc                                         >  2:3603555/1‑93 (MQ=255)

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

                                                   ggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCt                                      >  1:4509677/1‑93 (MQ=255)

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

                                                    gcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCtt                                     >  1:3165581/1‑93 (MQ=255)

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

                                                     cgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCttt                                    <  1:4590707/93‑1 (MQ=255)

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

                                                      ggAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTg                                   <  2:5972423/93‑1 (MQ=255)

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

                                                         aaaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCg                                <  1:2364397/93‑1 (MQ=255)

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

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

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

                                                            tatGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCTAAAAGATCCCTGGCGGTTTTCCATCGATGCAAACGCTTTGTCGAAg                             <  2:5250082/93‑1 (MQ=255)

                                                              tGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGGTTCCATCGATGCAAAAGCTTTGTCGAAGTa                           >  2:19675/1‑93 (MQ=255)

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

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

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

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

                                                                   gACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGaa                      >  1:5033124/1‑93 (MQ=255)

                                                                    aCCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAg                     >  1:2225383/1‑93 (MQ=255)

                                                                     ccGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGa                    <  1:5106262/93‑1 (MQ=255)

                                                                     ccGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGa                    <  1:6111291/93‑1 (MQ=255)

                                                                     ccGCTACGGCATCATCAGTCGCCAGGCGGTGATAACCGAAAATATCCCTGGCGGGCTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGa                    >  2:4898168/1‑93 (MQ=255)

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

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

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

                                                                            ggCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCggg             <  1:200655/93‑1 (MQ=255)

                                                                             gCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGc            <  1:1766966/93‑1 (MQ=255)

                                                                               atcatcAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGa          <  1:1090760/93‑1 (MQ=255)

                                                                                tcatcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGaa         >  1:1936173/1‑93 (MQ=255)

                                                                                tcatcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGaa         >  1:1396143/1‑93 (MQ=255)

                                                                                 catcaGTCGCTAGGCGGTGATAGCAGAAACGATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAt        <  1:5310399/93‑1 (MQ=255)

                                                                                 catcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAt        <  2:3496398/93‑1 (MQ=255)

                                                                                 catcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAt        <  2:1006450/93‑1 (MQ=255)

                                                                                  atcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAtt       >  2:2666509/1‑93 (MQ=255)

                                                                                  atcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAtt       >  2:1629924/1‑93 (MQ=255)

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

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

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

                                                                                     aGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATg    >  1:5481016/1‑93 (MQ=255)

                                                                                      gtcgccAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGc   <  1:955696/93‑1 (MQ=255)

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

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