BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,053,515	T→C	G619G (GGA→GGG)	torS ←	TorS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,515	0	T	C	98.9%	539.0 / ‑3.9	177	G619G (GGA→GGG)	torS	TorS
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (87/88); total (89/88)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 194 total aligned reads displayed.

CGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCT  >  NC_000913/1053426‑1053604

cGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGcccc                                                                                        <  2:4842291/93‑1 (MQ=255)

 gTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCa                                                                                       >  2:4784662/1‑93 (MQ=255)

 gTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCa                                                                                       >  2:1165544/1‑93 (MQ=255)

 gTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCa                                                                                       <  1:7335657/93‑1 (MQ=255)

  tAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAg                                                                                      >  1:1575221/1‑93 (MQ=255)

  tAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAg                                                                                      <  2:5579254/93‑1 (MQ=255)

   aaaCAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGc                                                                                     <  1:1796389/93‑1 (MQ=255)

    aaCAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGcc                                                                                    >  1:4305028/1‑93 (MQ=255)

     aCAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTcagccccagccc                                                                                   <  1:1407518/93‑1 (MQ=255)

      cAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGACAGCCACAGACCg                                                                                  >  1:107847/1‑93 (MQ=255)

       aaaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCgg                                                                                 >  1:6910058/1‑93 (MQ=255)

       aaaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCgg                                                                                 <  2:6204369/93‑1 (MQ=255)

        aaaaCATCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGt                                                                                <  2:2401321/93‑1 (MQ=255)

        aaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCCGCCCGGt                                                                                >  2:641393/1‑93 (MQ=255)

        aaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGt                                                                                >  2:6379309/1‑93 (MQ=255)

        aaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGt                                                                                >  1:5281251/1‑93 (MQ=255)

         aaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTg                                                                               <  1:2106673/93‑1 (MQ=255)

           aCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTgcc                                                                             <  1:680277/93‑1 (MQ=255)

            cAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTgccg                                                                            >  1:7275611/1‑93 (MQ=255)

             aGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTgccgc                                                                           >  1:2539125/1‑93 (MQ=255)

             aGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTgccgc                                                                           >  1:6002919/1‑93 (MQ=255)

             aGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTgccgc                                                                           <  1:7612079/93‑1 (MQ=255)

              gCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTgccgcc                                                                          >  2:4511319/1‑93 (MQ=255)

                ttCCAACCTCCGGCATGCGGGTGGCGCTCAGATCGCCGCCCATCGCCTGCGGCAGACTGCTACTGATAGGCAGACCCAGCCCCGTGCcgccgc                                                                        >  2:6410797/1‑93 (MQ=255)

                 tCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGCGCCGCcgcg                                                                       >  2:3312142/1‑93 (MQ=255)

                   cAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGtt                                                                     <  1:2768855/93‑1 (MQ=255)

                    aaCCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGttt                                                                    >  1:3290168/1‑93 (MQ=255)

                    aaCCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGttt                                                                    >  2:2234994/1‑93 (MQ=255)

                      ccTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGc                                                                  <  2:7123161/93‑1 (MQ=255)

                       cTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGcc                                                                 <  2:6795237/93‑1 (MQ=255)

                       cTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGcc                                                                 >  1:2866733/1‑93 (MQ=255)

                       cTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGcc                                                                 >  2:7501785/1‑93 (MQ=255)

                          cGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCt                                                              <  1:3373667/93‑1 (MQ=255)

                           ggCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCtt                                                             >  1:1098768/1‑93 (MQ=255)

                              ggggtggtgtCGCTCAGTGCGCCGCCGAGCGCCGGGTCCAGACGGCTACTGAGAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          <  1:2204559/89‑1 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          >  1:6529140/1‑93 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          >  1:6135501/1‑93 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          >  1:445879/1‑93 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          <  1:5419086/93‑1 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          <  1:397584/93‑1 (MQ=255)

                              gTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAcc                                                          <  1:513376/93‑1 (MQ=255)

                               tGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCt                                                         <  2:7052930/93‑1 (MQ=255)

                               tGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCt                                                         >  1:1766766/1‑93 (MQ=255)

                               tGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCt                                                         >  2:2591820/1‑93 (MQ=255)

                                 cTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGc                                                       >  2:4792931/1‑93 (MQ=255)

                                   ggtggCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGcac                                                     >  1:3826814/1‑93 (MQ=255)

                                   ggtggCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGcac                                                     >  1:558168/1‑93 (MQ=255)

                                      ggCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACaaa                                                  <  1:6677679/93‑1 (MQ=255)

                                      ggCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACaaa                                                  <  1:5852902/93‑1 (MQ=255)

                                       gcgcTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAAt                                                 <  2:4391960/93‑1 (MQ=255)

                                           tCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCt                                             >  2:6621510/1‑93 (MQ=255)

                                           tCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCt                                             >  2:2140247/1‑93 (MQ=255)

                                             aGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTgg                                           >  1:237933/1‑93 (MQ=255)

                                              gTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGa                                          >  1:6236696/1‑93 (MQ=255)

                                               tgccccgccCATCGCCTGGGGCAGACGGCTACGGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGaa                                         <  1:4639684/89‑1 (MQ=255)

                                                tCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAg                                        <  1:6804303/93‑1 (MQ=255)

                                                  gccgccCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGAt                                      <  2:5578846/93‑1 (MQ=255)

                                                  gccgccCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGAt                                      <  1:2831879/93‑1 (MQ=255)

                                                  gccgccCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGAt                                      <  1:6730808/93‑1 (MQ=255)

                                                   ccgccCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGAtt                                     <  1:7133208/93‑1 (MQ=255)

                                                     gccCATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTc                                   <  1:4683039/93‑1 (MQ=255)

                                                     ccccATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTc                                   <  1:1530154/92‑1 (MQ=255)

                                                     ccccATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTc                                   <  2:2451796/92‑1 (MQ=255)

                                                        cATCGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGc                                >  2:3265833/1‑93 (MQ=255)

                                                           cGCCTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAg                             <  2:2087622/93‑1 (MQ=255)

                                                             ccTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGtt                           <  2:876451/93‑1 (MQ=255)

                                                             ccTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGtt                           <  2:6679856/93‑1 (MQ=255)

                                                             ccTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGtt                           <  1:7151747/93‑1 (MQ=255)

                                                              cTGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGttt                          >  2:4002301/1‑93 (MQ=255)

                                                               tGGGCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTc                         >  1:3964313/1‑93 (MQ=255)

                                                                gggCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTcg                        >  2:7860908/1‑93 (MQ=255)

                                                                 ggCCAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTcgc                       >  2:3044808/1‑93 (MQ=255)

                                                                   ccAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCgg                     <  1:2809560/93‑1 (MQ=255)

                                                                   ccAGACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCgg                     <  2:1265465/93‑1 (MQ=255)

                                                                     aGACGGCTACTGATAGTCAGCCCCAGCCCGGGGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGa                   >  1:1782044/1‑93 (MQ=255)

                                                                      gACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGGCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGAt                  <  1:1327305/93‑1 (MQ=255)

                                                                      gACGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGAt                  >  1:5731996/1‑93 (MQ=255)

                                                                       aCGGCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACTTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATc                 >  2:7010818/1‑93 (MQ=255)

                                                                          gCTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAAt              <  2:3747297/93‑1 (MQ=255)

                                                                           cTACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATa             <  1:2663320/93‑1 (MQ=255)

                                                                           cTACTGAGAGGCAGCCCCAGCGCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATa             <  1:5263784/93‑1 (MQ=255)

                                                                            tACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATAc            >  1:893160/1‑93 (MQ=255)

                                                                            tACTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATAc            <  2:6081595/93‑1 (MQ=255)

                                                                             aCTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATAcc           <  2:2671953/93‑1 (MQ=255)

                                                                              cTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCg          <  1:5112611/93‑1 (MQ=255)

                                                                              cTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCg          >  2:1532471/1‑93 (MQ=255)

                                                                              cTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCg          <  1:2169876/93‑1 (MQ=255)

                                                                              cTGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTAGTCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCg          >  2:4957115/1‑93 (MQ=255)

                                                                               tGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGc         <  1:3640367/93‑1 (MQ=255)

                                                                               tGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGc         <  1:1019286/93‑1 (MQ=255)

                                                                               tGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGc         >  2:4242061/1‑93 (MQ=255)

                                                                               tGATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGc         <  1:1643288/93‑1 (MQ=255)

                                                                                gATAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCa        <  2:2853666/93‑1 (MQ=255)

                                                                                 aTAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAg       >  2:5859138/1‑93 (MQ=255)

                                                                                  tAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGc      <  2:4903159/93‑1 (MQ=255)

                                                                                  tAGTCAGCCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGc      <  2:1795669/93‑1 (MQ=255)

                                                                                   agtcagccCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGcc     <  1:76901/93‑1 (MQ=255)

                                                                                   agtcagccCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGcc     <  2:2285520/93‑1 (MQ=255)

                                                                                     tcagccCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGc   <  1:3717857/93‑1 (MQ=255)

                                                                                      cagccccagcccGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCt  <  1:1954401/93‑1 (MQ=255)

CGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCT  >  NC_000913/1053426‑1053604

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