BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,074,121	(T)_5→6	coding (82/966 nt)	yggF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,074,116	1	.	T	98.8%	942.5 / 1.6	251	coding (87/966 nt)	yggF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base T (135/112); total (136/115)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

TCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑T‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGACATCATT  >  NC_000913/3074024‑3074203

tCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta                                                                                           <  2:6602573/93‑1 (MQ=255)

tCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta                                                                                           <  2:4618459/93‑1 (MQ=255)

 cAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta‑t                                                                                         <  2:6685251/93‑1 (MQ=255)

   aTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta‑t‑tt                                                                                      >  1:3311476/1‑93 (MQ=255)

   aTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta‑t‑tt                                                                                      <  1:3447118/93‑1 (MQ=255)

    tttCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttta‑t‑ttt                                                                                     >  1:1038190/1‑93 (MQ=255)

     ttCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑tttt                                                                                    >  2:1764587/1‑93 (MQ=255)

     ttCTCCTTCGCCGATAACCACTCGCCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑tttt                                                                                    <  2:4751976/93‑1 (MQ=255)

      tctcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑ttttt                                                                                   <  2:6107688/93‑2 (MQ=255)

      tctcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑ttttt                                                                                   <  2:2894214/93‑2 (MQ=255)

         ccTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑tttttgtc                                                                                >  1:4643298/1‑89 (MQ=255)

         ccTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑t‑tttttgtc                                                                                <  2:2156253/93‑5 (MQ=255)

           ttCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGc                                                                               <  1:1750706/93‑1 (MQ=255)

              gCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCcac                                                                            >  2:3691189/1‑93 (MQ=255)

               ccGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCcaca                                                                           >  1:6533339/1‑93 (MQ=255)

                cGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACAt                                                                          <  2:95812/93‑1 (MQ=255)

                 gATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATc                                                                         <  2:5347696/93‑1 (MQ=255)

                    aaCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCgg                                                                      <  1:2619259/93‑1 (MQ=255)

                     aCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGt                                                                     >  2:1030481/1‑93 (MQ=255)

                     aCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGt                                                                     >  2:4141658/1‑93 (MQ=255)

                     aCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGt                                                                     >  1:2531096/1‑93 (MQ=255)

                      ccACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGtt                                                                    >  1:1134766/1‑93 (MQ=255)

                       cACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGttt                                                                   >  2:1333012/1‑93 (MQ=255)

                       cACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGttt                                                                   >  1:4803350/1‑93 (MQ=255)

                       cACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGttt                                                                   <  1:443432/93‑1 (MQ=255)

                        aCTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTg                                                                  >  1:166629/1‑93 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGc                                                                 >  2:490315/1‑93 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTAGT‑TTTTGTCGCCACATCCGGTTTGc                                                                 <  1:3237632/93‑1 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGGGACCGCCAGGCCATCACTTTTATT‑TTTTGTCGCCACATCCGGGTTGc                                                                 >  1:3729762/1‑93 (MQ=255)

                          ccGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCg                                                                <  1:1602621/92‑1 (MQ=255)

                             tCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGcc                                                             <  1:828500/93‑1 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAg                                                           >  2:4037160/1‑93 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAg                                                           <  2:2575346/93‑1 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAg                                                           >  1:5122573/1‑93 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAg                                                           >  1:4990508/1‑93 (MQ=255)

                                aCGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAgg                                                          >  1:706683/1‑93 (MQ=255)

                                   aaaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAg                                                       >  1:5232444/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGc                                                      >  1:4911387/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGc                                                      <  1:3612222/93‑1 (MQ=255)

                                      aGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCa                                                    <  2:4878838/93‑1 (MQ=255)

                                      aGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCa                                                    >  2:1584524/1‑93 (MQ=255)

                                       gCGACGTCGTTTAATGCCTGGCGCATCGCGGTGATCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCaa                                                   >  2:4354225/1‑93 (MQ=255)

                                        cGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAAc                                                  <  2:5813280/93‑1 (MQ=255)

                                         gACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACg                                                 >  2:3860724/1‑93 (MQ=255)

                                          aCGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCCACGc                                                >  1:1318149/1‑93 (MQ=255)

                                          aCGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGc                                                >  2:3601264/1‑93 (MQ=255)

                                          aCGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGc                                                >  2:5174220/1‑93 (MQ=255)

                                          aCGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGc                                                >  2:5211069/1‑93 (MQ=255)

                                          aCGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGc                                                <  1:5664152/93‑1 (MQ=255)

                                           cgtcgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCt                                               >  2:5700704/1‑93 (MQ=255)

                                              cgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGcc                                            <  1:6633786/93‑1 (MQ=255)

                                              cgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGcc                                            >  2:4803952/1‑93 (MQ=255)

                                               gtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCt                                           >  2:2444101/1‑93 (MQ=255)

                                               gtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCt                                           <  1:3653735/93‑1 (MQ=255)

                                                tttAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTg                                          >  1:285631/1‑93 (MQ=255)

                                                 ttAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGt                                         >  2:3621472/1‑93 (MQ=255)

                                                  tAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGtt                                        >  2:1875418/1‑93 (MQ=255)

                                                  tAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGtt                                        <  1:2722772/93‑1 (MQ=255)

                                                   aaTGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTc                                       <  1:6230330/93‑1 (MQ=255)

                                                   aaTGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTc                                       >  2:6331025/1‑93 (MQ=255)

                                                     tGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCg                                     >  2:5883693/1‑93 (MQ=255)

                                                      gCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGt                                    <  2:2738760/93‑1 (MQ=255)

                                                         tGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAAc                                 <  2:5144622/93‑1 (MQ=255)

                                                             gcATCGCGGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGa                             >  2:3422288/1‑93 (MQ=255)

                                                                 cgcgGTGCCCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaata                         <  1:2331595/93‑1 (MQ=255)

                                                                  gcgGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaataa                        >  1:6082335/1‑93 (MQ=255)

                                                                   cgGTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaataat                       <  1:2460887/93‑1 (MQ=255)

                                                                    ggTGACCGCGAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaataatt                      <  2:868865/93‑2 (MQ=255)

                                                                    ggTGACCGCCAGGCGATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATg                      <  1:2363944/93‑1 (MQ=255)

                                                                    ggTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATg                      <  1:6415063/93‑1 (MQ=255)

                                                                    ggTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATg                      <  1:3971704/93‑1 (MQ=255)

                                                                     gTGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATgg                     <  1:3063552/93‑1 (MQ=255)

                                                                      tGACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGc                    >  2:5553730/1‑93 (MQ=255)

                                                                       gACCGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGcc                   <  1:6381404/93‑1 (MQ=255)

                                                                       gACCGCCAGGCCATCAATTTTA‑TGTTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGcc                   <  1:2275947/93‑1 (MQ=255)

                                                                         ccGCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCCGTTCCGTAACGCGAAATAATGGCCAc                 >  1:2639061/1‑93 (MQ=255)

                                                                           gCCAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGCAATAATGGCCACGc               >  2:1435418/1‑93 (MQ=255)

                                                                            ccAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGcc              <  1:6348273/93‑1 (MQ=255)

                                                                             cAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCa             >  2:3496582/1‑93 (MQ=255)

                                                                             cAGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCa             <  2:2293982/93‑1 (MQ=255)

                                                                              aGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            >  1:3805733/1‑93 (MQ=255)

                                                                              aGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            <  1:4352131/93‑1 (MQ=255)

                                                                              aGGCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            >  1:1476003/1‑93 (MQ=255)

                                                                               ggCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgg           >  2:155936/1‑93 (MQ=255)

                                                                               ggCCATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgg           >  1:2183861/1‑93 (MQ=255)

                                                                               ggCCATCAATTTTA‑T‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAggg          <  1:1176552/93‑1 (MQ=255)

                                                                                  cATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgggcc        >  2:2987808/1‑91 (MQ=255)

                                                                                  cATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAc        >  2:74246/1‑93 (MQ=255)

                                                                                  cATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAc        <  1:4874411/93‑1 (MQ=255)

                                                                                  cATCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGACACGCCAGGGAc        >  2:6581786/1‑93 (MQ=255)

                                                                                    tCAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcat      <  1:5929427/93‑1 (MQ=255)

                                                                                     cAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatc     <  2:3948825/93‑1 (MQ=255)

                                                                                     cAATTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatc     <  2:3151669/93‑1 (MQ=255)

                                                                                      aaTTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    <  2:6004459/93‑1 (MQ=255)

                                                                                      aaTTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    >  1:5843287/1‑93 (MQ=255)

                                                                                      aaTTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    >  2:5648429/1‑93 (MQ=255)

                                                                                      aaTTTTATT‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    <  2:1664614/93‑1 (MQ=255)

                                                                                       attttatt‑ttTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatcat   >  1:5512253/1‑93 (MQ=255)

                                                                                        ttttatt‑ttTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTTCGTAACGCGGAATAATGGCCACCGCAAGGAcatccat  >  2:1166708/1‑90 (MQ=255)

                                                                                        ttttatt‑ttTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGACATCAtt  >  1:4326311/1‑93 (MQ=255)

TCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTA‑T‑TTTTGTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGACATCATT  >  NC_000913/3074024‑3074203

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