BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,141	T→C	S166P (TCA→CCA)	b1545 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,632,141	0	T	C	100.0%	408.5 / NA	140	S166P (TCA→CCA)	b1545	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (63/77); total (63/77)

Only 100 of 142 total aligned reads displayed.

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAA  >  NC_000913/1632051‑1632213

aTAGTAAGAGCCCGCGGCGCAGATAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCa                                                                        <  1:4008197/93‑1 (MQ=255)

  aGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAgg                                                                      <  1:2963852/93‑1 (MQ=255)

  aGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAgg                                                                      >  1:2420445/1‑93 (MQ=37)

   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAggg                                                                     >  1:5915152/1‑92 (MQ=255)

   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                     <  1:4108464/93‑1 (MQ=255)

   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                     >  1:2666838/1‑93 (MQ=37)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                    >  2:3614425/1‑93 (MQ=37)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                    >  1:1610621/1‑93 (MQ=37)

     aaGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtgt                                                                   <  2:4068914/93‑1 (MQ=255)

     aaGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtgt                                                                   <  1:2210715/93‑1 (MQ=255)

      agagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTa                                                                  <  1:6079815/93‑1 (MQ=37)

      agagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTa                                                                  >  2:3317532/1‑93 (MQ=37)

       gagCCCGCGGCGCAGTGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTaa                                                                 <  1:3615328/93‑1 (MQ=255)

       gagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAACCCAGGTGTaa                                                                 >  1:2520187/1‑93 (MQ=255)

         gCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGt                                                               >  1:4704499/1‑93 (MQ=37)

         gCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGt                                                               <  2:5854373/93‑1 (MQ=255)

          cccGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGta                                                              >  2:1217520/1‑93 (MQ=37)

           ccGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtat                                                             <  1:2120154/93‑1 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                            >  1:3011704/1‑93 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                            >  1:4129723/1‑93 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                            <  1:1024016/93‑1 (MQ=255)

              gggcgCAGTGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                          <  1:143600/92‑1 (MQ=255)

              cggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                          <  2:4183357/93‑1 (MQ=255)

              cggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                          <  1:4609742/93‑1 (MQ=255)

              cggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                          <  1:651569/93‑1 (MQ=255)

               ggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGt                                                         <  1:5206668/93‑1 (MQ=255)

               ggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGt                                                         <  1:1150480/93‑1 (MQ=255)

               ggcgCAGGGAAACGTTGTGGTCGACCACCTGTGTAAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGt                                                         <  2:1446102/93‑1 (MQ=255)

                gcgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTg                                                        <  2:677126/93‑1 (MQ=255)

                gcgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTg                                                        >  1:2239874/1‑93 (MQ=255)

                 cgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGc                                                       >  1:1897929/1‑93 (MQ=255)

                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                     <  1:977789/93‑1 (MQ=255)

                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                     <  1:2829913/93‑1 (MQ=255)

                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                     <  1:2384355/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                    <  2:2852493/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                    <  2:5315379/93‑1 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGTTGTAAGTATACGTGCCAtt                                                   <  2:6152707/93‑1 (MQ=11)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                   <  2:2568643/93‑1 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                   <  2:2407313/93‑1 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                   >  2:1436831/1‑93 (MQ=255)

                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                  >  1:3651944/1‑93 (MQ=255)

                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                  >  2:2563970/1‑93 (MQ=255)

                       gAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGc                                                 <  1:5343510/93‑1 (MQ=255)

                       gAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGc                                                 >  2:3721662/1‑93 (MQ=255)

                       gAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGc                                                 <  2:2111460/93‑1 (MQ=255)

                        aaaCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGcc                                                >  2:4644504/1‑93 (MQ=255)

                          aCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTTGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCg                                              >  2:1015707/1‑93 (MQ=255)

                           cGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCgg                                             >  1:199712/1‑93 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGa                                           >  2:5439490/1‑93 (MQ=255)

                              tttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGaa                                          <  2:1372279/93‑1 (MQ=255)

                                 ggTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTc                                       <  1:5937365/93‑1 (MQ=255)

                                 ggTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTc                                       <  2:2793223/93‑1 (MQ=37)

                                   tCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaa                                     >  2:3618728/1‑93 (MQ=255)

                                    cGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaa                                    <  1:816123/93‑1 (MQ=37)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                   >  1:5890108/1‑93 (MQ=255)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                   >  2:2825527/1‑93 (MQ=37)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                   <  1:1893136/93‑1 (MQ=37)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                   >  2:1211011/1‑93 (MQ=255)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                   <  1:2695356/93‑1 (MQ=255)

                                      acgcccTGTGTTAAATGAAGAACAGACACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                  <  1:155973/89‑1 (MQ=11)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                  <  1:4247651/93‑1 (MQ=255)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                  >  1:4364531/1‑93 (MQ=255)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                  >  2:6286875/1‑93 (MQ=255)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                  <  2:3532350/93‑1 (MQ=255)

                                       ccaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAcc                                 >  2:1639870/1‑93 (MQ=255)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                               >  2:6172540/1‑93 (MQ=25)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                               <  1:3099057/93‑1 (MQ=21)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                               <  2:5806784/93‑1 (MQ=21)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                               <  2:4555286/93‑1 (MQ=21)

                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                             <  1:3888314/93‑1 (MQ=18)

                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGGATAAGTGCCATTGCCCGGGAATTCAAAACGTcgc                             >  1:86078/1‑93 (MQ=255)

                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgcg                            >  2:5659825/1‑93 (MQ=21)

                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgcg                            <  1:93152/93‑1 (MQ=18)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                           >  1:1089910/1‑93 (MQ=21)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                           >  1:5457632/1‑93 (MQ=21)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                           >  2:4065662/1‑93 (MQ=18)

                                                 tAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCaaa                       <  2:3733866/93‑1 (MQ=18)

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                      <  2:2668754/93‑1 (MQ=18)

                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                    <  2:2914827/93‑1 (MQ=18)

                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                    <  2:4735294/93‑1 (MQ=18)

                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                    <  2:2196643/93‑1 (MQ=18)

                                                     tGAAGAACAGAAACAGGCGGTATTCGTACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAt                   <  2:4504847/93‑1 (MQ=16)

                                                      gaagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtt                  >  2:2608863/1‑93 (MQ=21)

                                                       aagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAttt                 <  2:5873621/93‑1 (MQ=18)

                                                       aagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAttt                 <  1:3134537/93‑1 (MQ=18)

                                                        agaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtttt                <  2:365749/93‑1 (MQ=18)

                                                          aaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAACTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTaa              <  2:1570987/93‑1 (MQ=16)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag             >  2:5201351/1‑93 (MQ=21)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag             <  1:2967922/93‑1 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga            >  1:5503664/1‑93 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga            <  2:5598192/93‑1 (MQ=18)

                                                             aGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAagag           >  1:2211937/1‑93 (MQ=21)

                                                              gAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGc          >  1:4230566/1‑93 (MQ=18)

                                                               aaaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGcc         >  2:5107255/1‑93 (MQ=21)

                                                                aaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCa        >  1:4756313/1‑93 (MQ=21)

                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa       <  1:1429535/93‑1 (MQ=18)

                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa      >  2:5343903/1‑93 (MQ=18)

                                                                   aGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAg     >  1:1459695/1‑93 (MQ=18)

                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAGGc    <  1:2889978/93‑1 (MQ=16)

                                                                     gCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCa   >  2:3454748/1‑93 (MQ=17)

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAA  >  NC_000913/1632051‑1632213

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