BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,686,056	(T)_6→5	coding (346/1647 nt)	fumA ←	fumarase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,686,051	0	T	.	95.9%	419.9 / 9.8	98	coding (351/1647 nt)	fumA	fumarase A
Reads supporting (aligned to +/- strand): ref base T (3/1); new base . (53/41); total (56/42)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 104 total aligned reads displayed.

GCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGC  >  NC_000913/1685960‑1686141

gCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCtt                                                                                           <  2:214771/93‑1 (MQ=255)

  gAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCtttt                                                                                         <  1:3215687/93‑1 (MQ=255)

   aGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCttttt                                                                                        <  2:4662427/93‑1 (MQ=255)

       gcgcAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCtttttacca                                                                                    <  1:2898506/93‑5 (MQ=255)

         gcAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACcaac                                                                                 <  1:4506572/93‑1 (MQ=255)

           aGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACcaacaa                                                                               >  2:641365/1‑93 (MQ=255)

           aGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACcaacaa                                                                               <  2:2662609/93‑1 (MQ=255)

           aGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACcaacaa                                                                               <  1:871445/93‑1 (MQ=255)

            taTTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaat                                                                              >  1:59272/2‑93 (MQ=255)

            gATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaat                                                                              >  1:4465747/1‑93 (MQ=255)

            gATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaat                                                                              <  1:537999/93‑1 (MQ=255)

              ttATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaataa                                                                            >  1:398303/1‑93 (MQ=255)

              ttATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaataa                                                                            >  1:1311167/1‑93 (MQ=255)

              ttATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaataa                                                                            >  1:2283072/1‑93 (MQ=255)

               tATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaataat                                                                           <  1:394257/93‑1 (MQ=255)

               tATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACaataat                                                                           <  1:3029960/93‑1 (MQ=255)

                aTCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATc                                                                          >  1:2080819/1‑93 (MQ=255)

                aTCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATc                                                                          >  1:4183276/1‑93 (MQ=255)

                 tCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                                                          >  1:3005719/1‑93 (MQ=255)

                   ttCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATcgcg                                                                       >  2:3404328/1‑93 (MQ=255)

                   ttCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATcgcg                                                                       >  1:3398655/1‑93 (MQ=255)

                    tCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGGCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCgg                                                                      >  1:1021395/1‑93 (MQ=255)

                     cGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGt                                                                     <  2:3639794/93‑1 (MQ=255)

                     cGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGt                                                                     <  2:4416166/93‑1 (MQ=255)

                      gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTg                                                                    >  1:4342631/1‑93 (MQ=255)

                       atatAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGc                                                                   >  1:612214/1‑93 (MQ=255)

                       atatAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGc                                                                   >  2:1192736/1‑93 (MQ=255)

                       atatAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGc                                                                   >  2:4536298/1‑93 (MQ=255)

                        tataAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGcc                                                                  >  2:3536522/1‑93 (MQ=255)

                         ataAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCg                                                                 >  1:1724909/1‑93 (MQ=255)

                          taAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCgg                                                                <  1:1156749/93‑1 (MQ=255)

                          taAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCgg                                                                >  1:1419594/1‑93 (MQ=255)

                             gtgtTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTAt                                                             >  2:2887603/1‑93 (MQ=255)

                             gtgtTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTAt                                                             >  2:2164963/1‑93 (MQ=255)

                              tgtTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATc                                                            >  1:20888/1‑93 (MQ=255)

                               gtTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATcc                                                           >  1:4251273/1‑93 (MQ=255)

                                ttATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                          >  2:317553/1‑93 (MQ=255)

                                ttATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                          <  1:1411893/93‑1 (MQ=255)

                                ttATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                          >  2:1699938/1‑93 (MQ=255)

                                 tataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTg                                                         >  1:425258/1‑93 (MQ=255)

                                  ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                        >  2:2434694/1‑93 (MQ=255)

                                  ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                        <  2:2378575/93‑1 (MQ=255)

                                  ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                        <  2:4493400/93‑1 (MQ=255)

                                  ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                        >  2:532197/1‑93 (MQ=255)

                                  ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                        >  1:2050153/1‑93 (MQ=255)

                                   taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGAc                                                       >  1:4369919/1‑93 (MQ=255)

                                   taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGAc                                                       >  2:717606/1‑93 (MQ=255)

                                   taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TCTTTACCAACAATAATCGCGGTGCCGGTATCCTGAc                                                       >  1:1660556/1‑93 (MQ=255)

                                      acacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAgg                                                    <  1:3257268/93‑1 (MQ=255)

                                        acCGCGCGCCTGCGCCGCTTCCTCACCACCACCGGGCAATACACGCTGCCC‑TTATTACTAACACTAATCGCTGTGCCGGCATCCCGACAggct                                                  >  2:918290/1‑91 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATCCACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  2:2524281/1‑93 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  2:1903250/1‑93 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  <  1:1559937/93‑1 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  1:194249/1‑93 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  1:3948471/1‑93 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  1:4288275/1‑93 (MQ=255)

                                        acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                                  >  1:2531188/1‑93 (MQ=255)

                                           gcgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGc                                               >  1:794270/1‑93 (MQ=255)

                                           gcgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGc                                               >  1:64001/1‑93 (MQ=255)

                                                ccAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAAc                                          >  1:2308717/1‑93 (MQ=255)

                                                  aGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                        <  1:1761939/93‑1 (MQ=255)

                                                   gcgcCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                                       <  2:111168/93‑1 (MQ=255)

                                                   gcgcCGCTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                                       <  1:339110/93‑1 (MQ=255)

                                                    ccccgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                      <  1:442421/91‑1 (MQ=255)

                                                     gccgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                     >  2:1836160/1‑93 (MQ=255)

                                                     gccgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                     >  2:844171/1‑93 (MQ=255)

                                                      ccgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCttt                                    >  1:1450498/1‑93 (MQ=255)

                                                       cgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                   <  1:4070808/93‑1 (MQ=255)

                                                       cgcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                   >  2:1806360/1‑93 (MQ=255)

                                                        gcTTCATCACCACCACCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcg                                  >  1:3234943/1‑93 (MQ=255)

                                                               caccaccacCGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGAt                           <  2:1155098/93‑1 (MQ=255)

                                                                 cccccaccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGTACGCCTTTCGCCGCGATGt                         <  2:4561395/90‑1 (MQ=25)

                                                                  caccaccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTc                        <  1:4199809/93‑1 (MQ=255)

                                                                  caccaccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTc                        <  2:4358461/93‑1 (MQ=255)

                                                                    ccaccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCgg                      <  2:3101095/93‑1 (MQ=255)

                                                                     caccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGa                     <  2:4188840/93‑1 (MQ=255)

                                                                       ccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                   <  2:208504/93‑1 (MQ=255)

                                                                       ccGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                   <  2:577979/93‑1 (MQ=255)

                                                                        cGGTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGtt                  >  2:3696578/1‑93 (MQ=255)

                                                                         ggTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTa                 <  1:269265/93‑1 (MQ=255)

                                                                         ggTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTa                 <  2:3984125/93‑1 (MQ=255)

                                                                          gTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTAc                <  1:409404/93‑1 (MQ=255)

                                                                          gTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTAc                <  1:1853466/93‑1 (MQ=255)

                                                                          gTCCATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTAc                <  2:3871693/93‑1 (MQ=255)

                                                                            ccATACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGc              >  1:3014952/1‑93 (MQ=255)

                                                                              aTACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg            <  1:3964014/93‑1 (MQ=255)

                                                                              aTACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg            <  2:411912/93‑1 (MQ=255)

                                                                              aTACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg            <  1:4655404/93‑1 (MQ=255)

                                                                               tACACGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg           >  1:2223452/1‑93 (MQ=255)

                                                                                acacGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGa          >  1:616616/1‑93 (MQ=255)

                                                                                 cacGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa         >  2:1793970/1‑93 (MQ=255)

                                                                                 cacGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa         <  2:165759/93‑1 (MQ=255)

                                                                                   cGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAtt       <  2:2197611/93‑1 (MQ=255)

                                                                                   cGCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAtt       <  1:4384216/93‑1 (MQ=255)

                                                                                    gCTGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg      <  2:2325440/93‑1 (MQ=255)

                                                                                     cTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg      >  2:2955795/1‑93 (MQ=255)

                                                                                      tGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACATGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa    >  2:1988052/1‑93 (MQ=255)

                                                                                      tGCCC‑TTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa    <  2:1040796/93‑1 (MQ=255)

                                                                                         ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAgc  >  1:2058489/1‑93 (MQ=255)

                                                                                         ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAgc  <  2:71694/93‑1 (MQ=255)

GCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGC  >  NC_000913/1685960‑1686141

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