BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,150,120	T→C	66.9%	F153F (TTT→TTC)	yghA →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,120	0	T	C	66.9%	183.9 / 128.7	157	F153F (TTT→TTC)	yghA	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (19/33); new base C (41/63); total (60/97)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 160 total aligned reads displayed.

AAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCA  >  NC_000913/3150023‑3150219

aaGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCa                                                                                                   <  1:886060/100‑1 (MQ=255)

  gcgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaaa                                                                                                 <  1:2181337/100‑1 (MQ=255)

   cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaaaa                                                                                                <  1:1898103/100‑1 (MQ=255)

   cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                                <  1:1064741/100‑1 (MQ=255)

   cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                                <  1:2656027/100‑1 (MQ=255)

   cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                                <  2:1219005/100‑1 (MQ=255)

   cgcACAAGGCGTTAGGCGGGCTGGATACTATGGCGCTGGTCGCTGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                                <  2:2325550/100‑1 (MQ=255)

        aaGGCGTTAGGCGGGCTGGATATTATGGCGGTGGGCGCCGGGAAACAGGATGCCATACCGGATATTGCAGAACTCACCAGCGAACAGTTTCAAAAGAc                                                                                             >  1:1551723/1‑98 (MQ=255)

           gCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTg                                                                                        <  2:397129/100‑1 (MQ=255)

             gTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGcc                                                                                      >  2:2123881/1‑100 (MQ=255)

              ttAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCa                                                                                     <  1:2522837/100‑1 (MQ=255)

               tAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCAt                                                                                    <  2:302101/100‑1 (MQ=255)

                 ggCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTa                                                                                  <  1:561684/100‑1 (MQ=255)

                 ggCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTa                                                                                  >  2:3299417/1‑100 (MQ=255)

                  gCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTaa                                                                                 >  1:110033/1‑100 (MQ=255)

                    gggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACg                                                                               <  2:479799/100‑1 (MQ=255)

                      gCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGtt                                                                             >  1:2664886/1‑100 (MQ=255)

                         ggATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACgagatc                                                                          >  2:627580/1‑95 (MQ=255)

                         ggATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTc                                                                          <  2:1635866/100‑1 (MQ=255)

                           atatTATGGCGCTGGTCGCCGGGATACAGGTCGCCATTCCGGATATTGCCGACCTCACCCGCGAACAGTTCCAAAATGCCTTTGCCATTAACGTTTTcg                                                                         <  1:619644/99‑1 (MQ=255)

                           atatTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcgc                                                                        >  1:2469538/1‑100 (MQ=255)

                            tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCGATTAACGTTTTcgcg                                                                       >  1:31066/1‑100 (MQ=255)

                            tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTcgcg                                                                       >  2:2959981/1‑100 (MQ=255)

                             attatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCgcgc                                                                      <  1:3299389/100‑1 (MQ=255)

                             attatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCgcgc                                                                      <  1:1235143/100‑1 (MQ=255)

                               tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTg                                                                    >  1:1733337/1‑100 (MQ=255)

                                atGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGt                                                                   <  1:935033/100‑1 (MQ=255)

                                  ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTc                                                                 >  2:3227043/1‑100 (MQ=255)

                                  ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGtt                                                                  >  2:3065922/1‑99 (MQ=255)

                                  ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTc                                                                 >  2:185245/1‑100 (MQ=255)

                                  ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTGCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTc                                                                 >  1:2290914/1‑100 (MQ=255)

                                    cgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTg                                                               >  2:1207871/1‑100 (MQ=255)

                                     gcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTgg                                                              >  2:3105541/1‑100 (MQ=255)

                                      cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                             <  1:180684/100‑1 (MQ=255)

                                      cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                             <  1:3377170/100‑1 (MQ=255)

                                      cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                             <  2:1668987/100‑1 (MQ=255)

                                       tGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCt                                                            <  2:561537/100‑1 (MQ=255)

                                           agCCCGGAAACAGGTTGCCATTACGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGACATTAACGTTTTCGCGCTGTTCTAGCTAAcc                                                        >  1:789985/2‑100 (MQ=255)

                                            gCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAccc                                                       <  1:2148153/100‑1 (MQ=255)

                                            gCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAccc                                                       >  1:2769331/1‑100 (MQ=255)

                                             ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                      <  2:2735735/100‑1 (MQ=255)

                                             ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                      <  1:2201888/100‑1 (MQ=255)

                                              cGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAg                                                     <  2:1541700/100‑1 (MQ=255)

                                              cGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTCGCCATTAACGTTTTCGCTCTGTTCTGGCTAAcc                                                        <  1:457080/97‑1 (MQ=255)

                                               gggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                    >  1:1301280/1‑100 (MQ=255)

                                               gggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                    >  2:3221023/1‑100 (MQ=255)

                                                   aaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGc                                                <  1:3067573/100‑1 (MQ=255)

                                                    aCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCg                                               >  2:828167/1‑100 (MQ=255)

                                                     cAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGa                                              <  1:176446/100‑1 (MQ=255)

                                                      aGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGAt                                             <  2:1092113/100‑1 (MQ=255)

                                                      aGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAACAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGAt                                             <  1:3045218/100‑1 (MQ=255)

                                                      aGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGAt                                             <  2:953640/100‑1 (MQ=255)

                                                       ggTTGCCATTCCTGATATTGCAGTCCTCACCCGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTACCCCAGGAAGCGATc                                            <  1:1343226/100‑1 (MQ=255)

                                                       ggTTGCCATTCCGGATATTGCAGCCCTCACCAGCGAACAGTTCCTAAAGACATTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATc                                            <  1:1622672/100‑1 (MQ=255)

                                                       ggTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATc                                            <  1:96633/100‑1 (MQ=255)

                                                       ggTTGCCACTCCGGATATTGCAGACCTCACCAGCGACCAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATc                                            <  2:2373785/100‑1 (MQ=255)

                                                             cATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTg                                      >  1:2884793/1‑100 (MQ=255)

                                                              aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                     <  2:795561/100‑1 (MQ=255)

                                                               ttCCGGATATTGCAGACCTCACCAGCGAACAGTTACAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                     <  2:1532032/99‑1 (MQ=255)

                                                                tCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTa                                   >  2:1796743/1‑100 (MQ=255)

                                                                 ccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAc                                  <  1:3142529/100‑1 (MQ=255)

                                                                 ccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAc                                  <  2:1945406/100‑1 (MQ=255)

                                                                  cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 <  2:1063534/100‑1 (MQ=255)

                                                                  cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 <  2:1047994/100‑1 (MQ=255)

                                                                  cGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 >  1:111450/1‑100 (MQ=255)

                                                                   ggATATTGCAGACCTCACCAGCTAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCg                                <  1:3337925/100‑1 (MQ=255)

                                                                    gATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGa                               <  2:2750299/100‑1 (MQ=255)

                                                                         tGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                          <  2:2350516/100‑1 (MQ=255)

                                                                          gCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                         >  2:753435/1‑100 (MQ=255)

                                                                          gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAAGTg                         >  1:147729/1‑100 (MQ=255)

                                                                           cAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGc                        >  1:3218639/1‑100 (MQ=255)

                                                                           cAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGc                        >  1:3293057/1‑100 (MQ=255)

                                                                            aGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAGAGGTGCa                       <  2:1778608/100‑1 (MQ=255)

                                                                            aGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCa                       <  2:1077721/100‑1 (MQ=255)

                                                                             gACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTATGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCaa                      >  2:749813/1‑100 (MQ=255)

                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:1371855/100‑1 (MQ=255)

                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:1372050/100‑1 (MQ=255)

                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:2941384/100‑1 (MQ=255)

                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:110033/100‑1 (MQ=255)

                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  1:2594752/100‑1 (MQ=255)

                                                                               ccTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGt                    >  1:1799436/1‑100 (MQ=255)

                                                                               ccTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGt                    >  2:1497058/1‑100 (MQ=255)

                                                                                 tCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTat                  >  1:656191/1‑100 (MQ=255)

                                                                                  cACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatc                 <  2:1911396/100‑1 (MQ=255)

                                                                                   aCCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatca                >  2:1390436/1‑100 (MQ=255)

                                                                                   aCCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatca                <  2:265424/100‑1 (MQ=255)

                                                                                    ccAGCGAAAAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcat               >  1:1001753/1‑100 (MQ=255)

                                                                                     cAGCGAACAGTTTCAACAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcatc              <  1:1984922/100‑1 (MQ=255)

                                                                                     cAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcatc              <  1:447280/100‑1 (MQ=255)

                                                                                     cAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcatc              <  1:2762388/100‑1 (MQ=255)

                                                                                      aGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTAtcatca             <  1:2274499/100‑1 (MQ=255)

                                                                                      aGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTAtcatca             <  1:762343/100‑1 (MQ=255)

                                                                                        cGAACAGTTCCAAAAGACCTTTGCCATTATCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc           <  1:530095/100‑1 (MQ=255)

                                                                                        cGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAgatcggaagag           >  2:2268676/1‑90 (MQ=255)

                                                                                        cGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc           <  2:2568440/100‑1 (MQ=255)

                                                                                         gAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacca          <  2:245885/100‑1 (MQ=255)

                                                                                            cAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACtt       >  1:1982793/1‑100 (MQ=255)

                                                                                             aGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTc      <  2:303157/100‑1 (MQ=255)

                                                                                              gTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg     >  2:684611/1‑100 (MQ=255)

                                                                                                ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTGCCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  2:3218667/100‑1 (MQ=255)

AAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCA  >  NC_000913/3150023‑3150219

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