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.1%	336.0 / NA	109	R1355R (CGT→CGC)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (44/64); minor base A (0/1); total (44/65)
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 115 total aligned reads displayed.

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

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

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

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

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

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

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

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

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

    ttATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCg                                                                                     <  1:1190040/93‑1 (MQ=255)

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

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

      aTTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGt                                                                                   >  2:4499458/1‑93 (MQ=255)

      aTTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGt                                                                                   >  1:846540/1‑93 (MQ=255)

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

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

         gCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGAt                                                                                >  2:2122597/1‑93 (MQ=255)

          gaGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATa                                                                               <  1:1648572/92‑1 (MQ=255)

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

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

           aGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAg                                                                              <  1:363723/93‑1 (MQ=255)

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

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

                gAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGaa                                                                         >  1:3645270/1‑93 (MQ=255)

                   ccACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAat                                                                      >  1:854088/1‑93 (MQ=255)

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

                     aCTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAatat                                                                    <  2:3269460/93‑1 (MQ=255)

                     aCTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAatat                                                                    <  1:4420703/93‑1 (MQ=255)

                     aCTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAatat                                                                    <  2:2535671/93‑1 (MQ=255)

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

                      cTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATc                                                                   >  1:2140672/1‑93 (MQ=255)

                       tAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGGATCATCATTCGCCAGGCGGGGATAGCCGAAAATATcc                                                                  >  2:4122018/1‑93 (MQ=255)

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

                        aaaCGATACCGAAAGGATTCTGGCGCTGGCGTAAAATATGCTCGACCGCTACGGCACCATCAGTCGCCAGGCGGTGATAGCCGAAAATAcccc                                                                 <  1:2513781/93‑1 (MQ=255)

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

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

                        aaaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATccc                                                                 >  2:388187/1‑93 (MQ=255)

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

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

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

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

                            gATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGc                                                             <  1:4342202/93‑1 (MQ=255)

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

                            gATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTCGc                                                             >  1:2333486/1‑93 (MQ=255)

                              tACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCgg                                                           <  1:619630/93‑1 (MQ=255)

                               aCCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCggg                                                          <  1:409651/93‑1 (MQ=255)

                               aCCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCggg                                                          <  2:4097242/93‑1 (MQ=255)

                                  gAAAGGATGCTGGCGCTGGCGAAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGttt                                                       <  2:1871491/93‑1 (MQ=255)

                                    aaGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTcc                                                     <  1:375800/93‑1 (MQ=255)

                                     aGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCa                                                    <  2:2745026/93‑1 (MQ=255)

                                     aGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGCTTCCa                                                    >  2:3862354/1‑93 (MQ=255)

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

                                        aTGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCg                                                 <  2:4196807/93‑1 (MQ=255)

                                        aTGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCg                                                 >  1:1609980/1‑93 (MQ=255)

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

                                         tGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGa                                                >  1:1907442/1‑93 (MQ=255)

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

                                           ctggcgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATg                                              >  2:4273133/1‑93 (MQ=255)

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

                                                gctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAAc                                         <  2:3204201/93‑1 (MQ=255)

                                                gctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAAc                                         <  1:4260219/93‑1 (MQ=255)

                                                ccgggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAAc                                         <  2:3817157/90‑1 (MQ=255)

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

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

                                                     gggAAAATATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCttt                                    <  1:3867838/92‑1 (MQ=255)

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

                                                          aaTATGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGa                               >  2:3124947/1‑93 (MQ=255)

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

                                                             atGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGt                            <  1:4441960/93‑1 (MQ=255)

                                                             atGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGt                            <  2:1540984/93‑1 (MQ=255)

                                                             atGCTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGt                            <  2:2835300/93‑1 (MQ=255)

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

                                                                cTCGACCGCTACGGCATCATCAGTCGCCCGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCTACGTATg                         >  2:262956/1‑93 (MQ=255)

                                                                cTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATg                         <  1:1135812/93‑1 (MQ=255)

                                                                cTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATg                         <  1:1060957/93‑1 (MQ=255)

                                                                cTCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATg                         <  2:1232536/93‑1 (MQ=255)

                                                                 tCGACCGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATgg                        >  1:1090928/1‑93 (MQ=255)

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

                                                                     ccGCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGCTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGa                    >  2:4449963/1‑93 (MQ=255)

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

                                                                      cGCTACGCCATCATCATTCGCCAGGCGGAGATAGCAGAAACTATCCCTGGCGGCTTTCCATCTATCCAAACGCTTTGTCGAAGTATGGAAGAc                   <  2:3469422/93‑1 (MQ=255)

                                                                      cGCTACGCCATCATCAGTCGCCAGGCGGTGATAGCTGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGAc                   <  2:1654801/93‑1 (MQ=255)

                                                                       gCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACt                  >  2:789647/1‑93 (MQ=255)

                                                                       gCTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACt                  >  1:3000566/1‑93 (MQ=255)

                                                                        cTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTc                 <  1:2343760/93‑1 (MQ=255)

                                                                        cTACGGCATCATCAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTc                 <  2:4338727/93‑1 (MQ=255)

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

                                                                         tACGGCATCATCAGTCGACAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTcc                <  1:2134548/93‑1 (MQ=255)

                                                                               atcatcAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGa          >  2:2201309/1‑93 (MQ=255)

                                                                               atcatcAGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGa          >  1:3779972/1‑93 (MQ=255)

                                                                                tcatcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGaa         <  1:4489334/93‑1 (MQ=255)

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

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

                                                                                 catcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAt        >  1:706070/1‑93 (MQ=255)

                                                                                  atcaGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGGAtt       >  2:2415108/1‑93 (MQ=255)

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

                                                                                    caGTCGCCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTAt     <  1:73113/93‑1 (MQ=255)

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

                                                                                      gtcgccAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATgg   >  1:265330/1‑92 (MQ=255)

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

                                                                                       tcgccAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCg  <  2:2808517/93‑1 (MQ=255)

GTCGTTATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCG  >  NC_000913/1731085‑1731264

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