BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,962,070	A→T	P335P (CCT→CCA)	flhA ←	FlhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,962,070	0	A	T	99.1%	374.6 / NA	110	P335P (CCT→CCA)	flhA	FlhA
Reads supporting (aligned to +/- strand): ref base A (1/0); major base T (58/50); minor base G (0/1); total (59/51)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAG  >  NC_000913/1961980‑1962158

cATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTgg                                                                                        <  1:656267/93‑1 (MQ=255)

 aTTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGTTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGgc                                                                                       <  1:3218179/93‑1 (MQ=255)

 aTTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGgc                                                                                       >  2:1753250/1‑93 (MQ=255)

 aTTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGgc                                                                                       >  1:3038518/1‑93 (MQ=255)

    cccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCt                                                                                    <  1:461684/93‑1 (MQ=255)

    cccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCt                                                                                    >  2:1644786/1‑93 (MQ=255)

    cccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCt                                                                                    <  2:3973633/93‑1 (MQ=255)

     ccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtt                                                                                   <  1:1380792/93‑1 (MQ=255)

     ccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtt                                                                                   >  1:3318228/1‑93 (MQ=255)

     ccAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtt                                                                                   >  2:4647168/1‑93 (MQ=255)

       agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtttt                                                                                 >  1:1145692/1‑93 (MQ=255)

       agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtttt                                                                                 <  2:3536120/93‑1 (MQ=255)

       agagAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCtttt                                                                                 >  1:2009243/1‑93 (MQ=255)

        gataaTGTTCGAGTTGTCCATCCTCCCACGTCGCTGCGACAACGGTATTATTCCCTCCCCTTTTTACCGGTTTGGGTTCGGCTTGCGCttttt                                                                                <  2:3552548/90‑1 (MQ=255)

           aaTCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGtt                                                                             >  2:3671003/1‑93 (MQ=255)

           aaTCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGtt                                                                             >  1:2676255/1‑93 (MQ=255)

               ttCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTcgcg                                                                         >  1:1276732/1‑93 (MQ=255)

               ttCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTcgcg                                                                         >  1:1618115/1‑93 (MQ=255)

                  cAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTcc                                                                      >  1:2139402/1‑93 (MQ=255)

                  cAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTcc                                                                      >  1:3352062/1‑93 (MQ=255)

                  cAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTcc                                                                      <  1:2218898/93‑1 (MQ=255)

                  cAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTcc                                                                      >  1:3199333/1‑93 (MQ=255)

                   aGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCcg                                                                     <  1:3378360/93‑1 (MQ=255)

                    gTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCcgc                                                                    >  1:2326214/1‑93 (MQ=255)

                     ttGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTGTACCGGTTTGGGTTCGGCGGGCGCTTTTTGTTCGCGTCcgcg                                                                   <  1:2019479/93‑1 (MQ=255)

                         aCATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATc                                                               <  1:1467330/93‑1 (MQ=255)

                            tCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATccac                                                            >  2:1220077/1‑93 (MQ=255)

                            tCGTTCCACGTCACTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATccac                                                            <  1:232984/93‑1 (MQ=255)

                             cGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATccacc                                                           <  1:2892888/93‑1 (MQ=255)

                              gTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATccacca                                                          >  2:2740231/1‑93 (MQ=255)

                               ttCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAg                                                         <  2:3933665/93‑1 (MQ=255)

                                 ccACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAggt                                                       >  1:3014082/1‑92 (MQ=255)

                                  cACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAggtc                                                      <  2:162142/93‑3 (MQ=255)

                                    cGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAg                                                    >  1:448144/1‑93 (MQ=255)

                                     gTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGc                                                   >  1:2497923/1‑93 (MQ=255)

                                     gTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGc                                                   >  2:2298819/1‑93 (MQ=255)

                                      tCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGcc                                                  >  1:2694544/1‑93 (MQ=255)

                                        gCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCg                                                >  1:1637074/1‑93 (MQ=255)

                                        gCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCg                                                <  1:4476622/93‑1 (MQ=255)

                                         cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTGAGCCCGa                                               >  2:1023847/1‑93 (MQ=255)

                                         cTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGa                                               <  2:3094384/93‑1 (MQ=255)

                                          ttCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAg                                              >  1:440192/1‑93 (MQ=255)

                                            cGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTTAGCCCGAGCa                                            >  2:3794642/1‑93 (MQ=255)

                                              aCAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAAt                                          >  1:35433/1‑93 (MQ=255)

                                              aCAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCATATCCACCAGGTCAGCCCGAGCaac                                          >  1:957488/1‑92 (MQ=255)

                                                aaCGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATcc                                        >  1:661895/1‑93 (MQ=255)

                                                  cGGTATTATTCTCTGCCGAATTTAACGGTTTGGATTCGGCTGGCGATTTTTGTTCGCGTCAGCGTATCCACCAAGTCAGCCCGAGCAATccag                                      <  1:4221452/93‑3 (MQ=255)

                                                  aggTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCgg                                      <  1:4276286/92‑1 (MQ=255)

                                                   ggTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGc                                     <  2:1212544/93‑1 (MQ=255)

                                                   ggTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGc                                     <  1:1234366/93‑1 (MQ=255)

                                                   cgTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGc                                     <  1:3140638/92‑1 (MQ=255)

                                                    gTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCa                                    >  2:4282449/1‑93 (MQ=255)

                                                     tattatTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAg                                   <  1:747819/93‑1 (MQ=255)

                                                     tattatTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAg                                   <  1:4613730/93‑1 (MQ=255)

                                                     tattatTCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAg                                   <  1:4613731/93‑1 (MQ=255)

                                                      attattCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGGCATCCGGCAGt                                  >  2:3501370/1‑93 (MQ=255)

                                                       ttattCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTg                                 >  2:4726627/1‑93 (MQ=255)

                                                        tattCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGGTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCACGCAGTGa                                >  1:4046547/1‑93 (MQ=255)

                                                        tattCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGa                                >  1:4046549/1‑93 (MQ=255)

                                                        tattCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGa                                >  1:4046548/1‑93 (MQ=255)

                                                          ttCTCTGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAAc                              >  1:1354504/1‑93 (MQ=255)

                                                            ctctGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAg                            <  2:3203712/93‑1 (MQ=255)

                                                            ctctGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAg                            <  2:3599473/93‑1 (MQ=255)

                                                            ctctGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAg                            >  2:1272598/1‑93 (MQ=255)

                                                             tctGCCATTTTTTCCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGc                           <  2:3551049/93‑1 (MQ=255)

                                                             tctGCCATTTTTACCAGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGc                           >  2:1364403/1‑93 (MQ=255)

                                                               tGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaa                         >  2:1692543/1‑93 (MQ=255)

                                                               tGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaa                         >  2:3821523/1‑93 (MQ=255)

                                                               tGCCATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaa                         >  2:703618/1‑93 (MQ=255)

                                                                 ccATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaaaa                       >  1:102800/1‑93 (MQ=255)

                                                                 ccATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaaaa                       >  1:4394876/1‑93 (MQ=255)

                                                                 ccATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaaaa                       >  1:4064495/1‑93 (MQ=255)

                                                                 ccATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCaaaa                       >  1:4042066/1‑93 (MQ=255)

                                                                  cATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaaaaa                      >  1:71837/1‑93 (MQ=255)

                                                                  cATTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCaaaaa                      <  1:1857700/93‑1 (MQ=255)

                                                                   aTTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAAt                     >  1:2780846/1‑93 (MQ=255)

                                                                   aTTTTTACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAAt                     >  1:3838375/1‑93 (MQ=255)

                                                                    tttttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATa                    <  1:4064282/93‑1 (MQ=255)

                                                                    tttttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATa                    >  2:2267989/1‑93 (MQ=255)

                                                                    tttttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATa                    <  1:2957892/93‑1 (MQ=255)

                                                                      tttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCCGGTCAGCCCGAGCAATCCGGCAGTGAACCGCAAAAATAcc                  >  1:2095180/1‑93 (MQ=255)

                                                                      tttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAcc                  <  2:4003172/93‑1 (MQ=255)

                                                                       ttACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCa                 <  1:449626/93‑1 (MQ=255)

                                                                        tACCGGTTTGGGTTCGGGTGGCGCTTTTTGTTCGCGTCCGCGTATTCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATAACAg                >  1:2175918/1‑93 (MQ=255)

                                                                        tACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAg                <  2:3099753/93‑1 (MQ=255)

                                                                        tACCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAg                <  2:665914/93‑1 (MQ=255)

                                                                         aCCGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAgg               <  1:870475/93‑1 (MQ=255)

                                                                          ccGGTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGt              >  1:2264251/1‑93 (MQ=255)

                                                                            ggTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTc            <  2:3868502/93‑1 (MQ=255)

                                                                            ggTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTc            <  1:3579622/93‑1 (MQ=255)

                                                                             gTTTGGGTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCg           <  2:2282960/93‑1 (MQ=255)

                                                                              tttGGGTTCGGCTGGCGCTTTTTTTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCACGTTCgg          <  2:4606691/93‑1 (MQ=255)

                                                                                 gggTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt       <  2:4483329/93‑1 (MQ=255)

                                                                                 gggTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt       <  1:1113528/93‑1 (MQ=255)

                                                                                 gggTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAt       <  1:3933007/93‑1 (MQ=255)

                                                                                  ggTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCAtt      >  1:1965433/1‑93 (MQ=255)

                                                                                   gTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTc     >  1:1868574/1‑93 (MQ=255)

                                                                                   gTTCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTc     <  1:939960/93‑1 (MQ=255)

                                                                                    ttCGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTcc    <  2:4130995/93‑1 (MQ=255)

                                                                                      cGGCTGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGTCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAg  <  1:1644786/93‑1 (MQ=255)

CATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCGGTTTGGGTTCGGCAGGCGCTTTTTGTTCGCGTCCGCGTATCCACCAGGCCAGCCCGAGCAATCCGGCAGTGAACAGCAAAAATACCAGGTTCGGCATTCCAG  >  NC_000913/1961980‑1962158

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