BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,365	A→G	L173P (CTT→CCT)	fecB ←	FecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,365	0	A	G	98.7%	494.9 / ‑3.1	156	L173P (CTT→CCT)	fecB	FecB
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (74/80); total (75/81)
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.89e-01

Only 100 of 165 total aligned reads displayed.

AGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATA  >  NC_000913/4511273‑4511454

  gCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgc                                                                                          <  1:265248/93‑1 (MQ=255)

  gCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgc                                                                                          <  2:2314626/93‑1 (MQ=255)

   ccAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgct                                                                                         <  2:2239813/93‑1 (MQ=255)

   ccAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgct                                                                                         <  1:4544424/93‑1 (MQ=255)

      gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggc                                                                                      <  1:141653/93‑1 (MQ=255)

      gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggc                                                                                      >  2:720592/1‑93 (MQ=255)

       cACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAggctggct                                                                                     <  2:2752538/93‑1 (MQ=255)

       cACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAggctggct                                                                                     >  2:4551993/1‑93 (MQ=255)

        aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                    >  2:3207570/1‑93 (MQ=255)

        aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                    <  2:4266900/93‑1 (MQ=255)

        aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                    >  2:869397/1‑93 (MQ=255)

        aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTGCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                    <  2:356385/93‑1 (MQ=255)

         cGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGctggctgg                                                                                   <  1:1529927/93‑1 (MQ=255)

             gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGGCCa                                                                               >  2:926632/1‑93 (MQ=255)

              cgggTGCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCAc                                                                              <  1:1738089/91‑1 (MQ=255)

                 gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                           >  1:4135112/1‑93 (MQ=255)

                 gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                           >  2:447649/1‑93 (MQ=255)

                 gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                           >  1:1615956/1‑93 (MQ=255)

                 gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                           >  2:3965169/1‑93 (MQ=255)

                   ccAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgcgc                                                                         <  1:2202496/93‑1 (MQ=255)

                    cAGGTCTCCTGAGTATGCAGATTGAATTGCTTTGCGCGTGATGTTCCAGAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGcc                                                                        <  2:3310561/93‑1 (MQ=255)

                     aGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCa                                                                       >  1:1427714/1‑93 (MQ=255)

                      ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCAt                                                                      >  2:1410957/1‑93 (MQ=255)

                       gTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATc                                                                     <  2:2494203/93‑1 (MQ=255)

                         atcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCct                                                                   <  1:272397/92‑1 (MQ=255)

                          tcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctc                                                                  <  2:728803/93‑1 (MQ=255)

                           ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGATGGCCCACTGCGCCTTCctct                                                                 >  2:1398151/1‑93 (MQ=255)

                            cTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctc                                                                >  2:1598021/1‑93 (MQ=255)

                            cTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctc                                                                >  1:1648923/1‑93 (MQ=255)

                             tGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctct                                                               >  1:2708244/1‑93 (MQ=255)

                              gAGTATGCGGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCtt                                                              <  2:1536961/93‑1 (MQ=255)

                              gAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCtt                                                              <  1:3169817/93‑1 (MQ=255)

                               aGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCttt                                                             <  1:1989828/93‑1 (MQ=255)

                               aGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCttt                                                             >  2:1089430/1‑93 (MQ=255)

                                gTATGCAGGTTGAATTGCTTTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTa                                                            <  2:1953154/93‑1 (MQ=255)

                                gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTa                                                            >  1:4658390/1‑93 (MQ=255)

                                    gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtg‑tt                                                       >  2:1925944/1‑93 (MQ=255)

                                      aGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtg‑ttgt                                                     >  2:3623832/1‑93 (MQ=255)

                                        gTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTc                                                   <  1:3826403/93‑1 (MQ=255)

                                        gTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTc                                                   >  1:3311158/1‑93 (MQ=255)

                                          tGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCa                                                 <  2:2378770/93‑1 (MQ=255)

                                           gAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAg                                                <  2:908261/93‑1 (MQ=255)

                                            aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGa                                               <  1:4079622/93‑1 (MQ=255)

                                            aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGa                                               >  1:3365210/1‑93 (MQ=255)

                                            aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGa                                               <  2:2074134/93‑1 (MQ=255)

                                             aTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGAc                                              >  1:689600/1‑93 (MQ=255)

                                             aTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGAc                                              >  1:3149645/1‑93 (MQ=255)

                                             aTTGCTGTTCGCGTGATGTGCCAAAGCCATCACGTGTCCCTTTGGTAGGCTGTGTGGCCCAATGCGCCGTCCTCTCTTTATG‑TTGTTCTAGAc                                              <  2:3142089/93‑1 (MQ=255)

                                             aTTGCTGTTCGCGTGAAGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTCATG‑TTGTTCCAGAc                                              >  2:3131185/1‑93 (MQ=255)

                                               tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGt                                            <  2:1452565/93‑1 (MQ=255)

                                               tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGt                                            >  1:1313849/1‑93 (MQ=255)

                                               tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGt                                            <  1:3172429/93‑1 (MQ=255)

                                                  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTTCAGACGTggc                                         >  2:3092437/1‑91 (MQ=255)

                                                  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGcc                                         >  1:1911994/1‑93 (MQ=255)

                                                  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGcc                                         <  1:2547331/93‑1 (MQ=255)

                                                  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGcc                                         >  2:2188010/1‑93 (MQ=255)

                                                  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGcc                                         >  1:4569255/1‑93 (MQ=255)

                                                  tGTTCGCGTGATGTGCAAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGcc                                         <  1:4749351/93‑1 (MQ=255)

                                                   gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCt                                        <  2:674458/93‑1 (MQ=255)

                                                   gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCt                                        >  1:2469684/1‑93 (MQ=255)

                                                    ttCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTg                                       <  1:4619136/93‑1 (MQ=255)

                                                     tCGCGGGATGTGCCAATGGCCTCGCGTGGCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGc                                      <  2:488975/93‑1 (MQ=255)

                                                      cgcgGGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCCCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCa                                     <  1:1677787/93‑1 (MQ=255)

                                                       gcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAt                                    <  2:3242781/93‑1 (MQ=255)

                                                        cgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtc                                   >  1:282807/1‑93 (MQ=255)

                                                         gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtct                                  >  1:1240017/1‑93 (MQ=255)

                                                          tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctc                                 >  2:3878773/1‑93 (MQ=255)

                                                          tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctc                                 <  2:3675807/93‑1 (MQ=255)

                                                          tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctc                                 <  2:2227460/93‑1 (MQ=255)

                                                            aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCTCTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctctc                               <  1:58304/93‑1 (MQ=255)

                                                            aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCCTCCTCTCTTTATGTTTGTTCCA‑ACGCGCCTGCAtctctg                               >  2:2600422/1‑92 (MQ=255)

                                                            aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctctc                               <  1:4667917/93‑1 (MQ=255)

                                                            aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCAtctctc                               <  1:3754290/93‑1 (MQ=255)

                                                             tgtgCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCg                              <  2:3375792/93‑1 (MQ=255)

                                                                  cAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCtttt                         <  2:3983839/93‑1 (MQ=255)

                                                                   aaaGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGGCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCCTATCTCGGtttta                        >  2:4546388/1‑92 (MQ=255)

                                                                    aaGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTa                       >  1:281038/1‑93 (MQ=255)

                                                                     aGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAc                      >  1:3270634/1‑93 (MQ=255)

                                                                     aGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGACTGCATCTCTCGCTTTTTAc                      >  1:2046073/1‑93 (MQ=255)

                                                                       gCCACGCGTGTCCCTTTGGGAGGATGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAccc                    <  2:2628584/93‑1 (MQ=255)

                                                                        ccACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAcccc                   >  2:1946858/1‑92 (MQ=255)

                                                                         cACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACccac                  <  2:633609/93‑1 (MQ=255)

                                                                          aCGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                 >  1:1695239/1‑93 (MQ=255)

                                                                            gcgTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAt               <  1:4000537/93‑1 (MQ=255)

                                                                            gcgTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAt               <  2:3594754/93‑1 (MQ=255)

                                                                             cgTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACATGCCTGCATCTCTCGCTTTTTACCCACCAtt              >  1:2189589/1‑93 (MQ=255)

                                                                              gtgtCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAttt             <  1:716196/93‑1 (MQ=255)

                                                                              gtgtCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAttt             >  1:2386773/1‑93 (MQ=255)

                                                                              gtgtCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAttt             <  2:1731267/93‑1 (MQ=255)

                                                                              gtgtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGGCTGCATCTCTCGCTTTTTACCCACCAttt             >  2:789994/1‑93 (MQ=255)

                                                                               tgtCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTc            >  2:1495089/1‑93 (MQ=255)

                                                                                 tCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGc          >  2:4198079/1‑93 (MQ=255)

                                                                                 tCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGc          >  1:2995642/1‑93 (MQ=255)

                                                                                  cccTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc         >  1:175139/1‑93 (MQ=255)

                                                                                   ccTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCg        >  2:856363/1‑93 (MQ=255)

                                                                                     tttGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgat      <  1:3506755/93‑1 (MQ=255)

                                                                                       tGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatga    <  1:862768/93‑1 (MQ=255)

                                                                                        gggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat   <  2:4265225/93‑1 (MQ=255)

                                                                                        gggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat   <  2:1724750/93‑1 (MQ=255)

                                                                                         ggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATa  <  2:2997662/93‑1 (MQ=255)

AGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATG‑TTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATA  >  NC_000913/4511273‑4511454

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