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.7%	924.3 / ‑5.3	291	V121A (GTG→GCG)	yacF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (161/129); total (161/130)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 298 total aligned reads displayed.

TCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112148‑112327

tCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGcc                                                                                         <  1:3218949/93‑1 (MQ=255)

 cgcgCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCa                                                                                        <  1:2271639/93‑1 (MQ=255)

  gcgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCag                                                                                       <  2:6831457/93‑1 (MQ=255)

  gcgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCag                                                                                       >  2:3457880/1‑93 (MQ=255)

   cgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCaga                                                                                      <  2:852695/93‑1 (MQ=255)

   cgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCaga                                                                                      <  1:4263789/93‑1 (MQ=255)

   cgcTGCGCCTGGGGAAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAATCTGCAACAGCCGCCTGGGATGCGCAGTCGCTGACGCGCCaga                                                                                      <  2:196747/93‑1 (MQ=255)

    gcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCagag                                                                                     <  2:7514687/93‑1 (MQ=255)

    gcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCagag                                                                                     >  1:1473157/1‑93 (MQ=255)

     cTGCGCGTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGc                                                                                    >  1:7750907/1‑93 (MQ=255)

     cTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGc                                                                                    <  1:2802803/93‑1 (MQ=255)

      tGCGCCTGGGGTAGATGCAGCCAAAAGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCa                                                                                   <  1:5936377/93‑1 (MQ=255)

        cgcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAAt                                                                                 <  1:4924061/93‑1 (MQ=255)

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

             tgggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGGCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  1:379306/92‑1 (MQ=255)

             ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            >  1:7319983/1‑93 (MQ=255)

              gggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACg                                                                           >  2:5890149/1‑93 (MQ=255)

              gggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACg                                                                           >  2:2463461/1‑93 (MQ=255)

                gTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGAt                                                                         <  1:2104477/93‑1 (MQ=255)

                gTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGAt                                                                         <  1:6769136/93‑1 (MQ=255)

                gTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGAt                                                                         <  1:3488473/93‑1 (MQ=255)

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

                 tAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATc                                                                        >  2:5907619/1‑93 (MQ=255)

                  aGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCt                                                                       >  1:1347973/1‑93 (MQ=255)

                  aGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCt                                                                       >  1:5201672/1‑93 (MQ=255)

                   gATGCAGCCAAATGTGCAATGTAGGCAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCtt                                                                      <  2:825659/93‑1 (MQ=255)

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

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

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

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

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

                          ccAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAg                                                               >  2:1170661/1‑93 (MQ=255)

                            aaaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaa                                                             >  2:5403741/1‑93 (MQ=255)

                            aaaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaa                                                             >  2:3525868/1‑93 (MQ=255)

                            aaaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaa                                                             >  1:1666501/1‑93 (MQ=255)

                            aaaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaa                                                             >  2:6068930/1‑93 (MQ=255)

                            aaaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaa                                                             >  2:444500/1‑93 (MQ=255)

                               tgtgCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAAtt                                                          <  1:7795205/93‑1 (MQ=255)

                                  gCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGcc                                                       >  2:6619869/1‑93 (MQ=255)

                                  gCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGcc                                                       >  2:4310279/1‑93 (MQ=255)

                                    aaTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCg                                                     <  1:4667592/93‑1 (MQ=255)

                                     aTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGa                                                    <  2:3247304/93‑1 (MQ=255)

                                     aTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGa                                                    <  1:6725890/93‑1 (MQ=255)

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

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

                                         aGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcg                                                >  2:3132063/1‑93 (MQ=255)

                                          ggTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgc                                               <  2:3913313/93‑1 (MQ=255)

                                          ggTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTACACGCAGAAATTGCCCGATAcgc                                               >  2:6165673/1‑93 (MQ=255)

                                           gTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcg                                              >  2:3345523/1‑93 (MQ=255)

                                             aaaTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggc                                            <  1:5296314/93‑1 (MQ=255)

                                             aaaTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggc                                            <  2:7737959/93‑1 (MQ=255)

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

                                               aTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgg                                          >  1:5864674/1‑92 (MQ=255)

                                               aTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgc                                          >  2:1899936/1‑93 (MQ=255)

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

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

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

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

                                                     gCTGCAAGAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAAt                                    <  2:1984291/93‑1 (MQ=255)

                                                     gCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAAt                                    <  1:1345526/93‑1 (MQ=255)

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

                                                         cagcAGCCACCTGGGAGGCTCAGTCGCTGCCGCGCCAGAGCAATCAAACGATCTTCACGGAGAAATCGCCCGATACGCGGCGAGGAAATTAAt                                <  2:3066577/90‑1 (MQ=255)

                                                         cAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAAt                                >  2:7600894/1‑93 (MQ=255)

                                                         cAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAAt                                >  2:2668771/1‑93 (MQ=255)

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

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

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

                                                          aaCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATa                               >  1:4561210/1‑93 (MQ=255)

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

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

                                                             aGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGc                            >  2:2141093/1‑93 (MQ=255)

                                                             aGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGc                            <  1:2400947/93‑1 (MQ=255)

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

                                                                cgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTcc                         >  1:2448957/1‑93 (MQ=255)

                                                                cgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTcc                         >  2:6713398/1‑93 (MQ=255)

                                                                 gccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTccc                        >  2:6583872/1‑93 (MQ=255)

                                                                 cccTGAGATGGTCAGTCGCTGACGCGCCAGAGCAATCAGACGATCTTCACGCAGAAATTGCCGGATACGCGGCGCGGAAATTAATACGCTccc                        <  2:4975122/92‑1 (MQ=255)

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

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

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

                                                                       gATGCTCAGTCGCTGCCGCGCCAGAGCAATAAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                  <  2:3741696/93‑1 (MQ=255)

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

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

                                                                         tGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCtt                >  1:5616045/1‑93 (MQ=255)

                                                                          gCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCttt               >  1:6870580/1‑93 (MQ=255)

                                                                           cTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCtttt              <  2:5205733/93‑1 (MQ=255)

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

                                                                            tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             <  2:2783191/93‑1 (MQ=255)

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

                                                                             cAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTaa            <  1:3169144/93‑1 (MQ=255)

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

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

                                                                               gTCGCTGACGCGCTAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGTGCGTAAATTTATACGCTCCCCGCTGCTTTTAAcg          >  2:6995206/1‑92 (MQ=255)

                                                                               gTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAct          >  1:5707991/1‑93 (MQ=255)

                                                                                tCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActg         >  1:704949/1‑93 (MQ=255)

                                                                                     gACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGaa    >  2:4063130/1‑93 (MQ=255)

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

                                                                                     gACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGaa    >  2:462875/1‑93 (MQ=255)

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

                                                                                      aCGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAAt   >  2:5778527/1‑93 (MQ=255)

TCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112148‑112327

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