BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,149	G→A	G75S (GGC→AGC)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,149	0	G	A	98.5%	937.4 / ‑0.6	266	G75S (GGC→AGC)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (3/1); new base A (130/132); total (133/133)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 269 total aligned reads displayed.

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

  cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                                       <  2:190113/93‑1 (MQ=255)

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

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      >  2:312123/1‑93 (MQ=255)

    cAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGa                                                                                     >  1:624011/1‑93 (MQ=255)

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

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

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

      gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACt                                                                                   >  2:3894321/1‑93 (MQ=255)

        ggtcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 <  1:2329539/90‑1 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 >  2:3925015/1‑93 (MQ=255)

         ggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGCATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaa                                                                                >  1:977740/1‑93 (MQ=255)

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

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

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

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

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            >  2:1954058/1‑93 (MQ=255)

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

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

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

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

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

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

                        tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 >  2:1684143/1‑93 (MQ=255)

                        tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 <  1:631282/93‑1 (MQ=255)

                         cAGGCAGCAAAGGCCAGAACGCATCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGc                                                                <  2:2580365/93‑1 (MQ=255)

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

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

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  1:3950068/93‑1 (MQ=255)

                                      ccAGAACGCACCAAAAGATCCACGTCTTGGCAGTAAAACACCTGTTCCATTGAGCGTGACTGCAAAAGTCACCAAACAGCACCCACAGAagag                                                   >  2:1156256/1‑93 (MQ=255)

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

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

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

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

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGACAAAGTCACCAAACAGCCCCAACCGAagaggg                                                 >  1:1355952/1‑91 (MQ=255)

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                                 >  1:1279293/1‑93 (MQ=255)

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                                 >  1:2298939/1‑93 (MQ=255)

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

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

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

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

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

                                              cACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAc                                           <  1:952741/93‑1 (MQ=255)

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

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

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

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

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

                                               aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                          >  1:1591582/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  1:4048549/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  1:4046791/93‑1 (MQ=255)

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

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

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

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGGCACCAAACAGCACAAAACGAAGAGTGAGAACCCTACGc                                     >  1:1110746/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                             gTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGc                            <  2:4542039/93‑1 (MQ=255)

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

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

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           <  1:4635548/93‑1 (MQ=255)

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

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

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

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

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

                                                                   gCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCgg                      <  2:4796/93‑1 (MQ=255)

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

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

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

                                                                       tAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGtt                  >  1:4131178/1‑93 (MQ=255)

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

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

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              >  2:1541715/1‑93 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  1:1255256/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                              ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           <  1:2675875/93‑1 (MQ=255)

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

                                                                                  ttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTgg       <  1:2836050/93‑1 (MQ=255)

                                                                                    ccATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaa     <  2:3979084/93‑1 (MQ=255)

                                                                                    ccATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaa     <  2:2947786/93‑1 (MQ=255)

                                                                                     cATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaaa    <  2:4503479/93‑1 (MQ=255)

                                                                                     cATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaaa    <  2:1861794/93‑1 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGCTGGAATTACTGGAAACg  >  2:3757922/1‑93 (MQ=255)

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

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