BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	2,587,316	T→C	Y356H (TAT→CAT)	acrD →	aminoglycoside/multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	2,587,316	0	T	C	100.0%	463.3 / NA	117	Y356H (TAT→CAT)	acrD	aminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (73/44); total (73/44)

Only 100 of 117 total aligned reads displayed.

TATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGTATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACTTTAC  >  AC_000091/2587229‑2587402

tATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCAt                                                                                      <  1:1447347/90‑1 (MQ=255)

tATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCAt                                                                                      >  1:2778713/1‑90 (MQ=255)

  tGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCt                                                                                    <  1:2912872/90‑1 (MQ=255)

   gAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTg                                                                                   >  2:1754156/1‑90 (MQ=255)

   gAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTg                                                                                   >  1:3112859/1‑90 (MQ=255)

   gAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTg                                                                                   <  1:548516/90‑1 (MQ=255)

   gAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTg                                                                                   <  2:2571676/90‑1 (MQ=255)

    aaaCCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGt                                                                                  >  1:3388938/1‑90 (MQ=255)

     aaCCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGtt                                                                                 <  1:3147033/90‑1 (MQ=255)

     aaCCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGtt                                                                                 <  2:729677/90‑1 (MQ=255)

      accaccTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTc                                                                                >  2:3465665/1‑90 (MQ=255)

      accaccTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTc                                                                                >  1:521602/1‑90 (MQ=255)

       ccaccTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTcc                                                                               <  1:2703798/90‑1 (MQ=255)

        caccTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCt                                                                              >  2:3220998/1‑90 (MQ=255)

         accTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTg                                                                             >  2:2362998/1‑90 (MQ=255)

          ccTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGc                                                                            <  2:844289/90‑1 (MQ=255)

               tttGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAAc                                                                       <  1:200259/90‑1 (MQ=255)

                ttGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACt                                                                      <  2:2181446/90‑1 (MQ=255)

                   ttAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTcc                                                                   >  1:2280896/1‑90 (MQ=255)

                   ttAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTcc                                                                   <  2:1183095/90‑1 (MQ=255)

                    tAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcg                                                                  <  1:3488781/90‑1 (MQ=255)

                    tAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcg                                                                  >  1:2124157/1‑90 (MQ=255)

                     aaaGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcgc                                                                 >  1:772575/1‑90 (MQ=255)

                     aaaGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcgc                                                                 >  1:3022580/1‑90 (MQ=255)

                      aaGCCTCCATTGAAGGCGTGGTGAAAACGCTGGTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcgcg                                                                >  2:246209/1‑90 (MQ=38)

                      aaGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCcgcg                                                                <  2:508346/90‑1 (MQ=255)

                       aGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCgcgc                                                               <  1:3230981/90‑1 (MQ=255)

                        gCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGGTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGcc                                                              >  2:333220/1‑90 (MQ=255)

                         ccTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCa                                                             >  2:570809/1‑90 (MQ=255)

                          cTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCAc                                                            >  1:154467/1‑90 (MQ=255)

                           tCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCATGCAAAACTTCCGCGCCACg                                                           <  1:248048/90‑1 (MQ=255)

                              aTTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTg                                                        >  1:2358087/1‑90 (MQ=255)

                                tGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGAt                                                      <  2:2306290/90‑1 (MQ=255)

                                 gAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAAGTTCCGCGCCACGCTGATa                                                     >  1:2597920/1‑90 (MQ=255)

                                 gAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATa                                                     >  1:3184557/1‑90 (MQ=255)

                                   aGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATAcc                                                   >  1:1521343/1‑90 (MQ=255)

                                    gACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATAccc                                                  >  2:105369/1‑90 (MQ=255)

                                    gACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATAccc                                                  >  1:97908/1‑90 (MQ=255)

                                      cGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCAc                                                <  2:196576/90‑1 (MQ=255)

                                      cGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCAc                                                >  2:2392499/1‑90 (MQ=255)

                                       gtggtgAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACt                                               >  2:1395242/1‑90 (MQ=255)

                                        tggtgAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACGTCCGCGCCACGCTGATACCCACTa                                              <  1:462525/90‑1 (MQ=255)

                                            gAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGc                                          >  1:758032/1‑90 (MQ=255)

                                             aaaaCGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGcc                                         >  2:3174538/1‑90 (MQ=255)

                                             aaaaCGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGcc                                         <  1:3096579/90‑1 (MQ=255)

                                             aaaaCGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGcc                                         >  2:353468/1‑90 (MQ=255)

                                             aaaaCGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGcc                                         <  2:344881/90‑1 (MQ=255)

                                                aCGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTg                                      <  2:2627464/90‑1 (MQ=255)

                                                    tgctgGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCgg                                  >  1:3402405/1‑90 (MQ=255)

                                                    tgctgGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCgg                                  >  2:412107/1‑90 (MQ=255)

                                                       tgGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCggtgg                               >  2:3404237/1‑90 (MQ=255)

                                                       tgGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCggtgg                               >  1:767267/1‑90 (MQ=255)

                                                       tgGAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCggtgg                               <  1:1724228/90‑1 (MQ=255)

                                                        ggAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCggtggt                              >  1:1205747/1‑90 (MQ=255)

                                                        ggAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCggtggt                              >  1:2133987/1‑90 (MQ=255)

                                                         gAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGgtggtg                             <  2:2646906/90‑1 (MQ=255)

                                                         gAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGgtggtg                             >  2:1894834/1‑90 (MQ=255)

                                                         gAAGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGgcggtg                             >  1:831207/1‑90 (MQ=255)

                                                           aGCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGtt                           >  2:3136943/1‑90 (MQ=255)

                                                            gCTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTg                          <  1:3431486/90‑1 (MQ=255)

                                                             cTATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCCCTATCGCCGTGCCGGTGGTGTTGa                         <  2:1922526/90‑1 (MQ=255)

                                                              tATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGAt                        >  2:1731042/1‑90 (MQ=255)

                                                              tATCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTGGAt                        >  2:897793/1‑90 (MQ=255)

                                                               aTCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATg                       >  2:1007825/1‑90 (MQ=255)

                                                               aTCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATg                       <  2:1108068/90‑1 (MQ=255)

                                                                tCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATgg                      >  1:2363698/1‑90 (MQ=255)

                                                                tCGCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATgg                      >  1:497397/1‑90 (MQ=255)

                                                                  gCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGa                    >  2:2973263/1‑90 (MQ=255)

                                                                  gCTCTGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGa                    >  2:3091432/1‑90 (MQ=255)

                                                                     ctGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAAcc                 >  2:3033127/1‑90 (MQ=255)

                                                                     ctGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAAcc                 >  2:322004/1‑90 (MQ=255)

                                                                      tGGTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCt                <  1:3347506/90‑1 (MQ=255)

                                                                       ggTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCtt               <  2:3015874/90‑1 (MQ=255)

                                                                       ggTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCtt               <  1:2388173/90‑1 (MQ=255)

                                                                       ggTTTTCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGCGGTGTTGATGGGAACCtt               <  2:3018482/90‑1 (MQ=255)

                                                                          tttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTtctc            >  1:716618/1‑90 (MQ=255)

                                                                          tttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTtctc            >  1:2946594/1‑90 (MQ=255)

                                                                          tttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTtctc            <  2:1504633/90‑1 (MQ=255)

                                                                          tttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTtctc            >  1:1957350/1‑90 (MQ=255)

                                                                           ttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTcc           >  2:2869105/1‑90 (MQ=255)

                                                                           ttCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTcc           >  2:3005778/1‑90 (MQ=255)

                                                                           ttCCTCGTTATGCATCTGCTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTcc           <  1:1422340/90‑1 (MQ=255)

                                                                            tCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCg          <  2:2121803/90‑1 (MQ=255)

                                                                            tCCTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCg          >  2:230466/1‑90 (MQ=255)

                                                                             ccTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGt         >  2:1793495/1‑90 (MQ=255)

                                                                             ccTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGt         <  1:2366031/90‑1 (MQ=255)

                                                                              cTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTa        >  1:774574/1‑90 (MQ=255)

                                                                              cTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTa        >  1:1485900/1‑90 (MQ=255)

                                                                              cTCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTa        >  1:1833297/1‑90 (MQ=255)

                                                                               tCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTAc       >  1:2490554/1‑90 (MQ=255)

                                                                               tCGTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTAc       >  1:978906/1‑90 (MQ=255)

                                                                                 gTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACtt     >  1:2797740/1‑90 (MQ=255)

                                                                                 gTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACtt     <  1:515991/90‑1 (MQ=255)

                                                                                 gTTATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGCGTTGATGGGAACCTTCTCCGTACtt     >  2:1335147/1‑90 (MQ=255)

                                                                                  ttATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGCTGGGAACCTTCTCCGTACttt    >  2:2592792/1‑90 (MQ=255)

                                                                                  ttATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACttt    >  2:1144922/1‑90 (MQ=255)

                                                                                  ttATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACttt    <  1:1733976/90‑1 (MQ=255)

                                                                                   tATGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACTTTa   >  2:2674904/1‑90 (MQ=255)

                                                                                    aTGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACTTTAc  <  2:1021728/90‑1 (MQ=255)

                                                                                    aTGCATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACTTTAc  >  2:3151193/1‑90 (MQ=255)

TATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGGAAGCTATCGCTCTGGTTTTCCTCGTTATGTATCTGTTCCTGCAAAACTTCCGCGCCACGCTGATACCCACTATCGCCGTGCCGGTGGTGTTGATGGGAACCTTCTCCGTACTTTAC  >  AC_000091/2587229‑2587402

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