BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,923,857	T→C	Q76R (CAG→CGG)	mioC ←	MioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,923,857	0	T	C	99.5%	1345.0 / ‑4.3	425	Q76R (CAG→CGG)	mioC	MioC
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (217/206); total (217/208)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 432 total aligned reads displayed.

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGG  >  NC_000913/3923767‑3923947

tCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                          <  2:5987863/93‑1 (MQ=255)

tCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                          >  1:3720758/1‑93 (MQ=255)

 cGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGtt                                                                                         >  1:2625798/1‑93 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                       >  1:1678900/1‑93 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                       <  2:6117281/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                       <  2:5359685/93‑1 (MQ=255)

    gTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCt                                                                                      <  1:5093245/93‑1 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                   >  1:4440379/1‑93 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                   >  2:4442111/1‑93 (MQ=255)

         tCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAAt                                                                                 <  1:4745172/93‑1 (MQ=255)

         tCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAAt                                                                                 >  2:454114/1‑93 (MQ=255)

         tCGATAGCCCCACAAAAGGAGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTGCCGTTCCTGCAAt                                                                                 >  1:6993589/1‑93 (MQ=255)

           gATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGc                                                                               <  1:5691234/93‑1 (MQ=255)

           gATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGc                                                                               <  1:738102/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                           <  1:6798952/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                           <  2:7735554/93‑1 (MQ=255)

                ccccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCAt                                                                          <  1:805474/93‑1 (MQ=255)

                ccccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCAt                                                                          <  1:4879423/93‑1 (MQ=255)

                  ccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATag                                                                        <  2:56817/93‑1 (MQ=255)

                    acaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaa                                                                      >  2:5399495/1‑93 (MQ=255)

                     caAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaa                                                                     <  1:3142683/93‑1 (MQ=255)

                     caAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaa                                                                     <  1:1556876/93‑1 (MQ=255)

                      aaaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaag                                                                    <  2:316162/93‑1 (MQ=255)

                      aaaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaag                                                                    >  2:2780467/1‑93 (MQ=255)

                       aaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaagg                                                                   <  1:1437109/93‑1 (MQ=255)

                        aaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaagga                                                                  >  1:3919756/1‑93 (MQ=255)

                         aGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggag                                                                 >  1:1897196/1‑93 (MQ=255)

                          ggTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggaga                                                                >  1:6539715/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               >  1:897098/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               <  2:3247526/93‑1 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               >  2:4783261/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               >  2:4849467/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               <  1:10688/93‑1 (MQ=255)

                           gtgtCAGATTCACGACTTCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                               <  2:1302304/93‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                             <  2:5235760/93‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                             >  1:6880607/1‑93 (MQ=255)

                              tCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaagg                                                            >  2:4994022/1‑93 (MQ=255)

                              tCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaagg                                                            <  1:580575/93‑1 (MQ=255)

                               cATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGt                                                           >  1:7207879/1‑93 (MQ=255)

                                atatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGtt                                                          <  2:7282789/93‑1 (MQ=255)

                                atatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGtt                                                          <  2:6574637/93‑1 (MQ=255)

                                 tatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTg                                                         >  1:6465305/1‑93 (MQ=255)

                                   ttCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTc                                                       >  1:7840693/1‑93 (MQ=255)

                                    tCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTcc                                                      >  1:4260183/1‑93 (MQ=255)

                                    tCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTcc                                                      >  1:1016800/1‑93 (MQ=255)

                                      aCGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCgg                                                    >  2:4244673/1‑93 (MQ=255)

                                      aCGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCgg                                                    <  1:1205530/93‑1 (MQ=255)

                                       cGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGa                                                   <  2:7499004/93‑1 (MQ=255)

                                        gACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGaa                                                  >  1:2575636/1‑93 (MQ=255)

                                         aCTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAat                                                 >  2:2694861/1‑93 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                                >  1:5192429/1‑93 (MQ=255)

                                           tGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAatat                                               <  1:4599593/93‑1 (MQ=255)

                                           tGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAatat                                               >  1:3936443/1‑93 (MQ=255)

                                            gCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtc                                              >  2:2434118/1‑93 (MQ=255)

                                            gCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtc                                              >  2:6500742/1‑93 (MQ=255)

                                            gCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtc                                              <  1:6034305/93‑1 (MQ=255)

                                            gCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtc                                              >  2:2948797/1‑93 (MQ=255)

                                             ccAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAatatat                                             >  1:1934781/1‑91 (MQ=255)

                                             ccAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtct                                             <  2:1642173/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                            <  2:81930/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                            <  2:1754727/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                            >  2:1014383/1‑93 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                           >  2:7691232/1‑93 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                           <  1:401791/93‑1 (MQ=255)

                                                    cGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAc                                      >  2:2162680/1‑93 (MQ=255)

                                                    cGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAc                                      >  1:2259894/1‑93 (MQ=255)

                                                     gATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAcc                                     >  2:6724564/1‑93 (MQ=255)

                                                      aTTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCg                                    >  1:3287989/1‑93 (MQ=255)

                                                      aTTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCg                                    <  1:6342790/93‑1 (MQ=255)

                                                      aTTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCg                                    >  2:1474703/1‑93 (MQ=255)

                                                       ttGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgt                                   <  2:5564487/93‑1 (MQ=255)

                                                           gcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggt                               >  2:1587372/1‑93 (MQ=255)

                                                            ccAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggtg                              >  2:2076544/1‑93 (MQ=255)

                                                             cAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggtgg                             <  2:3605858/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                            <  2:2364261/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                            <  1:1814993/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGGGGa                            >  1:3661801/1‑93 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGGTGTCCGGAATATCTCCGGCACCGTGGGTGGa                            >  2:6978601/1‑93 (MQ=255)

                                                               aaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAg                           <  2:2865023/93‑1 (MQ=255)

                                                                 gCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCt                         <  2:717443/93‑1 (MQ=255)

                                                                   ggACTGCAGAAAGATCGGGCTTCCGTTCCTGCATTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGa                       <  1:7631142/93‑1 (MQ=255)

                                                                    gACTGCAGAAAGATCAGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGAt                      >  1:3832759/1‑93 (MQ=255)

                                                                      cTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATaa                    <  1:3418216/93‑1 (MQ=255)

                                                                       tGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAc                   >  2:62910/1‑93 (MQ=255)

                                                                       tGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAc                   <  1:532181/93‑1 (MQ=255)

                                                                        gCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAcc                  >  1:7729095/1‑93 (MQ=255)

                                                                          aGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccag                >  2:3959116/1‑93 (MQ=255)

                                                                          aGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccag                >  1:3041295/1‑93 (MQ=255)

                                                                            aaaGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagcc              <  2:1973179/93‑1 (MQ=255)

                                                                            aaaGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagcc              >  1:3528626/1‑93 (MQ=255)

                                                                              aGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagccag            >  1:27662/1‑93 (MQ=255)

                                                                               gATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGa           >  1:2869606/1‑93 (MQ=255)

                                                                                   gggCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATccc       >  1:4819500/1‑93 (MQ=255)

                                                                                   gggCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATccc       <  1:6613539/93‑1 (MQ=255)

                                                                                    ggCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCt      <  2:1416053/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg     >  2:3462063/1‑93 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg     <  1:7506447/93‑1 (MQ=255)

                                                                                      cTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa    <  2:4060188/93‑1 (MQ=255)

                                                                                      cTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa    <  2:6827319/93‑1 (MQ=255)

                                                                                      cTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa    <  1:2952885/93‑1 (MQ=255)

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGG  >  NC_000913/3923767‑3923947

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