BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	641,051	T→C	M14V (ATG→GTG)	ybdQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	641,051	0	T	C	99.4%	1041.5 / ‑4.0	326	M14V (ATG→GTG)	ybdQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (149/175); total (149/177)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 330 total aligned reads displayed.

CTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640961‑641140

ctgACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCActt                                                                                         >  1:4201007/1‑93 (MQ=255)

   aCCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCaa                                                                                      >  1:841620/1‑93 (MQ=255)

    cccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCaaa                                                                                     <  2:6279974/93‑1 (MQ=255)

    cccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCaaa                                                                                     <  2:3062654/93‑1 (MQ=255)

     ccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAAt                                                                                    >  1:2373263/1‑93 (MQ=255)

      cGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATa                                                                                   <  2:5084957/93‑1 (MQ=255)

      cGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATa                                                                                   <  2:4139987/93‑1 (MQ=255)

       gggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATAc                                                                                  <  1:5748775/93‑1 (MQ=255)

        ggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACa                                                                                 >  2:822710/1‑93 (MQ=255)

        ggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACa                                                                                 >  1:6558291/1‑93 (MQ=255)

         gtagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACAt                                                                                >  2:2469025/1‑93 (MQ=255)

          tagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATc                                                                               >  2:6776786/1‑93 (MQ=255)

            gtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCaa                                                                             <  1:1092691/93‑1 (MQ=255)

             taCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            <  2:7431894/93‑1 (MQ=255)

             taCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            <  2:6522149/93‑1 (MQ=255)

             taCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            <  1:5479003/93‑1 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           >  1:3673544/1‑93 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           >  1:4841971/1‑93 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           <  1:2935355/93‑1 (MQ=255)

               cGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTg                                                                          <  2:5504058/93‑1 (MQ=255)

                gTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTgg                                                                         >  1:7059724/1‑93 (MQ=255)

                 ggAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATAAATCAACTGGc                                                                        <  2:1346829/92‑1 (MQ=255)

                  gAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCa                                                                       <  1:7897650/93‑1 (MQ=255)

                  gAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCa                                                                       <  1:755452/93‑1 (MQ=255)

                    aGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATa                                                                     <  1:7239422/93‑1 (MQ=255)

                    aGTAGATGAATAACTCCGTCATCATGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATa                                                                     <  2:4336120/93‑1 (MQ=255)

                     gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    <  1:1010387/93‑1 (MQ=255)

                     gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    <  2:1592502/93‑1 (MQ=255)

                        gATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGa                                                                 >  2:7373933/1‑93 (MQ=255)

                        gATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGa                                                                 >  2:3872911/1‑93 (MQ=255)

                         aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                >  1:1978436/1‑93 (MQ=255)

                         aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                >  1:1519878/1‑93 (MQ=255)

                          tGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAtt                                                               >  2:3370692/1‑93 (MQ=255)

                          tGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCCAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAtt                                                               <  1:4529537/93‑1 (MQ=255)

                            aaTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTgc                                                             >  2:3096001/1‑92 (MQ=255)

                            aaTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGt                                                             >  1:4416115/1‑93 (MQ=255)

                                aCTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTta                                                         <  1:5933301/93‑1 (MQ=255)

                                aCTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTta                                                         <  1:751429/93‑1 (MQ=255)

                                  tCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtata                                                       >  1:4009907/1‑93 (MQ=255)

                                  tCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtata                                                       <  1:5216732/93‑1 (MQ=255)

                                  tCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtata                                                       >  2:1936955/1‑93 (MQ=255)

                                    cggcATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     <  2:2204195/90‑1 (MQ=255)

                                    cGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     <  2:5683542/93‑1 (MQ=255)

                                    cGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     >  1:1866562/1‑93 (MQ=255)

                                    cGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     <  2:5311126/93‑1 (MQ=255)

                                      tCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATa                                                   <  1:461651/93‑1 (MQ=255)

                                      tCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATa                                                   >  1:3120991/1‑93 (MQ=255)

                                        aTCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAc                                                 <  2:6074825/93‑1 (MQ=255)

                                         tCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAcc                                                >  1:4066282/1‑93 (MQ=255)

                                          ccTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAccc                                               >  1:25765/1‑93 (MQ=255)

                                          ccTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAccc                                               <  1:2901657/93‑1 (MQ=255)

                                           cTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCt                                              <  1:4736569/93‑1 (MQ=255)

                                              ggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtc                                           <  2:4391269/93‑1 (MQ=255)

                                                ggAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtctc                                         <  1:6115917/92‑1 (MQ=255)

                                                cGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtctc                                         >  2:1741197/1‑93 (MQ=255)

                                                 gAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTcc                                        >  1:50396/1‑93 (MQ=255)

                                                 gAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTcc                                        >  1:5961441/1‑93 (MQ=255)

                                                  aGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTccc                                       >  1:5790908/1‑93 (MQ=255)

                                                  aGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTccc                                       <  2:2873406/93‑1 (MQ=255)

                                                   ggAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCt                                      >  1:1300878/1‑93 (MQ=255)

                                                   ggAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCt                                      <  2:3543956/93‑1 (MQ=255)

                                                     aaTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGt                                    >  2:6166194/1‑93 (MQ=255)

                                                     aaTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGt                                    >  2:270266/1‑93 (MQ=255)

                                                     aaTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGt                                    >  2:3678556/1‑93 (MQ=255)

                                                     aaTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGt                                    >  2:2317349/1‑93 (MQ=255)

                                                      aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   <  1:7695978/93‑1 (MQ=255)

                                                       ttCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTa                                  >  2:5994738/1‑93 (MQ=255)

                                                       ttCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTa                                  >  2:1083308/1‑93 (MQ=255)

                                                        tCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTaa                                 >  1:3289236/1‑93 (MQ=255)

                                                        tCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTaa                                 <  2:4489468/93‑1 (MQ=255)

                                                        tCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTaa                                 <  2:6939672/93‑1 (MQ=255)

                                                         cAGCGTGTCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAAt                                <  2:1416945/93‑1 (MQ=255)

                                                           gCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCa                              <  2:2962235/93‑1 (MQ=255)

                                                           gCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCa                              <  2:4025707/93‑1 (MQ=255)

                                                           gCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCa                              >  2:6878073/1‑93 (MQ=255)

                                                              tGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGa                           <  2:3783107/93‑1 (MQ=255)

                                                              tGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGa                           >  2:4740037/1‑93 (MQ=255)

                                                                 ggAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACa                        <  1:2266360/92‑1 (MQ=255)

                                                                   aaCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      >  1:6216549/1‑93 (MQ=255)

                                                                    aCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAAt                     >  2:7181168/1‑93 (MQ=255)

                                                                    aCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAAt                     >  1:6799042/1‑93 (MQ=255)

                                                                     cAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATc                    <  2:5610917/93‑1 (MQ=255)

                                                                      aGCTTTGTCGCTCAATGCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCa                   <  2:2038965/93‑1 (MQ=255)

                                                                       gCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAt                  >  2:3264058/1‑93 (MQ=255)

                                                                       gCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAt                  >  1:7787653/1‑93 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  1:5303396/93‑1 (MQ=255)

                                                                          ttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCg               >  2:536906/1‑93 (MQ=255)

                                                                            gTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGca             >  1:1193943/1‑92 (MQ=255)

                                                                             tCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCa            >  2:6133575/1‑93 (MQ=255)

                                                                             tCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCa            <  1:3880112/93‑1 (MQ=255)

                                                                              cGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCAt           <  1:7290515/93‑1 (MQ=255)

                                                                               gCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATg          >  1:5053793/1‑93 (MQ=255)

                                                                               gCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATg          >  2:7850188/1‑93 (MQ=255)

                                                                                 tCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAt        >  2:3426606/1‑93 (MQ=255)

                                                                                 tCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAt        >  1:6925437/1‑93 (MQ=255)

                                                                                 tCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAt        >  2:7721923/1‑93 (MQ=255)

                                                                                  cAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAtt       <  1:5630507/93‑1 (MQ=255)

                                                                                   aaTTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATta      <  1:964114/93‑1 (MQ=255)

                                                                                     ttCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATtata    <  2:5770680/93‑1 (MQ=255)

                                                                                      tCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTataa   <  2:6384054/93‑1 (MQ=255)

CTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640961‑641140

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