BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,959,168	(C)_6→7	coding (99/747 nt)	cutC ←	protein CutC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,959,162	1	.	C	99.6%	943.0 / ‑3.2	232	coding (105/747 nt)	cutC	protein CutC
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (107/124); total (108/124)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 269 total aligned reads displayed.

CGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑CCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAG  >  NC_000913/1959062‑1959257

cGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAg                                                                                                  >  1:3188157/1‑101 (MQ=255)

cGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAg                                                                                                  >  2:3188156/1‑101 (MQ=255)

 gCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑c                                                                                                >  1:3188161/1‑101 (MQ=255)

 gCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑c                                                                                                >  2:3188163/1‑101 (MQ=255)

 gCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑c                                                                                                <  1:3187786/101‑1 (MQ=255)

 gCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑c                                                                                                <  2:3187711/101‑1 (MQ=255)

   tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccc                                                                                              >  2:3188170/1‑101 (MQ=255)

   tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccc                                                                                              >  1:3188166/1‑101 (MQ=255)

       gCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccc                                                                                          <  2:3187800/101‑2 (MQ=255)

        cAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccct                                                                                         >  2:3188174/1‑99 (MQ=255)

         aaaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccctc                                                                                        <  2:3187576/101‑4 (MQ=255)

          aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccctct                                                                                       >  1:3188177/1‑97 (MQ=255)

          aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccctct                                                                                       <  1:3187876/101‑5 (MQ=255)

          aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑ccccccctct                                                                                       <  2:3187932/101‑5 (MQ=255)

            aaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCttt                                                                                      <  2:3187821/101‑1 (MQ=255)

            aaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCttt                                                                                      <  2:3187481/101‑1 (MQ=255)

                accaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTggg                                                                                  >  2:3188179/1‑101 (MQ=255)

                 ccaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTgggg                                                                                 <  2:3187691/101‑1 (MQ=255)

                 ccaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTgggg                                                                                 >  1:3188180/1‑101 (MQ=255)

                 ccaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTgggg                                                                                 <  2:3187651/101‑1 (MQ=255)

                     cgcgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgc                                                                             <  1:3187657/101‑1 (MQ=255)

                      gcgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcg                                                                            <  1:3187806/101‑1 (MQ=255)

                      gcgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcg                                                                            >  1:3188186/1‑101 (MQ=255)

                       cgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcgc                                                                           <  2:3187817/101‑1 (MQ=255)

                       cgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcgc                                                                           <  1:3187957/101‑1 (MQ=255)

                       cgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcgc                                                                           >  1:3188190/1‑101 (MQ=255)

                       cgTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTgcgc                                                                           >  2:3188189/1‑101 (MQ=255)

                        gTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCa                                                                          <  1:3187915/101‑1 (MQ=255)

                         tGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCAt                                                                         <  1:3187573/101‑1 (MQ=255)

                              gAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTc                                                                    >  1:3188193/1‑101 (MQ=255)

                              gAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTc                                                                    >  2:3188198/1‑101 (MQ=255)

                              gAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTc                                                                    >  1:3188197/1‑101 (MQ=255)

                              gAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTc                                                                    <  1:3187692/101‑1 (MQ=255)

                               aaTTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCa                                                                   <  2:3187323/101‑1 (MQ=255)

                                 ttATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAAc                                                                 >  1:3188202/1‑101 (MQ=255)

                                  tATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAAct                                                                >  1:3188205/1‑101 (MQ=255)

                                     cGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctg                                                             <  1:3187666/101‑1 (MQ=255)

                                     cGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctg                                                             <  1:3187394/101‑1 (MQ=255)

                                     cGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctg                                                             >  1:3188209/1‑101 (MQ=255)

                                     cGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctg                                                             <  1:3187845/101‑1 (MQ=255)

                                        aTGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctgtct                                                          >  2:3188211/1‑101 (MQ=255)

                                        aTGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctgtct                                                          >  1:3188212/1‑101 (MQ=255)

                                        aTGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctgtct                                                          <  2:3187611/101‑1 (MQ=255)

                                         tGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctgtctg                                                         <  1:3187890/101‑1 (MQ=255)

                                           cacaGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTgcg                                                       >  2:3188214/1‑101 (MQ=255)

                                            acaGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTgcgc                                                      <  2:3187539/101‑1 (MQ=255)

                                             caGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGcc                                                     >  1:3188216/1‑101 (MQ=255)

                                              aGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCg                                                    <  2:3187623/101‑1 (MQ=255)

                                                ggATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGtt                                                  >  1:3188220/1‑101 (MQ=255)

                                                 gATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGttt                                                 <  1:3187884/101‑1 (MQ=255)

                                                  aTCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGtttt                                                <  1:3187796/101‑1 (MQ=255)

                                                  aTCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGtttt                                                <  1:3187429/101‑1 (MQ=255)

                                                   tCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTtg                                               >  2:3188225/1‑101 (MQ=255)

                                                   tCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTtg                                               <  1:3187726/101‑1 (MQ=255)

                                                      tCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTgc                                                        >  1:3188227/1‑89 (MQ=255)

                                                      tCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTtgctg                                            <  1:3187723/101‑1 (MQ=255)

                                                       cACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTtgctg                                            <  2:3187784/100‑1 (MQ=255)

                                                         cccGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTgcgc                                                      <  1:3187738/88‑1 (MQ=255)

                                                         cccGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTgcgc                                         >  1:3188231/1‑101 (MQ=255)

                                                         cccGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTgcg                                          <  2:3187834/100‑1 (MQ=255)

                                                          ccGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGcc                                        >  1:3188233/1‑101 (MQ=255)

                                                             cTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGtt                                     >  1:3188236/1‑101 (MQ=255)

                                                              tGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTa                                    >  2:3188237/1‑101 (MQ=255)

                                                              tGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTa                                    <  2:3187758/101‑1 (MQ=255)

                                                                gcgcACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGt                                  <  1:3187649/101‑1 (MQ=255)

                                                                 cgcACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTg                                 >  2:3188239/1‑101 (MQ=255)

                                                                  gcACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGt                                  <  1:3187909/99‑1 (MQ=255)

                                                                  gcACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGc                                <  1:3187767/101‑1 (MQ=255)

                                                                  gcACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGc                                >  2:3188244/1‑101 (MQ=255)

                                                                   cACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGca                               >  1:3188246/1‑101 (MQ=255)

                                                                     cGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGcaca                             >  1:3188249/1‑101 (MQ=255)

                                                                      ggATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACAt                            >  2:3188250/1‑101 (MQ=255)

                                                                        aTTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTc                          >  2:3188254/1‑101 (MQ=255)

                                                                        aTTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTc                          >  1:3188255/1‑101 (MQ=255)

                                                                        aTTTCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTc                          <  2:3187745/101‑1 (MQ=255)

                                                                         tttCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTcc                         <  1:3187945/101‑1 (MQ=255)

                                                                           tCAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCAt                       >  1:3188260/1‑101 (MQ=255)

                                                                            cAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATg                      <  2:3187763/101‑1 (MQ=255)

                                                                            cAGTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATg                      >  2:3188261/1‑101 (MQ=255)

                                                                              gTACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCt                    <  2:3187771/101‑1 (MQ=255)

                                                                               tACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTg                   <  2:3188003/101‑1 (MQ=255)

                                                                               tACACCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTg                   <  1:3187863/101‑1 (MQ=255)

                                                                                acacacAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGt                  >  2:3188265/1‑101 (MQ=255)

                                                                                 cacCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACAtt                           <  2:3187793/91‑1 (MQ=255)

                                                                                 cacCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTa                 >  1:3188266/1‑101 (MQ=255)

                                                                                  acCCAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTaa                <  2:3187843/101‑1 (MQ=255)

                                                                                   cccAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACtctgtctg                                                         <  1:3187838/59‑1 (MQ=255)

                                                                                   cccAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTCAc               <  1:3187812/101‑1 (MQ=255)

                                                                                   cccAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAAc               <  1:3188055/101‑1 (MQ=255)

                                                                                     cAACGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAg             <  2:3187882/101‑1 (MQ=255)

                                                                                      aaCGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGc            <  1:3187968/101‑1 (MQ=255)

                                                                                      aaCGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGc            >  2:3188270/1‑101 (MQ=255)

                                                                                      aaCGACGGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGc            >  2:3188273/1‑101 (MQ=255)

                                                                                        cgacgGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCaa          >  1:3188274/1‑101 (MQ=255)

                                                                                         gacgGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTa                 <  1:3187742/93‑1 (MQ=255)

                                                                                         gacgGCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCaaa         >  1:3188277/1‑101 (MQ=255)

                                                                                            ggCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAAttt      >  1:3188281/1‑101 (MQ=255)

                                                                                             gCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTc     <  1:3187992/101‑1 (MQ=255)

                                                                                             gCGTTAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTc     >  2:3188282/1‑101 (MQ=255)

                                                                                                ttAAGCCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAg  >  1:3188286/1‑101 (MQ=255)

CGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAG‑CCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAG  >  NC_000913/1959062‑1959257

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