BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,333,152	A→G	L16P (CTA→CCA)	yjdB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,333,152	0	A	G	97.6%	634.2 / 1.1	206	L16P (CTA→CCA)	yjdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/1); new base G (109/92); total (113/93)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  2:3510626/93‑1 (MQ=255)

 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  1:324175/93‑1 (MQ=255)

 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        >  2:3063622/1‑93 (MQ=255)

   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      >  1:775148/1‑93 (MQ=255)

   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      <  2:3828867/93‑1 (MQ=255)

   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      <  1:3040828/93‑1 (MQ=255)

    gTTTAAAAAAGGCAATATTCGGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgc                                                                                     <  2:4384572/93‑1 (MQ=255)

        aaaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCtt                                                                                 >  1:3071754/1‑93 (MQ=255)

         aaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTc                                                                                >  2:2792340/1‑93 (MQ=255)

         aaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTc                                                                                >  1:4092700/1‑93 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  2:1273035/93‑1 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  2:100785/93‑1 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              >  2:2314731/1‑93 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              >  1:4084264/1‑93 (MQ=255)

              ggCAATTTTTAGGCAGAGAGAGATATAAAAAGCGGCCAACAATAGTCAGGTGAGTAAATTCAAAGAGGGTCTTGTTGTTAGGCGCTGCAACAt                                                                           >  2:4220439/1‑93 (MQ=255)

              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  1:3739492/1‑93 (MQ=255)

              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  2:3705957/1‑93 (MQ=255)

              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  1:2301195/1‑93 (MQ=255)

                cAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCa                                                                         >  2:3657452/1‑93 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  2:2979332/1‑93 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  1:2437583/1‑93 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  1:4504764/1‑93 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        <  2:2175710/93‑1 (MQ=255)

                  atatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACg                                                                       >  2:1301150/1‑93 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      <  1:3204881/93‑1 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGGCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  1:4540649/1‑93 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  2:2901737/93‑1 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  1:4631966/93‑1 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  1:4410738/93‑1 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    <  1:3657993/93‑1 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:482019/1‑93 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:325242/1‑93 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:2654711/1‑93 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:2867810/1‑93 (MQ=255)

                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                <  2:3395792/93‑1 (MQ=255)

                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                <  2:3159377/93‑1 (MQ=255)

                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                <  1:3788395/93‑1 (MQ=255)

                            aGATAGAGATATAGAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCAt                                                             >  1:3134921/1‑93 (MQ=255)

                             gATAGAGATATAAAAAGAGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCAta                                                            >  1:1351905/1‑92 (MQ=255)

                              aTAGAGNTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCg                                                           >  1:99798/1‑93 (MQ=255)

                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         >  1:4111421/1‑93 (MQ=255)

                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  2:2186684/93‑1 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        <  2:546721/93‑1 (MQ=255)

                                   gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaa                                                      >  2:728669/1‑93 (MQ=255)

                                   gaTATAAAAAGCGGACAACAATAGCCAGGCGAGAAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATAACGGTGTTTCCATCGAACaa                                                      >  2:3842630/1‑93 (MQ=255)

                                    atatAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaaa                                                     <  2:4564434/93‑1 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  2:4330068/93‑1 (MQ=255)

                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   <  1:1493407/93‑1 (MQ=255)

                                        aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgc                                                 <  2:4178495/93‑1 (MQ=255)

                                          aaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcgc                                               <  2:2576260/93‑1 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  1:118900/1‑93 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  1:3959079/1‑93 (MQ=255)

                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  1:1046079/93‑1 (MQ=255)

                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        <  1:3252254/93‑1 (MQ=255)

                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  1:2701177/93‑1 (MQ=255)

                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  2:4474347/93‑1 (MQ=255)

                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgc                                     <  1:4004833/93‑1 (MQ=255)

                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgc                                     >  1:2539331/1‑93 (MQ=255)

                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGGGCGAATATGCTCgca                                    >  2:4113046/1‑93 (MQ=255)

                                                      aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaa                                   <  2:1395476/93‑1 (MQ=255)

                                                      aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaa                                   >  2:3605508/1‑93 (MQ=255)

                                                      aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTGGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaa                                   <  2:3118064/93‑1 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  2:3361531/93‑1 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  2:4749658/93‑1 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  1:4315991/93‑1 (MQ=255)

                                                          gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagca                               >  2:816906/1‑93 (MQ=255)

                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            >  2:3576323/1‑93 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCCCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  2:489141/1‑93 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTTCATCGAACAAAGTGCGCATATGCTCGCACAGAAAAAt                           >  1:1263373/1‑93 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  1:2629509/1‑93 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  2:872573/1‑93 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  2:1607563/1‑93 (MQ=255)

                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        <  1:2129946/93‑1 (MQ=255)

                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        >  2:3854933/1‑93 (MQ=255)

                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  >  1:541498/1‑93 (MQ=255)

                                                                         tCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCt                >  2:208422/1‑93 (MQ=255)

                                                                         tCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCt                >  1:2250114/1‑93 (MQ=255)

                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               <  2:121487/93‑1 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             >  1:558464/1‑93 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  2:2934966/93‑1 (MQ=255)

                                                                             agagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAg            <  2:239752/93‑1 (MQ=255)

                                                                             agagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAg            <  1:1816231/93‑1 (MQ=255)

                                                                             agagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAg            <  1:4674092/93‑1 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGGTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:3301804/1‑93 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  1:2280518/1‑93 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  1:906320/93‑1 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  2:1827505/93‑1 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  2:4158923/1‑93 (MQ=255)

                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         >  1:53551/1‑93 (MQ=255)

                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         >  1:3542164/1‑93 (MQ=255)

                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       <  2:3426824/93‑1 (MQ=255)

                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATACGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       >  2:581428/1‑93 (MQ=255)

                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  1:2746408/93‑1 (MQ=255)

                                                                                    cTTTTTGGGAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  2:1004401/93‑1 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGTAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  1:3002706/1‑93 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  2:1598592/1‑93 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    <  2:4725425/93‑1 (MQ=255)

                                                                                       tttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:2906021/93‑1 (MQ=255)

                                                                                       tttGGTAGGCGCTTCAACATCACGATGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:4642714/93‑1 (MQ=255)

                                                                                       tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:3863414/93‑1 (MQ=255)

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

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