BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,199,347	T→C	R21G (AGA→GGA)	yjaI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,199,347	0	T	C	86.2%	850.7 / 122.7	354	R21G (AGA→GGA)	yjaI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (23/26); new base C (146/159); total (169/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-02

Only 100 of 355 total aligned reads displayed.

TCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTACAGCGCAGGG  >  NC_000913/4199257‑4199439

tCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCgg                                                                                            >  2:4260245/1‑93 (MQ=255)

  cGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCa                                                                                          >  1:2007238/1‑93 (MQ=255)

   gTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCAt                                                                                         <  1:6027942/93‑1 (MQ=255)

   gTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCAt                                                                                         <  1:2567224/93‑1 (MQ=255)

     ttCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAAGCCGGCATGa                                                                                       <  1:4296356/93‑1 (MQ=255)

      tCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGAt                                                                                      <  2:7303335/93‑1 (MQ=255)

      tCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGAt                                                                                      <  2:4585734/93‑1 (MQ=255)

      tCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGAt                                                                                      <  2:432347/93‑1 (MQ=255)

         tGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATaaa                                                                                   >  1:4959377/1‑93 (MQ=255)

            ttAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGa                                                                                >  1:203550/1‑93 (MQ=255)

              aaTCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGt                                                                              >  1:4244102/1‑93 (MQ=255)

              aaTCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGt                                                                              >  1:5279794/1‑93 (MQ=255)

                tcctccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTaa                                                                            <  2:2505264/93‑1 (MQ=255)

                tcctccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaa                                                                            >  2:2297751/1‑93 (MQ=255)

                tcctccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaa                                                                            <  2:3372708/93‑1 (MQ=255)

                  ctccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaa                                                                          >  2:2762675/1‑93 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         <  1:361030/93‑1 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         <  1:6534769/93‑1 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         >  1:5514288/1‑93 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         <  2:1900972/93‑1 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         >  1:2218579/1‑93 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         <  2:5587070/93‑1 (MQ=255)

                   tccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTaaaaa                                                                         <  2:1049846/93‑1 (MQ=255)

                    ccAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAAt                                                                        <  1:3002729/93‑1 (MQ=255)

                     cAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATg                                                                       <  2:5474781/93‑1 (MQ=255)

                     cAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATg                                                                       >  2:387574/1‑93 (MQ=255)

                      aGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGa                                                                      <  1:5741619/93‑1 (MQ=255)

                       gTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGAc                                                                     <  2:7505528/93‑1 (MQ=255)

                         ggTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTc                                                                   >  2:5989896/1‑93 (MQ=255)

                          gTTGTCTTACTTCGGGTATTGCATCTTCTGTGCCAACGATGAAACGTTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCg                                                                  >  1:3299748/1‑93 (MQ=255)

                           ttGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGc                                                                 >  1:4126716/1‑93 (MQ=255)

                           ttGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGc                                                                 <  1:3915604/93‑1 (MQ=255)

                            tGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGcc                                                                <  1:2749839/93‑1 (MQ=255)

                             gTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCct                                                               <  2:6401974/93‑1 (MQ=255)

                             gTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCct                                                               <  2:2188500/93‑1 (MQ=255)

                             gTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCct                                                               <  1:2791547/93‑1 (MQ=255)

                               cTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgc                                                             <  2:7215033/93‑1 (MQ=255)

                               cTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgc                                                             <  2:6116491/93‑1 (MQ=255)

                                ttACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgct                                                            >  2:2712232/1‑93 (MQ=255)

                                ttACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgct                                                            >  2:1766986/1‑93 (MQ=255)

                                 tACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCctgctg                                                           <  2:6765244/93‑1 (MQ=255)

                                 tACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgctg                                                           >  1:5794339/1‑93 (MQ=255)

                                 tACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCctgctg                                                           <  2:7166476/93‑1 (MQ=255)

                                  aCTTCGGGTATTGCATCTTCCGTGCCAACGATTAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCtgctgc                                                          <  1:69198/93‑1 (MQ=255)

                                    ttCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCgg                                                        <  1:5168205/93‑1 (MQ=255)

                                     tCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCggg                                                       >  1:5940899/1‑93 (MQ=255)

                                     tCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCggg                                                       <  2:3471337/93‑1 (MQ=255)

                                      cGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGt                                                      <  2:6282550/93‑1 (MQ=255)

                                       gggTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTa                                                     >  1:1057228/1‑93 (MQ=255)

                                       gggTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTa                                                     <  1:2336257/93‑1 (MQ=255)

                                         gTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGc                                                   <  2:3749681/93‑1 (MQ=255)

                                         gTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGc                                                   >  2:5232567/1‑93 (MQ=255)

                                          tATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCg                                                  >  2:6152541/1‑93 (MQ=255)

                                          tATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCg                                                  >  1:7313109/1‑93 (MQ=255)

                                          tATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCg                                                  >  2:4847466/1‑93 (MQ=255)

                                          tATTGCATCTTCCCTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCg                                                  >  1:4522801/1‑93 (MQ=255)

                                           aTTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGa                                                 <  1:1825324/93‑1 (MQ=255)

                                            ttGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAg                                                >  1:794555/1‑93 (MQ=255)

                                            ttGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAg                                                >  1:4586785/1‑93 (MQ=255)

                                             tGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGt                                               >  2:3671054/1‑93 (MQ=255)

                                               cATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCa                                             >  2:5187570/1‑93 (MQ=255)

                                                aTCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTGGCGAGTCAt                                            <  1:2642943/93‑1 (MQ=255)

                                                aTCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAt                                            >  1:3205480/1‑93 (MQ=255)

                                                 tCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAtt                                           <  2:2631799/93‑1 (MQ=255)

                                                 tCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAtt                                           >  2:1405646/1‑93 (MQ=255)

                                                   ttCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAtttt                                         <  2:5435218/93‑1 (MQ=255)

                                                   ttCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAtttt                                         >  2:2264317/1‑93 (MQ=255)

                                                   ttCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCAtttt                                         >  1:1901516/1‑93 (MQ=255)

                                                      cGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTAc                                      <  1:4770121/93‑1 (MQ=255)

                                                      cGTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTAc                                      >  2:4028911/1‑93 (MQ=255)

                                                       gTGCCAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACt                                     <  1:630709/93‑1 (MQ=255)

                                                        tGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTc                                    >  1:2272483/1‑93 (MQ=255)

                                                           cAACGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCAtt                                 <  1:7189980/93‑1 (MQ=255)

                                                             aCGATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGa                               >  1:3175496/1‑93 (MQ=255)

                                                               gATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGaaa                             <  2:851328/93‑1 (MQ=255)

                                                               gATGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGaaa                             >  1:5862348/1‑93 (MQ=255)

                                                                aTGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAAc                            >  1:5594988/1‑93 (MQ=255)

                                                                 tGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCCTTGAAACt                           >  1:1433795/1‑93 (MQ=255)

                                                                 tGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACt                           >  2:4522631/1‑93 (MQ=255)

                                                                 tGAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACt                           <  1:6067605/93‑1 (MQ=255)

                                                                  gAAACGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACtt                          >  2:5700687/1‑93 (MQ=255)

                                                                    aaCGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTaa                        <  1:352139/93‑1 (MQ=255)

                                                                    aaCGCTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTaa                        <  1:95697/93‑1 (MQ=255)

                                                                        cTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCt                    <  2:3167006/93‑1 (MQ=255)

                                                                        cTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCt                    <  1:3704995/93‑1 (MQ=255)

                                                                        cTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCt                    >  2:652692/1‑93 (MQ=255)

                                                                        cTGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCt                    >  2:328659/1‑93 (MQ=255)

                                                                         tGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCtt                   <  1:4352900/93‑1 (MQ=255)

                                                                         tGATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCtt                   >  2:1585316/1‑93 (MQ=255)

                                                                         tGATATGACGGGTAATCCGGCATGATAAACGAATAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCtt                   >  1:6768280/1‑93 (MQ=255)

                                                                          gATATGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCttt                  <  1:6580264/93‑1 (MQ=255)

                                                                           atatGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTtg                 >  2:591279/1‑93 (MQ=255)

                                                                           atatGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTtg                 <  2:1664777/93‑1 (MQ=255)

                                                                            tatGACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTtgt                <  1:1703598/93‑1 (MQ=255)

                                                                             atGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTtgtg               >  2:2996027/1‑93 (MQ=255)

                                                                               gACGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGtt             <  1:3033115/93‑1 (MQ=255)

                                                                                aCGGGTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTa            >  1:6413252/1‑93 (MQ=255)

                                                                                   ggTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTACAg         <  2:70804/93‑1 (MQ=255)

                                                                                    gTAATCCGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTACAgc        >  2:1333334/1‑93 (MQ=255)

                                                                                       aTCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTACAGCGCa     >  1:119999/1‑93 (MQ=255)

TCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGTATTGCATCTTCCGTGCCAACGATGAAACGCTGATATGACGGGTAATCTGGCATGATAAACGAGTAAAAATGACTCGCCTGCTGCGGGTAGCGAGTCATTTTTACTCATTGAAACTTAAGCCTTTGTGTTACAGCGCAGGG  >  NC_000913/4199257‑4199439

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