BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	464,356	A→T	15.5%	intergenic (+48/‑46)	ybaV → / → fadM	putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	464,356	0	A	T	15.5%	430.6 / 67.4	168	intergenic (+48/‑46)	ybaV/fadM	putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand): ref base A (65/77); new base T (10/16); total (75/93)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 169 total aligned reads displayed.

GGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTC  >  NC_000913/464258‑464455

gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAg                                                                                                    <  1:611434/100‑1 (MQ=255)

gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAg                                                                                                    <  2:3139873/100‑1 (MQ=255)

   aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                 >  2:3488437/1‑100 (MQ=255)

    tGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                                >  2:2637306/1‑100 (MQ=255)

      ggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                              >  1:1310817/1‑100 (MQ=255)

      ggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                              >  2:2828279/1‑100 (MQ=255)

        cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                            <  2:796077/100‑1 (MQ=255)

        cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                            <  2:2408484/100‑1 (MQ=255)

          aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACctc                                                                                          <  1:476834/100‑1 (MQ=255)

           ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  1:3133818/100‑1 (MQ=255)

           ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  2:2791418/100‑1 (MQ=255)

           ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  1:1211387/100‑1 (MQ=255)

           ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  2:3647614/100‑1 (MQ=255)

            tCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctctt                                                                                        <  1:1234577/100‑2 (MQ=39)

             cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                       >  1:2067266/1‑100 (MQ=255)

             cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                       >  1:2522906/1‑100 (MQ=255)

              gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                      <  1:1133156/100‑1 (MQ=255)

              gCTAGTTGAACATAATCTGGCCATATGAACCCATGAATTAATTTCCATACTGGCAATTTTTTGCCAGACTGAAGAAGTCATAACTGCTATGACctcttca                                                                                      <  2:1322071/100‑4 (MQ=255)

               cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttcag                                                                                     <  1:3215526/100‑5 (MQ=38)

                tggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                    <  2:1001244/100‑1 (MQ=255)

                 ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                   >  1:649878/1‑100 (MQ=255)

                    ggAACGTAATCTGGCGGTATTGACCCTGTAATTAATTTGTATAGTGGCTATTTTTTGCCAGACTGAAGAGGTAATCCCAGTTATGACCTGTGTACTtata                                                                                <  2:3064425/100‑1 (MQ=255)

                    ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:88366/100‑1 (MQ=255)

                    ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:2686723/100‑1 (MQ=255)

                    ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:194478/100‑1 (MQ=255)

                      aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaac                                                                              <  1:6197/100‑1 (MQ=255)

                      aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaac                                                                              <  2:1937972/100‑1 (MQ=255)

                      aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                               <  1:909861/99‑1 (MQ=255)

                       aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             <  1:1084223/100‑1 (MQ=255)

                       aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             <  1:1402059/100‑1 (MQ=255)

                       aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             >  1:2796837/1‑100 (MQ=255)

                         gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaac                                                                           <  2:1426099/100‑1 (MQ=255)

                          tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                          <  1:3411111/100‑1 (MQ=255)

                          tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                          >  1:70518/1‑100 (MQ=255)

                               tGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTa                                                                     <  2:1621169/100‑1 (MQ=255)

                                ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                    >  2:3538718/1‑100 (MQ=255)

                                ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                    >  2:1848857/1‑100 (MQ=255)

                                 gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                   <  1:1848867/100‑1 (MQ=255)

                                  cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                  >  1:2435321/1‑100 (MQ=255)

                                     tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATCCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                               <  2:1867324/100‑1 (MQ=255)

                                     tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTAGCCAGAATGAAGAGGTAATACCAGATATGACCTCTGTACTTATAACAACAACGTacgctt                                                                >  2:1764771/1‑94 (MQ=255)

                                        tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGATCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                            <  2:2538044/100‑1 (MQ=255)

                                          aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                          <  2:2053694/100‑1 (MQ=255)

                                          aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                          >  2:460353/1‑100 (MQ=255)

                                          aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                          >  1:3434664/1‑100 (MQ=255)

                                             cTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                       <  1:3797534/100‑1 (MQ=255)

                                              tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGACGAGGTCATACCAGTTATGACCTCTTTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                      <  2:307248/100‑1 (MQ=255)

                                              tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                       <  1:80005/99‑1 (MQ=255)

                                              tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                      <  1:3679383/100‑1 (MQ=255)

                                                 aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   <  1:1863975/100‑1 (MQ=255)

                                                 aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   <  2:3122018/100‑1 (MQ=255)

                                                 aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   >  1:3344727/1‑100 (MQ=255)

                                                  attaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   >  1:3186893/1‑99 (MQ=255)

                                                    taattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                <  2:3194703/100‑1 (MQ=255)

                                                    taattTGCATAGTGGCAATTTTTTGCCAGACTGAACAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                <  2:784249/100‑1 (MQ=255)

                                                       tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                             >  2:2662168/1‑100 (MQ=255)

                                                       tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                             >  2:675077/1‑100 (MQ=255)

                                                         tGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                           >  2:1114287/1‑100 (MQ=255)

                                                          gCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAt                                          >  1:3811728/1‑100 (MQ=255)

                                                           cATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATc                                         >  1:2804607/1‑100 (MQ=255)

                                                            aTAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCa                                        >  1:259720/1‑100 (MQ=255)

                                                              aGTGGCAATGTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaaa                                      <  2:592218/100‑1 (MQ=255)

                                                               gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                     <  1:2206308/100‑1 (MQ=255)

                                                                tggcaaaaaattGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGt                                    >  2:2201346/11‑100 (MQ=255)

                                                                 ggCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                   >  2:1263847/1‑100 (MQ=255)

                                                                  gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAtagatcgga                                  >  1:1766486/1‑92 (MQ=255)

                                                                      ttttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                              <  1:3605956/100‑1 (MQ=255)

                                                                      ttttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                              <  2:3507251/100‑1 (MQ=255)

                                                                      tcttctGCCAGACTGAAGAGGTCATACCAGTTATTACCTCTGTACTTATAACCACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                              <  2:2309432/95‑1 (MQ=255)

                                                                       tttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGGTCGTGGa                             <  1:2001815/100‑1 (MQ=255)

                                                                        ttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                            <  2:307400/100‑1 (MQ=255)

                                                                         tttGCCAGACTGAAGAGGTTATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGata                           <  1:3085243/100‑1 (MQ=255)

                                                                            gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCCAAGTTCGTGGATatca                        >  2:3353245/1‑100 (MQ=255)

                                                                            gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                        <  2:3425003/100‑1 (MQ=255)

                                                                            gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                        >  1:3461408/1‑100 (MQ=255)

                                                                             ccAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcat                       >  1:2462925/1‑100 (MQ=255)

                                                                                gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                    >  1:1196547/1‑100 (MQ=255)

                                                                                gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                    >  1:3583383/1‑100 (MQ=255)

                                                                                 aCTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCg                   >  2:1622196/1‑100 (MQ=255)

                                                                                  gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                  >  2:3641837/5‑100 (MQ=38)

                                                                                    gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                >  2:774391/1‑100 (MQ=255)

                                                                                    gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                <  1:1959207/100‑1 (MQ=255)

                                                                                    acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGGTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                >  1:1281328/3‑100 (MQ=38)

                                                                                     cagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt               >  1:3485520/2‑100 (MQ=39)

                                                                                      agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGtt              >  2:3560450/1‑100 (MQ=255)

                                                                                       gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGttt             >  1:560448/1‑100 (MQ=255)

                                                                                         ggTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc           >  1:1077953/1‑100 (MQ=255)

                                                                                         ggTCATACCAGTTATGACCTCTGTACGTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc           >  2:754129/1‑100 (MQ=255)

                                                                                          gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc          <  1:2893614/100‑1 (MQ=255)

                                                                                           tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa         <  2:2716050/100‑1 (MQ=255)

                                                                                            cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg        >  2:3605996/1‑100 (MQ=255)

                                                                                              tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  2:2956608/100‑1 (MQ=255)

                                                                                              tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  1:2624421/100‑1 (MQ=255)

                                                                                              tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      >  2:3739118/1‑100 (MQ=255)

                                                                                              tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      >  2:1766901/1‑100 (MQ=255)

                                                                                              tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  2:3438783/100‑1 (MQ=255)

                                                                                                ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg    >  1:3820075/1‑100 (MQ=255)

                                                                                                ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg    <  1:1482720/100‑1 (MQ=255)

                                                                                                 cAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGt   >  1:749585/1‑100 (MQ=255)

                                                                                                  aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTc  <  1:155229/100‑1 (MQ=255)

GGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTC  >  NC_000913/464258‑464455

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: `‑`