BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,132,875	T→C	L360P (CTG→CCG)	flgE →	hook protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,132,875	0	T	C	98.4%	395.6 / ‑2.8	127	L360P (CTG→CCG)	flgE	hook protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (56/69); total (57/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 133 total aligned reads displayed.

ACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAAC  >  NC_000913/1132785‑1132965

aCAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCgg                                                                                          <  1:3576849/93‑1 (MQ=255)

  aaCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGaa                                                                                        <  2:931207/93‑1 (MQ=255)

  aaCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGaa                                                                                        <  2:1936599/93‑1 (MQ=255)

  aaCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGaa                                                                                        >  1:2335952/1‑93 (MQ=255)

   aCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAg                                                                                       >  2:1952609/1‑93 (MQ=255)

   aCGTCTGGTCTGCGAAGCAATCTTCTGGCTTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACACTGACCAACGGAGCGCCGGaac                                                                                       >  2:1952607/1‑92 (MQ=255)

      tctggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGTCCAACGGTGCGCCGGAAGGGt                                                                                    >  1:3502387/1‑93 (MQ=255)

      tctggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGt                                                                                    >  2:4577317/1‑93 (MQ=255)

      tctggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGt                                                                                    >  2:2627787/1‑93 (MQ=255)

          gtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCaa                                                                                <  1:3721025/93‑1 (MQ=255)

          gtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCaa                                                                                >  2:3704322/1‑93 (MQ=255)

           tctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAAc                                                                               >  1:3842388/1‑93 (MQ=255)

           tctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTCGCACCCTGACCAACGCTGCGCCGGAAGCGTCCAAc                                                                               >  1:3842390/1‑93 (MQ=255)

             tgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGt                                                                             <  2:1211873/93‑1 (MQ=255)

              gCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTc                                                                            >  2:3553776/1‑93 (MQ=255)

                gACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGa                                                                          <  1:1123373/93‑1 (MQ=255)

                gACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACAAACGGTGCGCCGGAAGCGTCCAACGTCgc                                                                          >  1:275994/1‑92 (MQ=255)

                 aCGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAt                                                                         <  1:4660422/93‑1 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCCCCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        <  2:2787424/93‑1 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        <  2:2787423/93‑1 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        >  2:214331/1‑93 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        >  2:2030500/1‑93 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        <  1:1420755/93‑1 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        <  2:1917501/93‑1 (MQ=255)

                  cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                        <  2:1083361/93‑1 (MQ=255)

                    cAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtctc                                                                      <  1:2309295/93‑1 (MQ=255)

                    cAATCTGCTGGCGTGGCGCTGTTAAGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtctc                                                                      <  1:4211364/93‑1 (MQ=255)

                     aaTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCa                                                                     <  2:1362414/93‑1 (MQ=255)

                     aaTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGCGCGCCGGAAGCGTCCAACGTCGATCTCa                                                                     >  2:1549833/1‑93 (MQ=255)

                      aTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAg                                                                    >  1:267771/1‑93 (MQ=255)

                       tcttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGt                                                                   >  2:4347030/1‑93 (MQ=255)

                       tcttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGt                                                                   <  2:1104646/93‑1 (MQ=255)

                       tcttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGt                                                                   >  2:3914325/1‑93 (MQ=255)

                         ttctGGCGTTGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCCCCCTGACCAACGGTGCTCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                 <  1:2908235/93‑1 (MQ=255)

                         ttctGGCGTGGCGCTGTTGGTGACAGCCGGGACGGGAAACTTTGGGACCCTGACCAACGGAGCGCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                 >  2:4572292/1‑93 (MQ=255)

                         ttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                 >  1:1066576/1‑93 (MQ=255)

                         ttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                 >  1:1434096/1‑93 (MQ=255)

                           ctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAg                                                               <  2:2771438/93‑1 (MQ=255)

                            tggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGa                                                              >  1:570544/1‑93 (MQ=255)

                             ggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGaa                                                             <  1:4248841/93‑1 (MQ=255)

                              gcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAAc                                                            <  1:3010354/93‑1 (MQ=255)

                                gtggcgCTGTGGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTg                                                          <  2:1140498/93‑1 (MQ=255)

                                 tggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTgg                                                         >  2:3945330/1‑93 (MQ=255)

                                 tggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTgg                                                         <  1:3784813/93‑1 (MQ=255)

                                  ggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGt                                                        <  1:350463/93‑1 (MQ=255)

                                      cTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAat                                                    <  1:1832047/93‑1 (MQ=255)

                                       tGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAata                                                   <  2:1739985/93‑1 (MQ=255)

                                       tGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAata                                                   <  2:2428335/93‑1 (MQ=255)

                                       tGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAata                                                   <  2:2471277/93‑1 (MQ=255)

                                        gTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAatat                                                  >  2:3334419/1‑93 (MQ=255)

                                          tGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGa                                                <  1:4622275/93‑1 (MQ=255)

                                          tGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGa                                                <  1:4622276/93‑1 (MQ=255)

                                          tGGGGACAGCCAGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGa                                                <  1:2114122/93‑1 (MQ=255)

                                           ggggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGAt                                               >  1:2773871/1‑93 (MQ=255)

                                           ggggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGAt                                               <  1:3873962/93‑1 (MQ=255)

                                           ggggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGAt                                               >  1:3298332/1‑93 (MQ=255)

                                            gggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATc                                              <  1:3682197/93‑1 (MQ=255)

                                              gACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGt                                            <  2:288647/93‑1 (MQ=255)

                                              gACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGt                                            <  2:4163810/93‑1 (MQ=255)

                                              gACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGt                                            >  1:1212946/1‑93 (MQ=255)

                                               aCAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGtt                                           >  1:3920821/1‑93 (MQ=255)

                                               aCAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGGCAATATGATCGtc                                           >  1:2939387/1‑92 (MQ=255)

                                                 aGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGc                                         <  1:864609/93‑1 (MQ=255)

                                                  gCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGcc                                        >  1:104461/1‑93 (MQ=255)

                                                      ggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGc                                    >  1:3056267/1‑93 (MQ=255)

                                                      ggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGc                                    >  2:49898/1‑93 (MQ=255)

                                                      ggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGc                                    <  1:1989899/93‑1 (MQ=255)

                                                         cgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTa                                 >  1:777958/1‑93 (MQ=255)

                                                          gggaCACTTTGGCACCGTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                                <  2:547112/93‑1 (MQ=255)

                                                          gggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                                >  2:103390/1‑93 (MQ=255)

                                                          gggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                                >  2:4458701/1‑93 (MQ=255)

                                                          gggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                                <  2:232426/93‑1 (MQ=255)

                                                           ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                               <  2:830180/93‑1 (MQ=255)

                                                           ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                               >  1:4609535/1‑93 (MQ=255)

                                                              aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                            >  2:2726271/1‑93 (MQ=255)

                                                              aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                            <  2:2766664/93‑1 (MQ=255)

                                                              aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                            <  2:4636734/93‑1 (MQ=255)

                                                              aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                            <  1:3188937/93‑1 (MQ=255)

                                                               aCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGCCAATATGATCGTTGCCCAGCGTAACTATc                           <  2:391210/93‑1 (MQ=255)

                                                                cTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCa                          <  1:2928798/93‑1 (MQ=255)

                                                                 tttGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAg                         >  2:2660806/1‑93 (MQ=255)

                                                                   tGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAGCGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTc                       <  2:1274410/93‑1 (MQ=255)

                                                                   tGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTc                       <  2:3260654/93‑1 (MQ=255)

                                                                      cACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTaa                    >  2:1667/1‑93 (MQ=255)

                                                                       aCCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAAc                   <  1:3907787/93‑1 (MQ=255)

                                                                          cTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGcc                >  1:4403855/1‑93 (MQ=255)

                                                                           tGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGccc               <  2:2814700/93‑1 (MQ=255)

                                                                           tGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGccc               <  2:496022/93‑1 (MQ=255)

                                                                           tGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGccc               >  1:330256/1‑93 (MQ=255)

                                                                           tGACCAACGGGGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAAAGccc               >  1:330258/1‑93 (MQ=255)

                                                                            gACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCa              <  2:554680/93‑1 (MQ=255)

                                                                             aCCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAg             <  1:1483342/93‑1 (MQ=255)

                                                                              ccAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGa            >  1:1808469/1‑93 (MQ=255)

                                                                               cAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAc           <  2:1288493/93‑1 (MQ=255)

                                                                                aaCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAcc          >  1:2050267/1‑93 (MQ=255)

                                                                                 aCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCa         <  2:3587676/93‑1 (MQ=255)

                                                                                 aCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCa         <  1:2381728/93‑1 (MQ=255)

                                                                                 aCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCa         <  2:2760972/93‑1 (MQ=255)

                                                                                    gTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCa      <  1:655483/93‑1 (MQ=255)

                                                                                     tGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCaa     >  1:1800949/1‑93 (MQ=255)

ACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAAC  >  NC_000913/1132785‑1132965

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