BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,225,456	T→C	T5T (ACT→ACC)	ygjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,225,456	0	T	C	99.3%	884.3 / ‑3.3	273	T5T (ACT→ACC)	ygjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (187/84); total (188/85)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 279 total aligned reads displayed.

ATAGTTTCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACTGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACAC  >  NC_000913/3225366‑3225545

aTAGTTTCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGc                                                                                         >  1:5128310/1‑93 (MQ=255)

aTAGTTTCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGc                                                                                         >  1:1860671/1‑93 (MQ=255)

 tAGTTTCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCa                                                                                        >  1:3670282/1‑93 (MQ=255)

    tttCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTGTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCa                                                                                     <  2:4327007/93‑1 (MQ=255)

     ttCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCAt                                                                                    <  2:6006893/93‑1 (MQ=255)

       cTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGc                                                                                  >  2:1254517/1‑93 (MQ=255)

       cTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGc                                                                                  <  1:6465401/93‑1 (MQ=255)

        tGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCa                                                                                 >  1:5271278/1‑93 (MQ=255)

        tGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCa                                                                                 >  1:3311367/1‑93 (MQ=255)

         gCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCag                                                                                >  1:2476077/1‑93 (MQ=255)

         gCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCaa                                                                                >  2:627704/1‑92 (MQ=255)

         gCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACAATGctg                                                                                >  1:808277/1‑91 (MQ=255)

          ccTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCaga                                                                               >  2:5718431/1‑93 (MQ=255)

            tttATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGc                                                                             >  2:6391706/1‑93 (MQ=255)

             ttATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCa                                                                            >  1:3423677/1‑93 (MQ=255)

               atatTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCAtt                                                                          >  2:2059864/1‑93 (MQ=255)

               atatTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCAtt                                                                          >  2:3680468/1‑93 (MQ=255)

               atatTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCAtt                                                                          >  1:3812030/1‑93 (MQ=255)

                     ttATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGc                                                                    >  2:55940/1‑93 (MQ=255)

                     ttATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGc                                                                    <  1:4631200/93‑1 (MQ=255)

                     ttATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGc                                                                    >  1:3309240/1‑93 (MQ=255)

                      tataATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGct                                                                   <  1:4017748/93‑1 (MQ=255)

                      tataATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGct                                                                   >  2:6159266/1‑93 (MQ=255)

                       ataATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGctc                                                                  >  1:1555555/1‑93 (MQ=255)

                       ataATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGctc                                                                  >  2:139265/1‑93 (MQ=255)

                        taATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGctct                                                                 <  2:4769773/93‑1 (MQ=255)

                        taATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGctct                                                                 <  2:3071960/93‑1 (MQ=255)

                        taATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGctct                                                                 >  1:2655857/1‑93 (MQ=255)

                         aaTTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTg                                                                >  2:392196/1‑93 (MQ=255)

                          aTTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGc                                                               >  2:5340050/1‑93 (MQ=255)

                          aTTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGc                                                               <  2:2157323/93‑1 (MQ=255)

                           ttACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGcc                                                              <  2:68682/93‑1 (MQ=255)

                           ttACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGcc                                                              >  1:5755943/1‑93 (MQ=255)

                              cAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGtt                                                           <  1:5577568/93‑1 (MQ=255)

                              cAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGtt                                                           <  2:3489880/93‑1 (MQ=255)

                              cAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGtt                                                           <  1:6002324/93‑1 (MQ=255)

                               aaCGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGttt                                                          >  1:1647041/1‑93 (MQ=255)

                               aaCGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGttt                                                          <  2:4132949/93‑1 (MQ=255)

                                  gATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGc                                                       <  1:3009991/93‑1 (MQ=255)

                                  gATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGc                                                       >  2:2333674/1‑93 (MQ=255)

                                   aTAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCt                                                      <  1:6272267/93‑1 (MQ=255)

                                    tAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTa                                                     >  1:5381765/1‑93 (MQ=255)

                                    tAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTa                                                     >  1:2094241/1‑93 (MQ=255)

                                      aaaGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACg                                                   >  2:2776933/1‑93 (MQ=255)

                                      aaaGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACg                                                   >  1:4929614/1‑93 (MQ=255)

                                      aaaGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACg                                                   >  1:4160366/1‑93 (MQ=255)

                                      aaaGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACg                                                   >  2:4727107/1‑93 (MQ=255)

                                        aGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGcc                                                 >  2:4437325/1‑93 (MQ=255)

                                        aGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGcc                                                 >  1:6396861/1‑93 (MQ=255)

                                        aGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGcc                                                 >  1:2363240/1‑93 (MQ=255)

                                          gCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCtt                                               <  2:3428599/93‑1 (MQ=255)

                                            tGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCtttt                                             <  1:1025377/93‑1 (MQ=255)

                                            tGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCtttt                                             >  1:4365926/1‑93 (MQ=255)

                                              tACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCt                                           <  2:1748085/93‑1 (MQ=255)

                                              tACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCt                                           >  2:5529027/1‑93 (MQ=255)

                                              tACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCt                                           >  2:6088979/1‑93 (MQ=255)

                                              tACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCt                                           >  2:2975676/1‑93 (MQ=255)

                                               aCTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTg                                          >  1:222915/1‑93 (MQ=255)

                                               aCTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTg                                          <  2:266736/93‑1 (MQ=255)

                                                cTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGc                                         >  1:1017362/1‑93 (MQ=255)

                                                 tttttCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCg                                        >  2:2557578/1‑93 (MQ=255)

                                                     tCTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAg                                    <  2:1818037/93‑1 (MQ=255)

                                                      cTTTAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAgg                                   <  1:5340481/93‑1 (MQ=255)

                                                       tttAGCTCGTGGATTAACAAAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGga                                  >  1:6571985/1‑93 (MQ=255)

                                                       tttAGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGga                                  >  2:6660326/1‑93 (MQ=255)

                                                          aGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGagaa                               >  1:4688398/1‑93 (MQ=255)

                                                          aGCTCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGagaa                               <  2:3177890/93‑1 (MQ=255)

                                                             tCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAgaagaa                            >  1:3554849/1‑93 (MQ=255)

                                                             tCATGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAgaagaa                            <  1:2729364/93‑1 (MQ=255)

                                                               aTGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGc                          >  1:3695983/1‑93 (MQ=255)

                                                                tagATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCa                         <  1:6104880/91‑1 (MQ=255)

                                                                tGGATTAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCa                         >  2:1065209/1‑93 (MQ=255)

                                                                   aTTAACAGAATGAAATTAATCACCCCACAATCCAGAGCATGACTTGCTCGGCCGTTTTGCTACGCCTGTTCTGCGGCAGGAGAATAATCACgc                      >  1:406428/1‑92 (MQ=255)

                                                                    ttAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTc                     <  2:3321619/93‑1 (MQ=255)

                                                                     tAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTcc                    <  2:614324/93‑1 (MQ=255)

                                                                     tAACACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTcc                    <  2:4630293/93‑1 (MQ=255)

                                                                      aaCACAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCg                   >  1:2006575/1‑93 (MQ=255)

                                                                       acacAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCgg                  <  1:5917520/93‑1 (MQ=255)

                                                                       acacAATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGAACGTCCgg                  >  1:2597723/1‑93 (MQ=255)

                                                                        cacaATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGc                 <  2:5841138/93‑1 (MQ=255)

                                                                         acaATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGGAGCACGTCCGGCa                >  1:5623054/1‑93 (MQ=255)

                                                                         acaATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCa                >  2:6555429/1‑93 (MQ=255)

                                                                         acaATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAAAAGAAGCACGTCCGGCa                >  1:1278579/1‑93 (MQ=255)

                                                                          caATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAg               >  2:429443/1‑93 (MQ=255)

                                                                          caATGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAg               >  1:582309/1‑93 (MQ=255)

                                                                            aTGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGaa             >  2:5637251/1‑93 (MQ=255)

                                                                            aTGAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGaa             >  1:3175354/1‑93 (MQ=255)

                                                                              gAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACa           <  1:1811411/93‑1 (MQ=255)

                                                                              gAAATTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACa           >  1:89527/1‑93 (MQ=255)

                                                                               aaaTTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACAt          <  2:2465283/93‑1 (MQ=255)

                                                                               aaaTTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACAt          <  2:745777/93‑1 (MQ=255)

                                                                                 aTTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATga        <  2:515962/93‑1 (MQ=255)

                                                                                 aTTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATga        >  2:2860099/1‑93 (MQ=255)

                                                                                 aTTAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATga        >  1:177790/1‑93 (MQ=255)

                                                                                  ttAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATgac       <  1:3066170/93‑1 (MQ=255)

                                                                                   tAATCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATgacg      >  2:1006613/1‑93 (MQ=255)

                                                                                     aTCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATgacgac    >  1:6484214/1‑93 (MQ=255)

                                                                                     aTCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATgacgac    >  2:1530987/1‑93 (MQ=255)

                                                                                      tCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGaca   >  2:4995534/1‑93 (MQ=255)

                                                                                      tCACCGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACCGCCGGCAGAACATGACGaca   >  1:505631/1‑93 (MQ=255)

ATAGTTTCTGCCTTTATATTTTTATAATTACAACGATAAAAGGCTGTACTTTTTCTTTAGCTCATGGATTAACACAATGAAATTAATCACTGCACCATGCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACAC  >  NC_000913/3225366‑3225545

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