BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	112,238	A→G	V121A (GTG→GCG)	yacF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	112,238	0	A	G	99.6%	872.7 / ‑4.9	265	V121A (GTG→GCG)	yacF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (134/130); total (135/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 269 total aligned reads displayed.

CGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112149‑112327

  cgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCaga                                                                                      >  2:921679/1‑93 (MQ=255)

      gcgcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCaa                                                                                  <  1:1269250/93‑1 (MQ=255)

        gcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATc                                                                                <  2:3656142/93‑1 (MQ=255)

        gcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATc                                                                                <  2:3331801/93‑1 (MQ=255)

         ccTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCa                                                                               >  2:4448353/1‑93 (MQ=255)

         ccTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCa                                                                               >  2:275380/1‑93 (MQ=255)

          cTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaa                                                                              <  2:5929357/93‑1 (MQ=255)

          cTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaa                                                                              >  2:5344714/1‑93 (MQ=255)

           tGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaaa                                                                             <  2:4888673/93‑1 (MQ=255)

           tGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaaa                                                                             >  1:221290/1‑93 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  2:1176011/93‑1 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  2:6629193/93‑1 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  1:4622521/93‑1 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  1:2507653/93‑1 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  1:1426892/93‑1 (MQ=255)

               gTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGTAATCAAACGAt                                                                         <  2:1502791/93‑1 (MQ=255)

                tAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATc                                                                        <  1:269848/93‑1 (MQ=255)

                tAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATc                                                                        <  2:311736/93‑1 (MQ=255)

                tAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATc                                                                        <  2:3844480/93‑1 (MQ=255)

                 aGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCt                                                                       <  1:2607579/93‑1 (MQ=255)

                  gATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCtt                                                                      >  1:1015370/1‑93 (MQ=255)

                   aTGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTc                                                                     >  2:5521763/1‑93 (MQ=255)

                   aTGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTc                                                                     <  2:2510119/93‑1 (MQ=255)

                    tGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCa                                                                    <  1:740164/93‑1 (MQ=255)

                    tGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCa                                                                    >  1:6071642/1‑93 (MQ=255)

                     gCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCACCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCAc                                                                   <  2:3410153/93‑1 (MQ=255)

                       aGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGc                                                                 <  2:1238563/93‑1 (MQ=255)

                       aGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGc                                                                 >  2:3538938/1‑93 (MQ=255)

                        gCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCa                                                                >  1:489887/1‑93 (MQ=255)

                        gCCAAATGTGCAATGTAGGCAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGACAGAGCAATCAAACGATCTTCACGCa                                                                >  1:5206880/1‑93 (MQ=255)

                          cAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGa                                                              >  2:6250510/1‑93 (MQ=255)

                            aaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaaa                                                            >  2:4763373/1‑93 (MQ=255)

                             aTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAAt                                                           <  1:89679/93‑1 (MQ=255)

                             aTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAAt                                                           >  1:6713108/1‑93 (MQ=255)

                                tgCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGc                                                        >  1:1753959/1‑93 (MQ=255)

                                   aaTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCg                                                     >  1:4981440/1‑93 (MQ=255)

                                   aaTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCg                                                     >  2:6183639/1‑93 (MQ=255)

                                   aaTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCg                                                     >  2:6100458/1‑93 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   >  2:4992955/1‑93 (MQ=255)

                                      gTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATa                                                  >  1:5930152/1‑93 (MQ=255)

                                      gTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATa                                                  <  1:3878932/93‑1 (MQ=255)

                                          gTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcg                                              <  2:2203206/93‑1 (MQ=255)

                                           tAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcgg                                             >  1:5877122/1‑93 (MQ=255)

                                             aaTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcg                                           <  1:6144609/93‑1 (MQ=255)

                                             aaTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcg                                           >  1:5429362/1‑93 (MQ=255)

                                               tCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgcg                                         >  1:6166697/1‑93 (MQ=255)

                                               tCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgcg                                         >  1:6490102/1‑93 (MQ=255)

                                               acAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgcg                                         <  1:4445081/92‑1 (MQ=255)

                                                cAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATCcgcggcgcgg                                        >  2:1267244/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  1:1729628/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       <  1:2002298/93‑1 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  1:1106901/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  1:3117594/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  1:362909/1‑93 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCTAGAGCAATCAAACCATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      >  2:168003/1‑93 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      <  2:86891/93‑1 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      >  2:2362317/1‑93 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      >  2:4838491/1‑93 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      >  2:3484094/1‑93 (MQ=255)

                                                   aGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaaa                                     >  2:466105/1‑93 (MQ=255)

                                                     cTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGGCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAAtt                                   >  2:6645798/1‑93 (MQ=255)

                                                      tGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTa                                  <  2:1576511/93‑1 (MQ=255)

                                                      tGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTa                                  >  1:4286432/1‑93 (MQ=255)

                                                       gCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTaa                                 <  1:6404475/93‑1 (MQ=255)

                                                       gCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTaa                                 >  2:3196242/1‑93 (MQ=255)

                                                         aaCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATa                               <  1:5057843/93‑1 (MQ=255)

                                                          acacCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              <  1:2498759/93‑1 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              >  2:2835695/1‑93 (MQ=255)

                                                             gccgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCt                           >  1:3005042/1‑93 (MQ=255)

                                                             gccgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCt                           <  2:5502878/93‑1 (MQ=255)

                                                               cgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTcc                         <  1:17939/93‑1 (MQ=255)

                                                                gccTGGGAAGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTccc                        >  1:6532038/1‑93 (MQ=255)

                                                                 ccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTcccc                       <  1:482751/93‑1 (MQ=255)

                                                                 ccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTcccc                       >  2:5006762/1‑93 (MQ=255)

                                                                  cTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccg                      >  1:1553960/1‑93 (MQ=255)

                                                                  cTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccg                      >  2:5175681/1‑93 (MQ=255)

                                                                  cTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccg                      >  2:5125782/1‑93 (MQ=255)

                                                                   tGGGATGCTCAGTCGCTGAGGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCAATACGCGGCGCGGAAATTAATACGCTCCccgc                     <  1:1451323/93‑1 (MQ=255)

                                                                   tGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgc                     >  1:6273402/1‑93 (MQ=255)

                                                                    gggATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgcc                    >  1:2846477/1‑93 (MQ=255)

                                                                    gggATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgcc                    >  2:4027888/1‑93 (MQ=255)

                                                                    gggATGCTCAGTCGCTGACGCGACAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCTGCGCGGAAAATGATACGCTCCccgcc                    >  1:3755970/1‑93 (MQ=255)

                                                                       aTGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCt                 <  2:5319973/93‑1 (MQ=255)

                                                                         gCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCttt               <  1:5507118/93‑1 (MQ=255)

                                                                         gCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCttt               <  1:4350504/93‑1 (MQ=255)

                                                                         gCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCttt               <  2:5599388/93‑1 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  1:4494817/1‑93 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  2:4499872/1‑93 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  2:6391477/1‑93 (MQ=255)

                                                                             aGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAc           <  2:5063926/93‑1 (MQ=255)

                                                                             aGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAc           <  1:3346499/93‑1 (MQ=255)

                                                                             aGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAc           >  1:3299711/1‑93 (MQ=255)

                                                                                cGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgc        <  2:6689478/93‑1 (MQ=255)

                                                                                cGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgc        >  2:1792520/1‑93 (MQ=255)

                                                                                cGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgc        <  2:2019801/93‑1 (MQ=255)

                                                                                 gCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgct       >  2:4276185/1‑93 (MQ=255)

                                                                                 gCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgct       <  1:2144505/93‑1 (MQ=255)

                                                                                    gACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGaa    <  2:2641714/93‑1 (MQ=255)

                                                                                     aCGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAAt   <  1:4377422/93‑1 (MQ=255)

                                                                                     aCGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAAt   <  1:5684448/93‑1 (MQ=255)

CGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112149‑112327

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