BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,140	A→G	I72V (ATT→GTT)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,140	0	A	G	98.8%	839.6 / ‑1.2	259	I72V (ATT→GTT)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (123/133); total (125/134)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 262 total aligned reads displayed.

GGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTA  >  NC_000913/4049052‑4049229

gggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTcc                                                                                       >  1:799492/1‑93 (MQ=255)

gggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTcc                                                                                       >  1:4491062/1‑93 (MQ=255)

gggCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTcc                                                                                       <  1:336201/93‑1 (MQ=255)

   cAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCAtt                                                                                    <  1:588664/93‑1 (MQ=255)

   cAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCAtt                                                                                    <  2:2626122/93‑1 (MQ=255)

     gCCGCGCAGTGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTga                                                                                  <  2:3338480/93‑2 (MQ=255)

     gCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCCACAGATCCACGTATTGGCAGTACAACCCCTGTTCCATTgc                                                                                  >  1:1978402/1‑92 (MQ=255)

     gCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTga                                                                                  <  1:1849321/93‑2 (MQ=255)

         cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGACCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                              <  1:485020/93‑1 (MQ=255)

            aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                           <  1:2755814/93‑1 (MQ=255)

             gcgggcggAAACACCACGTCAGGCAGCAAAGGCCTGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                          <  1:3809782/93‑1 (MQ=255)

               gtgcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                        <  2:3906839/91‑1 (MQ=255)

                ggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaa                                                                       >  1:3939953/1‑93 (MQ=255)

                ggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaa                                                                       >  1:977739/1‑93 (MQ=255)

                 gcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaa                                                                      <  2:1280244/93‑1 (MQ=255)

                  cggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaa                                                                     <  1:2990866/93‑1 (MQ=255)

                   ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    <  2:1625903/93‑1 (MQ=255)

                   ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                    >  2:135640/1‑93 (MQ=255)

                     aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGt                                                                  >  2:147303/1‑93 (MQ=255)

                      aaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTc                                                                 >  2:1214136/1‑93 (MQ=255)

                      aaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTc                                                                 >  1:363786/1‑93 (MQ=255)

                         accacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCAcc                                                              <  2:4094470/93‑1 (MQ=255)

                          ccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCa                                                             <  2:3946388/93‑1 (MQ=255)

                             cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                          <  2:1867119/93‑1 (MQ=255)

                              gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                         <  2:2747781/93‑1 (MQ=255)

                                cAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGc                                                       >  1:4047945/1‑93 (MQ=255)

                                  ggCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGcac                                                     >  1:306110/1‑93 (MQ=255)

                                     agcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaaa                                                  <  2:2506398/93‑1 (MQ=255)

                                        aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                               >  1:4577391/1‑93 (MQ=255)

                                        aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                               >  2:1767631/1‑93 (MQ=255)

                                        aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                               <  1:1652170/93‑1 (MQ=255)

                                         aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                              >  1:2380554/1‑93 (MQ=255)

                                         aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                              >  2:1060182/1‑93 (MQ=255)

                                          aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                             <  2:1951251/93‑1 (MQ=255)

                                            gCCAGAACGCACCTACAGATCCGCGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                           <  1:4225682/93‑1 (MQ=255)

                                              cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                         >  1:3158329/1‑93 (MQ=255)

                                              cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                         <  1:1587033/93‑1 (MQ=255)

                                              cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                         <  1:3552091/93‑1 (MQ=255)

                                               aGAACGCACCAAAAGATCCGCGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAagagcg                                        >  1:1728400/1‑91 (MQ=255)

                                               aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                        <  1:1270776/93‑1 (MQ=255)

                                               aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                        <  1:1924991/93‑1 (MQ=255)

                                                gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCATCAAACAGCACAAACCGAAGAGTga                                       >  2:4338148/1‑93 (MQ=255)

                                                gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                       >  1:403883/1‑93 (MQ=255)

                                                gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                       >  2:3804672/1‑93 (MQ=255)

                                                gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                       <  1:2312438/93‑1 (MQ=255)

                                                 aaCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgag                                      >  2:3871745/1‑93 (MQ=255)

                                                 aaCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgag                                      >  1:2018786/1‑93 (MQ=255)

                                                    gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                   >  1:1113699/1‑93 (MQ=255)

                                                      aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  1:2255014/1‑93 (MQ=255)

                                                      aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  2:3413158/1‑93 (MQ=255)

                                                      aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 >  1:2502773/1‑93 (MQ=255)

                                                      aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                 <  2:3000215/93‑1 (MQ=255)

                                                       ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                >  1:1378540/1‑93 (MQ=255)

                                                       ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                <  1:1224256/93‑1 (MQ=255)

                                                       ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                <  2:499894/93‑1 (MQ=255)

                                                       ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                >  1:3686985/1‑93 (MQ=255)

                                                        cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                               >  1:3776505/1‑93 (MQ=255)

                                                          aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTTAGCGTGGCTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                             >  2:682474/1‑93 (MQ=255)

                                                          aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                             <  1:2828535/93‑1 (MQ=255)

                                                          aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                             <  1:2948106/93‑1 (MQ=255)

                                                          aaaGATCAACGTATTGGCAGTAAATCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                             <  2:4164225/93‑1 (MQ=255)

                                                           aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                            <  2:1581769/93‑1 (MQ=255)

                                                            aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAGGAGTGAGAAACCTATGCt                           >  1:2676302/1‑93 (MQ=255)

                                                            aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCt                           <  1:2105326/93‑1 (MQ=255)

                                                            aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCt                           >  1:4169773/1‑93 (MQ=255)

                                                             gATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt                          <  2:3418504/93‑1 (MQ=255)

                                                              aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                         <  1:4348466/93‑1 (MQ=255)

                                                               tCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTc                        >  2:2262204/1‑93 (MQ=255)

                                                               tCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTc                        <  1:4080268/93‑1 (MQ=255)

                                                               tCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTc                        <  1:2293372/93‑1 (MQ=255)

                                                                ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                       <  2:3610428/93‑1 (MQ=255)

                                                                ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                       >  2:3488538/1‑93 (MQ=255)

                                                                 cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGTAGAGTGAGAAACCTATGCTTTCAc                      <  2:4554912/93‑1 (MQ=255)

                                                                  aCGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAcc                     <  2:1077749/93‑1 (MQ=255)

                                                                  aCGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAcc                     >  2:4505493/1‑93 (MQ=255)

                                                                   cGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCg                    >  2:3533919/1‑93 (MQ=255)

                                                                    gTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGc                   >  1:714433/1‑93 (MQ=255)

                                                                     tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                  >  1:944087/1‑93 (MQ=255)

                                                                     tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                  >  1:1078517/1‑93 (MQ=255)

                                                                       ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                >  1:575203/1‑93 (MQ=255)

                                                                       ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                >  1:2148999/1‑93 (MQ=255)

                                                                        tGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGc               >  2:3069504/1‑93 (MQ=255)

                                                                         ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg              >  2:859680/1‑93 (MQ=255)

                                                                         ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg              <  1:3271772/93‑1 (MQ=255)

                                                                          gCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCgg             >  2:4446717/1‑93 (MQ=255)

                                                                            aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg           >  1:434157/1‑93 (MQ=255)

                                                                            aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg           >  2:4189632/1‑93 (MQ=255)

                                                                             gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt          >  1:46824/1‑93 (MQ=255)

                                                                             gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt          <  1:2539842/93‑1 (MQ=255)

                                                                              tAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGtt         <  2:3153674/93‑1 (MQ=255)

                                                                               aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg        <  1:3847547/93‑1 (MQ=255)

                                                                                 aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa      >  2:401793/1‑93 (MQ=255)

                                                                                 aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa      >  2:2069984/1‑93 (MQ=255)

                                                                                   ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt    <  1:630834/93‑1 (MQ=255)

                                                                                   ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt    >  2:4661455/1‑93 (MQ=255)

                                                                                   ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt    <  1:1298511/93‑1 (MQ=255)

                                                                                    cccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAtt   <  1:2629328/93‑1 (MQ=255)

                                                                                    cccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAtt   >  2:62854/1‑93 (MQ=255)

                                                                                     ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa  >  2:2966156/1‑93 (MQ=255)

                                                                                     ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGGTGGAATTa  >  1:1009462/1‑93 (MQ=255)

GGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTA  >  NC_000913/4049052‑4049229

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