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.6%	865.1 / ‑4.8	262	M14V (ATG→GTG)	ybdQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (120/141); total (120/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 265 total aligned reads displayed.

TGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640962‑641140

tgACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTc                                                                                        <  2:1017056/93‑1 (MQ=255)

tgACCCGGGTAGTACGTGAAGCAGATGAATAACTCCGGCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTc                                                                                        >  1:1096211/1‑93 (MQ=255)

 gACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCa                                                                                       >  2:4191933/1‑93 (MQ=255)

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

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

      gggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATAc                                                                                  >  2:4435106/1‑93 (MQ=255)

       ggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACa                                                                                 <  1:6194919/93‑1 (MQ=255)

       ggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACa                                                                                 <  2:1763848/93‑1 (MQ=255)

        gtagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACAt                                                                                <  2:3896251/93‑1 (MQ=255)

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

         tagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATc                                                                               <  2:3148631/93‑1 (MQ=255)

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

           gtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCaa                                                                             >  1:3388324/1‑93 (MQ=255)

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

            gaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            <  1:6668891/92‑1 (MQ=255)

              cGTGAAGTAGATGAATACCTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTg                                                                          <  2:5660424/93‑1 (MQ=255)

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

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

               gTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTgg                                                                         <  2:273025/93‑1 (MQ=255)

               gTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTgg                                                                         <  1:4873333/93‑1 (MQ=255)

                tGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGc                                                                        <  2:6268178/93‑1 (MQ=255)

                tGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGc                                                                        >  1:594059/1‑93 (MQ=255)

                tGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGc                                                                        >  2:4691394/1‑93 (MQ=255)

                 gAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCa                                                                       >  1:100668/1‑93 (MQ=255)

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

                  aaGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATATATCAACTGGCAt                                                                      >  2:6142279/1‑93 (MQ=255)

                  aaGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCAt                                                                      >  1:4269149/1‑93 (MQ=255)

                  aaGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCAt                                                                      <  1:4683547/93‑1 (MQ=255)

                   aGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGCCGCTCAATTCCACTTCACATACATCAACTGGCATa                                                                     >  2:2831802/1‑93 (MQ=255)

                    gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    >  2:1952239/1‑93 (MQ=255)

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

                    gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    >  1:2839569/1‑93 (MQ=255)

                    gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    >  1:6364936/1‑93 (MQ=255)

                    ctAGATGAATAACTCCGTCATCCTGGGCGAGGAATGCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    <  1:5363672/92‑1 (MQ=255)

                       gATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGa                                                                 <  1:5217479/93‑1 (MQ=255)

                        aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                <  2:1216408/93‑1 (MQ=255)

                        aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                <  2:4779515/93‑1 (MQ=255)

                        aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                <  2:5255626/93‑1 (MQ=255)

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

                             tAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCt                                                           <  2:6168437/93‑1 (MQ=255)

                             tAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCt                                                           >  2:6478480/1‑93 (MQ=255)

                               aCTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTta                                                         >  2:6635944/1‑93 (MQ=255)

                                 tCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtatc                                                       <  1:6436541/93‑2 (MQ=255)

                                  ccGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTACACTTCAAATACATCAACTGGCATAATGATTGTCTTATAc                                                      >  1:1568657/1‑93 (MQ=255)

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

                                    gTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACAt                                                    >  1:1315895/1‑93 (MQ=255)

                                    gTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACAt                                                    >  2:4547974/1‑93 (MQ=255)

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

                                      cATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATaa                                                  <  2:2507471/93‑1 (MQ=255)

                                       aTCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAc                                                 >  2:1561306/1‑93 (MQ=255)

                                            gggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCttt                                            >  1:5104752/1‑93 (MQ=255)

                                            gggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCttt                                            >  1:6563499/1‑93 (MQ=255)

                                            gggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCttt                                            <  2:3399243/93‑1 (MQ=255)

                                             ggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtc                                           >  2:2550175/1‑93 (MQ=255)

                                              gCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtct                                          <  1:3434818/93‑1 (MQ=255)

                                                 aGGAATTCAGCGTGGTGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTccc                                       <  2:732483/93‑1 (MQ=255)

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

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

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

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

                                                   gAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTg                                     >  2:3484907/1‑93 (MQ=255)

                                                   gAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTg                                     <  1:681106/93‑1 (MQ=255)

                                                     aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   >  2:89961/1‑93 (MQ=255)

                                                     aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   >  1:2100073/1‑93 (MQ=255)

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

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

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

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

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

                                                           cGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCAt                             >  2:2653256/1‑93 (MQ=255)

                                                           cGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCAt                             >  1:3645881/1‑93 (MQ=255)

                                                            gTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATg                            <  2:5764316/93‑1 (MQ=255)

                                                            gTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATg                            >  1:390852/1‑93 (MQ=255)

                                                              ggCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGaa                          <  2:3995919/93‑1 (MQ=255)

                                                              ggCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGaa                          >  1:2639423/1‑93 (MQ=255)

                                                                 gAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaa                       >  1:5058200/1‑93 (MQ=255)

                                                                  aaCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      <  2:3016090/93‑1 (MQ=255)

                                                                  aaCAGCTTTGTCGCTCAATTCCACTCCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      <  1:2998450/93‑1 (MQ=255)

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

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

                                                                    cAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATc                    >  2:544210/1‑93 (MQ=255)

                                                                     aGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCa                   <  2:3328122/93‑1 (MQ=255)

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

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

                                                                       cTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAtt                 >  2:4035186/1‑93 (MQ=255)

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

                                                                        tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                >  2:6535943/1‑93 (MQ=255)

                                                                        cttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  1:4491838/92‑1 (MQ=255)

                                                                         ttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCg               <  2:4807941/93‑1 (MQ=255)

                                                                             cGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCAt           >  1:5826794/1‑93 (MQ=255)

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

                                                                               cTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGa         <  1:1226737/93‑1 (MQ=255)

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

                                                                                 cAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAtt       >  2:4108920/1‑93 (MQ=255)

                                                                                 cAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAtt       <  1:5517359/93‑1 (MQ=255)

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

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

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

                                                                                     tCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATCCGCCATGATTataa   <  1:4064851/93‑1 (MQ=255)

                                                                                      ccACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTataat  >  2:2155151/1‑93 (MQ=255)

TGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640962‑641140

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