BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,519	T→C	T122A (ACG→GCG)	yddE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,519	0	T	C	99.6%	795.1 / ‑5.1	249	T122A (ACG→GCG)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (130/118); total (131/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 258 total aligned reads displayed.

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608

aGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAc                                                                                         >  2:3879428/1‑93 (MQ=255)

  ccTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCt                                                                                       <  2:5172532/93‑1 (MQ=255)

   cTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCtt                                                                                      >  1:4893861/1‑93 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    >  1:4522130/1‑93 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    <  2:4944633/93‑1 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    <  2:8458222/93‑1 (MQ=255)

      gCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGtt                                                                                   >  1:179136/1‑93 (MQ=255)

        aGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTcc                                                                                 <  1:6277272/93‑1 (MQ=255)

         gAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCa                                                                                >  2:7157806/1‑93 (MQ=255)

          aaTATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAg                                                                               <  1:5609027/93‑1 (MQ=255)

          aaTATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAg                                                                               <  2:8317780/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCTCCTTCCAGCGGTGGCTCAAATCCCGGCGCACCTTGTTCCAGc                                                                              <  2:2533018/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGc                                                                              <  1:923109/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGc                                                                              >  1:2510511/1‑93 (MQ=255)

            tatCGTCCTCAGTTAGATGATGCTCGTTGATAATCTCCGCTCGTGTTTCACCTTCCAGTGGTGGCTCAAAGCCCGTCGCACCTTGTTCCATCg                                                                             >  1:8940645/1‑93 (MQ=255)

            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                             >  1:3701996/1‑93 (MQ=255)

               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                          <  2:904770/93‑1 (MQ=255)

                gTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAAt                                                                         >  1:8667279/1‑93 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                        <  2:4457307/93‑1 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                        <  1:9657714/93‑1 (MQ=255)

                  ccTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAc                                                                       <  2:2516948/93‑1 (MQ=255)

                   cTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACg                                                                      >  1:6550177/1‑93 (MQ=255)

                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     <  2:4152026/93‑1 (MQ=255)

                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGAt                                                                    <  2:1723221/93‑1 (MQ=255)

                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGCGAAATACGAt                                                                    <  1:3650947/93‑1 (MQ=255)

                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:6633193/1‑93 (MQ=255)

                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:6695418/1‑93 (MQ=255)

                       gTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                  <  2:4385253/93‑1 (MQ=255)

                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 <  2:1188644/93‑1 (MQ=255)

                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 >  2:7818327/1‑93 (MQ=255)

                         gagaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACTATAatc                                                                >  1:7504659/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               >  2:4262074/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               <  1:2564877/93‑1 (MQ=255)

                           gaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcat                                                              >  2:5373350/1‑93 (MQ=255)

                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            <  2:9363140/93‑1 (MQ=255)

                              taaGCGCGTTGATAATCTCCGCCTGTGTTTCACCGTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCGGCGAAATACGATAATCATCGt                                                           <  1:3474488/92‑1 (MQ=255)

                              gAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGt                                                           <  2:9902022/93‑1 (MQ=255)

                               aaGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGtt                                                          <  2:7233092/93‑1 (MQ=255)

                               aaGCGCGTTGATAATCGCCCCACGGGCTTCACCTTCCAGCGGGGGCTCAAAGCCGGGGGCACCTTGTTCACGCGAAATACGTTATTCATAGtt                                                          >  2:848418/1‑93 (MQ=255)

                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         <  1:7656166/93‑1 (MQ=255)

                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         >  1:3842309/1‑93 (MQ=255)

                                aGCGCGTTAATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         <  2:6696697/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                      <  2:5304499/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                      <  1:7083039/93‑1 (MQ=255)

                                    cgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCt                                                     <  2:6555757/93‑1 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  1:2056802/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  2:5751460/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  2:7363229/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   <  1:5928081/93‑1 (MQ=255)

                                        gATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                 >  1:7882945/1‑93 (MQ=255)

                                        gATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                 <  2:7636255/93‑1 (MQ=255)

                                          tAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCTCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCgt                                               >  1:5274776/1‑92 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGCTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                              >  1:6867062/1‑93 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:6987802/93‑1 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:7454279/93‑1 (MQ=255)

                                               cccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                          <  1:6840248/92‑1 (MQ=255)

                                                 cgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAc                                        >  2:2692894/1‑93 (MQ=255)

                                                    aCGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgcg                                     >  1:727894/1‑93 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgcca                                    >  2:1939619/1‑91 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:1605234/93‑1 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:8912745/93‑1 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    >  2:4731064/1‑93 (MQ=255)

                                                      gtgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAt                                   >  2:5412496/1‑93 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  >  1:6397984/1‑93 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  <  2:7957714/93‑1 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  >  1:8782451/1‑93 (MQ=255)

                                                        gtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGt                                 >  2:1214351/1‑93 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                <  1:9798985/93‑1 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                <  2:3123535/93‑1 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                >  1:4406163/1‑93 (MQ=255)

                                                              ccTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCt                           <  1:2308866/93‑1 (MQ=255)

                                                               cTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTg                          >  2:9692743/1‑93 (MQ=255)

                                                                ttCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGc                         >  2:3690960/1‑93 (MQ=255)

                                                                 tCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                        >  1:2798815/1‑93 (MQ=255)

                                                                 tCCAGCGGGGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                        <  1:2695225/93‑1 (MQ=255)

                                                                   cAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAg                      >  1:3211488/1‑93 (MQ=255)

                                                                     gCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                    >  2:664221/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   >  1:7155764/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   >  1:9112649/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   <  2:5882691/93‑1 (MQ=255)

                                                                       ggTGGCTCAAAGCCCTGCGCGCCTTGTTCCAGCAAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAt                  <  2:8188586/93‑1 (MQ=255)

                                                                       ggTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAt                  >  1:1647252/1‑93 (MQ=255)

                                                                        gTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATg                 >  1:6774947/1‑93 (MQ=255)

                                                                        gTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATg                 >  1:3370984/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  2:2795712/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  2:1442592/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  1:467850/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGgtt              >  2:2405970/1‑93 (MQ=255)

                                                                              cAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGcc           >  1:9563859/1‑93 (MQ=255)

                                                                              cAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGcc           >  2:3620577/1‑93 (MQ=255)

                                                                               aaaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCa          <  2:344021/93‑1 (MQ=255)

                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg         >  1:9486087/1‑93 (MQ=255)

                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg         <  1:829393/93‑1 (MQ=255)

                                                                                 aGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGa        <  1:3878966/93‑1 (MQ=255)

                                                                                    ccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCg     >  1:4407808/1‑93 (MQ=255)

                                                                                    ccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCg     <  1:3663094/93‑1 (MQ=255)

                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt    >  1:839197/1‑93 (MQ=255)

                                                                                      ggCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg   <  1:5107383/93‑1 (MQ=255)

                                                                                      cgcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg   <  2:3195599/92‑1 (MQ=255)

                                                                                       gcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGc  <  1:7912434/93‑1 (MQ=255)

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608

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