BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,019	T→C	V142A (GTG→GCG)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,019	0	T	C	99.6%	812.4 / ‑4.9	261	V142A (GTG→GCG)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (113/147); total (114/147)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 271 total aligned reads displayed.

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

ccATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGgcgc                                                                                         <  2:2275973/93‑1 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        <  2:4991489/93‑1 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        <  1:4258759/93‑1 (MQ=255)

   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTg                                                                                      >  2:5687410/1‑93 (MQ=255)

     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    >  2:1600687/1‑93 (MQ=255)

     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    >  1:1810368/1‑93 (MQ=255)

     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    >  1:3284228/1‑93 (MQ=255)

       tCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCg                                                                                  >  2:7361058/1‑93 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 >  2:3142764/1‑93 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTACGc                                                                                 >  1:355610/1‑93 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGATTTCTGCGAGCAGACTGGCGTACTGACGGGCTGGGAGGCGCTTGTTCgg                                                                                 >  1:1939674/1‑92 (MQ=255)

         tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCt                                                                                <  2:4402665/93‑1 (MQ=255)

             gTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGa                                                                            <  1:3003754/93‑1 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           >  1:4908354/1‑93 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           <  2:5542341/93‑1 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           >  2:7612542/1‑93 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           >  2:692000/1‑93 (MQ=255)

                ggAAAACCACGAATTCAAGCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACGGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  1:76098/93‑1 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  2:6846122/93‑1 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  1:5536606/93‑1 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  1:83268/93‑1 (MQ=255)

                 gAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACAt                                                                        >  1:7379466/1‑93 (MQ=255)

                 gAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACAt                                                                        <  2:5512988/93‑1 (MQ=255)

                  aaaaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATc                                                                       <  2:624000/93‑1 (MQ=255)

                  aaaaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATc                                                                       >  1:5116668/1‑93 (MQ=255)

                     aCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATccgc                                                                    >  1:4165728/1‑91 (MQ=255)

                     aCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcac                                                                    <  1:6164152/93‑1 (MQ=255)

                     aCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcac                                                                    <  1:991435/93‑1 (MQ=255)

                      ccACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcaca                                                                   <  1:7800290/93‑1 (MQ=255)

                      ccACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcaca                                                                   <  1:3816149/93‑1 (MQ=255)

                        aCGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaaa                                                                 <  1:4480498/93‑1 (MQ=255)

                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAAc                                                                <  2:1455461/93‑1 (MQ=255)

                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAAc                                                                <  2:3508670/93‑1 (MQ=255)

                          gAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACg                                                               >  1:700230/1‑93 (MQ=255)

                           aaTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACgg                                                              <  1:201197/93‑1 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             >  2:7292228/1‑93 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             >  1:3505152/1‑93 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             <  1:3073833/93‑1 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACAGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             <  2:2657656/93‑1 (MQ=255)

                              tCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGaa                                                           <  2:5126320/93‑1 (MQ=255)

                                   ccGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                                      >  1:5168402/1‑93 (MQ=255)

                                    cGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                                     <  1:6820271/93‑1 (MQ=255)

                                    cGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                                     >  1:2579301/1‑93 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATTCACAAACGGGAATTATccc                                                    >  1:2429694/1‑93 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    <  2:7213450/93‑1 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    >  1:5793698/1‑93 (MQ=255)

                                      tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCa                                                   <  1:3898157/93‑1 (MQ=255)

                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCg                                                <  1:6877408/93‑1 (MQ=255)

                                          tCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCgg                                               >  1:259563/1‑93 (MQ=255)

                                           ccAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGa                                              <  2:107112/93‑1 (MQ=255)

                                             aaCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATc                                            >  1:7568889/1‑93 (MQ=255)

                                             aaCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATc                                            <  2:5854971/93‑1 (MQ=255)

                                               ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAg                                          >  1:3738777/1‑93 (MQ=255)

                                                cGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGt                                         >  2:6500195/1‑93 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       >  2:3948729/1‑93 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       <  2:931596/93‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       <  1:4980012/93‑1 (MQ=255)

                                                   ttttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTa                                      <  2:3717158/93‑1 (MQ=255)

                                                    tttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAt                                     <  1:1508154/93‑1 (MQ=255)

                                                      tCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTc                                   <  2:5373281/93‑1 (MQ=255)

                                                        tGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCct                                 >  1:4076237/1‑93 (MQ=255)

                                                        tGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCct                                 <  1:5237955/93‑1 (MQ=255)

                                                         gcgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctc                                <  2:6491670/93‑1 (MQ=255)

                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  2:3245614/93‑1 (MQ=255)

                                                          cgcgcAGACGGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  2:5938561/93‑1 (MQ=255)

                                                          cccgcgGACTGGCGTACTGCCGGGCTGTGAGGCGCTTGCCCTCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATACctct                               <  2:7044760/91‑1 (MQ=255)

                                                            cgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTgg                             <  1:3143224/93‑1 (MQ=255)

                                                             gcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGc                            <  1:3428697/93‑1 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                          <  1:4754415/93‑1 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGGTTATTCCTCTGGCgg                          >  2:4452915/1‑93 (MQ=255)

                                                                gACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                         <  1:133008/93‑1 (MQ=255)

                                                                 aCTGGTGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg                        <  2:6256037/93‑1 (MQ=255)

                                                                 aCTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg                        >  2:1231529/1‑93 (MQ=255)

                                                                  cTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGt                       <  2:6443899/93‑1 (MQ=255)

                                                                  cTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGt                       <  2:4884591/93‑1 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      >  2:2890954/1‑93 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      <  1:1187028/93‑1 (MQ=255)

                                                                    ggCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGtaa                     <  1:945563/93‑3 (MQ=255)

                                                                      cGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATc                   <  2:3019281/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  2:6554618/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  2:1988354/93‑1 (MQ=255)

                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGCATTATCCCACCGGATCAGTTTATTCCTCTGGCGGGGACATCCg                 >  1:6036068/1‑93 (MQ=255)

                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCg                 <  1:7142899/93‑1 (MQ=255)

                                                                         aCTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCgg                >  2:808017/1‑93 (MQ=255)

                                                                          cTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGt               <  2:6989932/93‑1 (MQ=255)

                                                                          cTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGt               >  1:1053669/1‑93 (MQ=255)

                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              >  2:2604925/1‑93 (MQ=255)

                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              >  2:2153074/1‑93 (MQ=255)

                                                                            gACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCt             <  2:1746299/93‑1 (MQ=255)

                                                                             aCGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATTCGGTCtt            <  2:7761731/93‑1 (MQ=255)

                                                                              cGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTa           <  1:7417027/93‑1 (MQ=255)

                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          <  2:6034210/93‑1 (MQ=255)

                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         <  2:5770649/93‑1 (MQ=255)

                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         <  1:5735606/93‑1 (MQ=255)

                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        >  1:2194529/1‑93 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      >  1:4041520/1‑93 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      <  1:6367546/93‑1 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:4290689/1‑93 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  1:1537118/1‑93 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  <  2:4237648/93‑1 (MQ=255)

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

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