BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	2,857,875	G→A	W709* (TGG→TGA)	mutS →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	2,857,875	0	G	A	100.0%	613.1 / NA	156	W709* (TGG→TGA)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (76/80); total (76/80)

Only 100 of 165 total aligned reads displayed.

TTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAA  >  AC_000091/2857788‑2857961

tttACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgagc                                                                                      >  2:2705597/1‑90 (MQ=255)

    cATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTgc                                                                                  <  2:888955/90‑1 (MQ=255)

     aTAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTgcg                                                                                 >  2:175915/1‑90 (MQ=255)

      tAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTgcgc                                                                                >  2:1754325/1‑90 (MQ=255)

       aaCGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGcgcg                                                                               >  2:2774633/1‑90 (MQ=255)

       aaCGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGcgcg                                                                               <  1:1835354/90‑1 (MQ=255)

       aaCGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGcgcg                                                                               <  1:757368/90‑1 (MQ=255)

       aaCGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGcgcg                                                                               <  1:2962479/90‑1 (MQ=255)

        aCGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCgg                                                                              <  2:3014807/90‑1 (MQ=255)

          gCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGaa                                                                            <  2:3277235/90‑1 (MQ=255)

             aCCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAAt                                                                         <  2:1028490/90‑1 (MQ=255)

               cGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCt                                                                       <  2:48632/90‑1 (MQ=255)

               cGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCt                                                                       >  2:2140126/1‑90 (MQ=255)

               cGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCt                                                                       <  1:1736035/90‑1 (MQ=255)

                gAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTg                                                                      <  1:1157413/90‑1 (MQ=255)

                gAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTg                                                                      >  2:2546819/1‑90 (MQ=255)

                 aaTACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTgg                                                                     <  2:1972246/90‑1 (MQ=255)

                  aTACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGc                                                                    >  1:2218154/1‑90 (MQ=255)

                   tACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCg                                                                   <  1:2236195/90‑1 (MQ=255)

                    aCAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGa                                                                  <  1:1704826/90‑1 (MQ=255)

                    aCAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGa                                                                  <  2:2451299/90‑1 (MQ=255)

                    aCAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGa                                                                  <  1:1531234/90‑1 (MQ=255)

                     cAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGaa                                                                 >  1:3415861/1‑90 (MQ=255)

                      aGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAAt                                                                <  1:1134826/90‑1 (MQ=255)

                      aGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAAt                                                                >  2:1924305/1‑90 (MQ=255)

                      aGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAAt                                                                >  2:1762984/1‑90 (MQ=255)

                       gTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATa                                                               <  2:3233213/90‑1 (MQ=255)

                       gTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATa                                                               >  2:1004674/1‑90 (MQ=255)

                        tCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATaa                                                              <  1:1999084/90‑1 (MQ=255)

                         cTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAg                                                             >  2:3143836/1‑90 (MQ=255)

                          tGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGa                                                            <  1:219110/90‑1 (MQ=255)

                            gtgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGAtt                                                          >  2:3164011/1‑90 (MQ=255)

                            gtgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGAtt                                                          <  1:3129860/90‑1 (MQ=255)

                             tgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTa                                                         <  1:2167911/90‑1 (MQ=255)

                             tgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTa                                                         >  1:2373849/1‑90 (MQ=255)

                             tgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTa                                                         <  1:2691035/90‑1 (MQ=255)

                             tgtTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTa                                                         <  1:2973811/90‑1 (MQ=255)

                               ttAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAATATTAAg                                                       <  2:1741446/90‑1 (MQ=255)

                                tAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAgg                                                      >  2:1029656/1‑90 (MQ=255)

                                 aaTGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGc                                                     >  2:568014/1‑90 (MQ=255)

                                 aaTGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGc                                                     >  2:1424750/1‑90 (MQ=255)

                                 aaTGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGc                                                     >  2:3378653/1‑90 (MQ=255)

                                  aTGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                                    <  1:2016797/90‑1 (MQ=255)

                                    ggATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGTCAtt                                                  >  2:852184/1‑90 (MQ=255)

                                    ggATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                                                  <  1:477490/90‑1 (MQ=255)

                                    ggATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                                                  <  1:1344535/90‑1 (MQ=255)

                                    ggATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                                                  >  2:3360213/1‑90 (MQ=255)

                                      aTGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGa                                                >  1:2447742/1‑90 (MQ=255)

                                       tGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGAc                                               <  1:835188/90‑1 (MQ=255)

                                        gagaTCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACg                                              <  2:3383022/90‑1 (MQ=255)

                                         agaTCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGt                                             >  2:1212398/1‑90 (MQ=255)

                                           aTCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTa                                           <  2:63925/90‑1 (MQ=255)

                                           aTCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTa                                           <  2:776255/90‑1 (MQ=255)

                                            tCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAt                                          >  2:84906/1‑90 (MQ=255)

                                            tCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAt                                          <  2:1120701/90‑1 (MQ=255)

                                            tCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAt                                          <  1:490847/90‑1 (MQ=255)

                                            tCGGGCGTGGAACGTCCACCTACGATGGTCCGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAt                                          >  1:476702/1‑90 (MQ=255)

                                              gggCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAttt                                        >  1:1248251/1‑90 (MQ=255)

                                               ggCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTg                                       <  1:2312586/90‑1 (MQ=255)

                                                 cGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCt                                     <  1:498928/90‑1 (MQ=255)

                                                  gTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTa                                    <  2:2077828/90‑1 (MQ=255)

                                                    ggAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                                  >  2:1873729/1‑90 (MQ=255)

                                                     gAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAccc                                 <  2:3023117/90‑1 (MQ=255)

                                                     gAACGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAccc                                 <  2:3272266/90‑1 (MQ=255)

                                                       aCGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc                               >  1:3054116/1‑90 (MQ=255)

                                                        cGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACt                              <  1:1460316/90‑1 (MQ=255)

                                                        cGTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACt                              <  1:343298/90‑1 (MQ=255)

                                                         gTCCACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTa                             >  2:1680969/1‑90 (MQ=255)

                                                           ccACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTAtt                           >  2:513351/1‑90 (MQ=255)

                                                           ccACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTAtt                           >  2:2422540/1‑90 (MQ=255)

                                                            cACCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGCCGTTATTTGCTACCCACTAttt                          <  2:1375863/90‑1 (MQ=255)

                                                             aCCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTc                         <  2:2171322/90‑1 (MQ=255)

                                                             aCCTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTc                         <  2:3437019/90‑1 (MQ=255)

                                                              ccTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCg                        <  2:1791512/90‑1 (MQ=255)

                                                              ccTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCg                        >  1:2895258/1‑90 (MQ=255)

                                                               cTACGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGa                       >  2:1599700/1‑90 (MQ=255)

                                                                 aCGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGc                     <  1:3222889/90‑1 (MQ=255)

                                                                  cGATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCt                    >  2:2303302/1‑90 (MQ=255)

                                                                   gATGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTg                   >  1:445929/1‑90 (MQ=255)

                                                                    aTGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGa                  >  1:2391731/1‑90 (MQ=255)

                                                                    aTGGTCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGa                  >  1:3154879/1‑90 (MQ=255)

                                                                        tCTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCa              <  1:2800764/90‑1 (MQ=255)

                                                                         cTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAg             <  2:1727615/90‑1 (MQ=255)

                                                                         cTGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAg             >  1:3339156/1‑90 (MQ=255)

                                                                          tGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGt            <  1:196368/90‑1 (MQ=255)

                                                                          tGTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGt            >  1:1184924/1‑90 (MQ=255)

                                                                           gTCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGtt           <  1:1934668/90‑1 (MQ=255)

                                                                            tCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTa          >  2:2054059/1‑90 (MQ=255)

                                                                            tCGCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTa          <  2:1748282/90‑1 (MQ=255)

                                                                              gCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTAcc        >  1:1856762/1‑90 (MQ=255)

                                                                              gCTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTAcc        >  1:3367030/1‑90 (MQ=255)

                                                                               cTGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCg       <  2:2180003/90‑1 (MQ=255)

                                                                                tGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCgg      >  1:208766/1‑90 (MQ=255)

                                                                                tGGCGTGAGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCgg      >  2:2708978/1‑90 (MQ=255)

                                                                                 ggcgtgagcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGga     >  2:62573/1‑90 (MQ=255)

                                                                                 ggcgtgagcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGga     >  1:1643500/1‑90 (MQ=255)

                                                                                 ggcgtgagcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGga     >  1:785328/1‑90 (MQ=255)

                                                                                   cgtgagcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGTCGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGgaga   <  2:1519513/90‑1 (MQ=255)

                                                                                   cgtgagcatgCGCGGAAAATCTAGCGAATAAGATAAAGGCAGTGACGTTATTTGCTACCAACTATTTCGAGCTGACCCAGTTACCGgaga   >  2:1246020/1‑90 (MQ=255)

                                                                                    gtgagcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGaa  >  1:3361346/1‑90 (MQ=255)

TTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAA  >  AC_000091/2857788‑2857961

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