BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,850,117	C→T	T427T (ACC→ACT)	sppA →	protease IV, a signal peptide peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,850,117	0	C	T	100.0%	1131.0 / NA	259	T427T (ACC→ACT)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (120/139); total (120/139)

Only 100 of 262 total aligned reads displayed.

GCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCACCGC  >  NC_000913/1850019‑1850214

    gTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGAcc                                                                                             <  2:2809914/101‑1 (MQ=255)

     tAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCg                                                                                            <  1:2809867/101‑1 (MQ=255)

     tAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCg                                                                                            >  1:2810318/1‑101 (MQ=255)

     tAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCg                                                                                            >  2:2810319/1‑101 (MQ=255)

          cTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTATGGCTAACCCCAGCACTCTGACCGGTTCt                                                                                       <  1:2809846/101‑1 (MQ=255)

            gtgGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTAt                                                                                     >  2:2810324/1‑101 (MQ=255)

             tgGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATc                                                                                    >  1:2810326/1‑101 (MQ=255)

               gTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCgg                                                                                  <  2:2809886/101‑1 (MQ=255)

                ttGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGt                                                                                 <  1:2809975/101‑1 (MQ=255)

                 tGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTa                                                                                <  2:2810042/101‑1 (MQ=255)

                 tGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTa                                                                                >  1:2810330/1‑101 (MQ=255)

                  gTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTAt                                                                               <  1:2810081/101‑1 (MQ=255)

                   tATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATc                                                                              <  1:2809667/101‑1 (MQ=255)

                   tATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATc                                                                              >  2:2810331/1‑101 (MQ=255)

                   tATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATc                                                                              <  1:2809797/101‑1 (MQ=255)

                      cGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTc                                                                           >  1:2810335/1‑101 (MQ=255)

                       gATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCg                                                                          >  1:2810339/1‑101 (MQ=255)

                       gATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCg                                                                          >  2:2810337/1‑101 (MQ=255)

                       gATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCg                                                                          <  1:2810125/101‑1 (MQ=255)

                        aTGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCgg                                                                         <  2:2810013/101‑1 (MQ=255)

                         tGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGc                                                                        <  1:2809947/101‑1 (MQ=255)

                         tGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGc                                                                        >  1:2810341/1‑101 (MQ=255)

                         tGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCAGTTCTATCGGTATCTTCGGc                                                                        <  2:2809957/101‑1 (MQ=255)

                          gggCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCg                                                                       <  1:2809894/101‑1 (MQ=255)

                           ggcggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGt                                                                      <  1:2809891/101‑1 (MQ=255)

                            gcggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCCCTCTGACCGGTTCTATCGGTATCTTCGGCGTg                                                                     <  2:2809936/101‑1 (MQ=255)

                             cggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGa                                                                    <  2:2809872/101‑1 (MQ=255)

                             cggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGa                                                                    >  1:2810345/1‑101 (MQ=255)

                             cggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGa                                                                    <  1:2809929/101‑1 (MQ=255)

                              ggcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGAt                                                                   >  2:2810347/1‑101 (MQ=255)

                               gcatggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATc                                                                  <  2:2810080/101‑1 (MQ=255)

                                catggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATca                                                                 >  1:2810351/1‑101 (MQ=255)

                                 atggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATcac                                                                >  2:2810352/1‑101 (MQ=255)

                                   ggcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATcacca                                                              >  2:2810354/1‑101 (MQ=255)

                                    gcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATcaccac                                                             <  2:2810010/101‑1 (MQ=255)

                                    gcggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATcaccac                                                             <  1:2809939/101‑1 (MQ=255)

                                      ggcatCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCg                                                           <  1:2809734/101‑1 (MQ=255)

                                        catCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGCTCACCACCGTa                                                         <  1:2809806/101‑1 (MQ=255)

                                         atCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTGGGCGTGATCACCACCGTAg                                                        <  1:2810199/101‑1 (MQ=255)

                                         atCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAg                                                        >  1:2810359/1‑101 (MQ=255)

                                         atCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAg                                                        <  2:2810001/101‑1 (MQ=255)

                                           cTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGAt                                                                   >  1:2810363/1‑88 (MQ=255)

                                           cTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGaa                                                      <  2:2810112/101‑1 (MQ=255)

                                           cTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGaa                                                      >  2:2810361/1‑101 (MQ=255)

                                             ggtggtTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGaaaa                                                    >  1:2810366/1‑101 (MQ=255)

                                              gtggtTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAAt                                                   <  1:2809979/101‑1 (MQ=255)

                                              gtggtTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAAt                                                   <  1:2810021/101‑1 (MQ=255)

                                                 gtTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGt                                                >  2:2810369/1‑101 (MQ=255)

                                                  ttACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTc                                               >  1:2810372/1‑101 (MQ=255)

                                                  ttACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTc                                               <  1:2810086/101‑1 (MQ=255)

                                                   tACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCt                                              >  2:2810375/1‑101 (MQ=255)

                                                    aCTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTg                                             <  1:2810064/101‑1 (MQ=255)

                                                    aCTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTg                                             >  1:2810377/1‑101 (MQ=255)

                                                      tGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGa                                           <  1:2809963/101‑1 (MQ=255)

                                                      tGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGa                                           <  2:2810091/101‑1 (MQ=255)

                                                       ggATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGAt                                          <  2:2809849/101‑1 (MQ=255)

                                                       ggATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGAt                                          >  1:2810380/1‑101 (MQ=255)

                                                       ggATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGAt                                          <  2:2810103/101‑1 (MQ=255)

                                                        gATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGAtt                                         <  1:2809997/101‑1 (MQ=255)

                                                        gATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGAtt                                         >  2:2810381/1‑101 (MQ=255)

                                                         aTTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTc                                        >  2:2810385/1‑101 (MQ=255)

                                                         aTTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTc                                        >  1:2810386/1‑101 (MQ=255)

                                                         aTTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTc                                        <  1:2809777/101‑1 (MQ=255)

                                                          tttCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCg                                       <  2:2810120/101‑1 (MQ=255)

                                                           ttCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCGCCACCGTAGAAAATAGTCTGGATTCGa                                      <  2:2810038/101‑1 (MQ=255)

                                                            tCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGAt                                     <  2:2809990/101‑1 (MQ=255)

                                                             ccACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCg                                                                                            >  2:2810390/1‑45 (MQ=255)

                                                              cACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTg                                   >  1:2810392/1‑101 (MQ=255)

                                                              cACGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTg                                   <  1:2810071/101‑1 (MQ=255)

                                                               aCGCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTgg                                  >  2:2810395/1‑101 (MQ=255)

                                                                 gCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGgtg                                >  1:2810398/1‑101 (MQ=255)

                                                                 gCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGgtg                                >  2:2810399/1‑101 (MQ=255)

                                                                 gCCAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGgtg                                <  2:2810187/101‑1 (MQ=255)

                                                                  ccAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGgtgt                               >  2:2810402/1‑101 (MQ=255)

                                                                   cAGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGtt                              <  1:2809926/101‑1 (MQ=255)

                                                                    aGCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTc                             >  2:2810406/1‑101 (MQ=255)

                                                                     gCTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCa                            <  1:2810213/101‑1 (MQ=255)

                                                                      cTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTgg                                  >  2:2810408/1‑94 (MQ=255)

                                                                      cTAATTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCAt                           <  1:2809833/101‑1 (MQ=255)

                                                                        aaTTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATAc                         >  1:2810410/1‑101 (MQ=255)

                                                                         aTTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACt                        >  1:2810414/1‑101 (MQ=255)

                                                                         aTTACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACt                        <  2:2809713/101‑1 (MQ=255)

                                                                          ttACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATAc                         >  2:2810416/1‑99 (MQ=255)

                                                                          ttACATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTg                       <  2:2810165/101‑1 (MQ=255)

                                                                            aCATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGAt                     >  2:2810418/1‑101 (MQ=255)

                                                                             cATTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATg                    <  1:2810044/101‑1 (MQ=255)

                                                                              aTTGTGGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATgg                   >  1:2810420/1‑101 (MQ=255)

                                                                                tgtgGCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCAACGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATgg                   >  2:2810423/1‑99 (MQ=255)

                                                                                    gCTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGtctc             <  2:2810067/101‑1 (MQ=255)

                                                                                     cTAACCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCa            >  2:2810427/1‑101 (MQ=255)

                                                                                       aaCCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAAc          <  1:2810094/101‑1 (MQ=255)

                                                                                        aCCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACt         <  2:2809795/101‑1 (MQ=255)

                                                                                        aCCCCAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACt         >  2:2810429/1‑101 (MQ=255)

                                                                                         ccccAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACtt        <  2:2810285/101‑1 (MQ=255)

                                                                                          cccAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTc       <  1:2810012/101‑1 (MQ=255)

                                                                                           ccAGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCa      >  2:2810432/1‑101 (MQ=255)

                                                                                             aGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCAcc    <  2:2809982/101‑1 (MQ=255)

                                                                                             aGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCAcc    >  1:2810436/1‑101 (MQ=255)

                                                                                             aGCACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCAcc    <  1:2810148/101‑1 (MQ=255)

                                                                                               cACTCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCACCGc  >  1:2810438/1‑101 (MQ=255)

GCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGGTATCTTCGGCGTGATCACCACCGTAGAAAATAGTCTGGATTCGATTGGTGTTCATACTGATGGTGTCTCAACTTCACCGC  >  NC_000913/1850019‑1850214

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