BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,854,698	T→C	Y445C (TAC→TGC)	yidJ ←	putative sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,854,698	0	T	C	99.7%	1002.2 / ‑5.2	282	Y445C (TAC→TGC)	yidJ	putative sulfatase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (124/157); total (125/157)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 286 total aligned reads displayed.

ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGAT  >  NC_000913/3854608‑3854789

 cTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAAc                                                                                          <  2:3424986/93‑1 (MQ=255)

 cTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAAc                                                                                          >  2:2923600/1‑93 (MQ=255)

  tGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCCTTGGCAACt                                                                                         >  1:1029846/1‑93 (MQ=255)

   ggCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTg                                                                                        >  1:243797/1‑93 (MQ=255)

       aaTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGa                                                                                    >  2:1005649/1‑93 (MQ=255)

        aTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGaa                                                                                   >  1:380191/1‑93 (MQ=255)

         tAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAAc                                                                                  <  2:771103/93‑1 (MQ=255)

           gCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACgg                                                                                <  2:937948/93‑1 (MQ=255)

           gCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACgg                                                                                <  2:6093517/93‑1 (MQ=255)

            ccATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGa                                                                               >  2:3648183/1‑93 (MQ=255)

             cATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGAt                                                                              <  2:2462625/93‑1 (MQ=255)

             cATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGAt                                                                              <  2:6144613/93‑1 (MQ=255)

              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                             <  2:6207440/93‑1 (MQ=255)

              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                             >  1:1723933/1‑93 (MQ=255)

                cTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgc                                                                           <  1:742995/93‑1 (MQ=255)

                 ttGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgcg                                                                          <  1:5765436/93‑1 (MQ=255)

                  tgtgGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGa                                                                         >  2:1864113/1‑93 (MQ=255)

                     ggacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAtt                                                                      <  2:6323297/93‑1 (MQ=255)

                      gacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAttt                                                                     <  1:2852564/93‑1 (MQ=255)

                      gacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAttt                                                                     <  1:3289005/93‑1 (MQ=255)

                      gacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAttt                                                                     <  1:519838/93‑1 (MQ=255)

                           ggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCCAATTTTGTcc                                                                >  2:310656/1‑93 (MQ=255)

                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                               >  2:3939797/1‑93 (MQ=255)

                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                               <  2:2016376/93‑1 (MQ=255)

                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                               >  1:2492927/1‑93 (MQ=255)

                              cgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATg                                                             <  1:3927901/93‑1 (MQ=255)

                              cgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATg                                                             <  1:1197316/93‑1 (MQ=255)

                              cgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATg                                                             <  2:1154606/93‑1 (MQ=255)

                                aaaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTTCATGTa                                                           >  1:5643361/1‑93 (MQ=255)

                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                          <  2:598946/93‑1 (MQ=255)

                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                          <  2:5360154/93‑1 (MQ=255)

                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                          >  2:6396026/1‑93 (MQ=255)

                                   cGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATc                                                        <  2:3538144/93‑1 (MQ=255)

                                   cGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATc                                                        <  1:4431673/93‑1 (MQ=255)

                                      ccccATCCAGCTCGGCCGTGCATCTTTACGCCGCGGACTCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                     <  2:529577/93‑1 (MQ=255)

                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                     >  2:663638/1‑93 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                    <  2:2562580/93‑1 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                    <  1:2337045/93‑1 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                    >  1:4976664/1‑93 (MQ=255)

                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                  >  2:5003733/1‑93 (MQ=255)

                                          aTCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAg                                                 <  2:275818/93‑1 (MQ=255)

                                          aTCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAg                                                 <  1:4120909/93‑1 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                                <  2:5893062/93‑1 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                                <  1:1652675/93‑1 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                                >  1:4176794/1‑93 (MQ=255)

                                            ccAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGc                                               >  2:3849871/1‑93 (MQ=255)

                                            ccAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGc                                               >  2:4450020/1‑93 (MQ=255)

                                            ccAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGc                                               <  1:1350592/93‑1 (MQ=255)

                                            ccAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGc                                               <  2:2355755/93‑1 (MQ=255)

                                              aGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGt                                             <  1:3530751/93‑1 (MQ=255)

                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                            <  2:4091373/93‑1 (MQ=255)

                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                            >  2:4167804/1‑93 (MQ=255)

                                                cgcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCa                                           >  1:358667/1‑93 (MQ=255)

                                                cgcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCa                                           <  1:4236577/93‑1 (MQ=255)

                                                 gcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCAt                                          <  1:5479164/93‑1 (MQ=255)

                                                  cgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATg                                         <  1:917196/93‑1 (MQ=255)

                                                   ggCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGc                                        >  1:4407241/1‑93 (MQ=255)

                                                    gCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCa                                       >  1:10910/1‑93 (MQ=255)

                                                    gCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCa                                       <  1:1899465/93‑1 (MQ=255)

                                                    gCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCa                                       <  2:4832824/93‑1 (MQ=255)

                                                      cGTGCATCTTTACGCCACGGACTCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtt                                     <  2:5126823/93‑1 (MQ=255)

                                                      cGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtt                                     <  1:6337371/93‑1 (MQ=255)

                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                    <  2:2159248/93‑1 (MQ=255)

                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                    <  2:5731274/93‑1 (MQ=255)

                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                   <  2:324143/93‑1 (MQ=255)

                                                         gCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtg                                  <  1:2699437/93‑1 (MQ=255)

                                                          cATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgc                                 <  2:4687538/93‑1 (MQ=255)

                                                           aTCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgct                                <  1:5839477/93‑1 (MQ=255)

                                                             cTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctgc                              <  2:457907/93‑1 (MQ=255)

                                                              tttACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCg                             >  1:4889054/1‑93 (MQ=255)

                                                              tttACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCg                             >  1:2606731/1‑93 (MQ=255)

                                                                tACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGaa                           >  1:1461084/1‑93 (MQ=255)

                                                                 aCGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAAc                          <  2:2697456/93‑1 (MQ=255)

                                                                   gCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGt                        <  1:1518/93‑1 (MQ=255)

                                                                   gCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGt                        >  2:5813652/1‑93 (MQ=255)

                                                                   gCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGt                        <  2:5549825/93‑1 (MQ=255)

                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                       <  1:4483247/93‑1 (MQ=255)

                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                       <  1:4644784/93‑1 (MQ=255)

                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                       >  1:1247842/1‑93 (MQ=255)

                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                       <  1:6448537/93‑1 (MQ=255)

                                                                     cACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCt                      <  1:6061231/93‑1 (MQ=255)

                                                                      aCGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTg                     <  1:2908574/93‑1 (MQ=255)

                                                                       cGGACACAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGc                    <  1:5701563/93‑1 (MQ=255)

                                                                        ggACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCa                   >  1:4558618/1‑93 (MQ=255)

                                                                          aCGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCaaa                 >  1:2957246/1‑93 (MQ=255)

                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGa            <  2:3799662/93‑1 (MQ=255)

                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGa            <  1:6232957/93‑1 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat           >  1:4845519/1‑93 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat           <  1:6157552/93‑1 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat           >  2:5672024/1‑93 (MQ=255)

                                                                                  tCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGatat         >  2:4632658/1‑93 (MQ=255)

                                                                                  tCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGatat         >  2:4511333/1‑93 (MQ=255)

                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCTTGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACCGATatc        >  2:4287348/1‑93 (MQ=255)

                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc        <  2:1530426/93‑1 (MQ=255)

                                                                                     aTTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcat      <  2:3987122/93‑1 (MQ=255)

                                                                                      ttGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc     >  1:2076119/1‑93 (MQ=255)

                                                                                      ttGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc     >  1:5299748/1‑93 (MQ=255)

                                                                                       tGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAACCGGATATCATCg    <  1:4814395/93‑1 (MQ=255)

                                                                                       tGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCg    <  1:3082869/93‑1 (MQ=255)

                                                                                         gTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGAt  >  2:5974105/1‑93 (MQ=255)

ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGAT  >  NC_000913/3854608‑3854789

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