BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,958	G→A	G439S (GGC→AGC)	yjeF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,958	0	G	A	99.5%	1330.0 / ‑6.1	364	G439S (GGC→AGC)	yjeF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (175/187); total (175/189)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 390 total aligned reads displayed.

GCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392868‑4393047

gCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCggcagc                                                                                         >  2:5915233/1‑93 (MQ=255)

 cGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGgcagca                                                                                        <  2:1695740/93‑1 (MQ=255)

 cGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGgcagca                                                                                        >  2:4017898/1‑93 (MQ=255)

  ggTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCAt                                                                                       <  1:6221166/93‑1 (MQ=255)

   gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      >  1:632281/1‑93 (MQ=255)

   gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      >  2:6288725/1‑93 (MQ=255)

   gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      >  2:1600962/1‑93 (MQ=255)

   gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      <  2:1937977/93‑1 (MQ=255)

    tGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATgg                                                                                     >  2:5017756/1‑93 (MQ=255)

     gCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATggg                                                                                    <  1:7374950/93‑1 (MQ=255)

     gCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATggg                                                                                    <  2:1132358/93‑1 (MQ=255)

       tGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCg                                                                                  >  1:7111311/1‑93 (MQ=255)

        gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGTATGGCGAGCGGCAGCATGGGCGa                                                                                 <  2:3316036/93‑1 (MQ=255)

        gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 >  1:5955744/1‑93 (MQ=255)

        gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 <  1:3764812/93‑1 (MQ=255)

        gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 <  2:264277/93‑1 (MQ=255)

         aaaGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAt                                                                                <  2:482740/93‑1 (MQ=255)

         aaaGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAt                                                                                <  1:4772336/93‑1 (MQ=255)

            ggTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAtgtg                                                                             >  2:6400354/1‑93 (MQ=255)

              tGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGct                                                                           <  2:4128331/93‑1 (MQ=255)

               gCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTATGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctc                                                                          <  1:4086358/93‑1 (MQ=255)

               gCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctc                                                                          <  1:489496/93‑1 (MQ=255)

               gCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctc                                                                          <  1:2853733/93‑1 (MQ=255)

                   gAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTg                                                                      <  2:6856290/93‑1 (MQ=255)

                    aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTgg                                                                     <  2:6626912/93‑1 (MQ=255)

                    aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTgg                                                                     <  2:7224797/93‑1 (MQ=255)

                     aCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGt                                                                    <  2:212153/93‑1 (MQ=255)

                       cGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtat                                                                  >  2:3605406/1‑93 (MQ=255)

                       cGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtat                                                                  >  1:7447494/1‑93 (MQ=255)

                          ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtattat                                                               <  2:291752/93‑1 (MQ=255)

                          ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtattat                                                               <  2:812670/93‑1 (MQ=255)

                           gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTattatt                                                              >  1:1190636/1‑93 (MQ=255)

                           gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTattatt                                                              <  2:5249098/93‑1 (MQ=255)

                            tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTg                                                             <  2:7500575/93‑1 (MQ=255)

                            tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTg                                                             <  1:2087974/93‑1 (MQ=255)

                            tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTg                                                             >  2:3847582/1‑93 (MQ=255)

                            gcgccgcCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTg                                                             <  1:676879/92‑1 (MQ=255)

                             cgccgcCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTgg                                                            >  1:862083/1‑93 (MQ=255)

                               ccgccCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGgcg                                                          <  1:2623276/93‑1 (MQ=255)

                                cgccCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGTATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGgcgc                                                         <  2:6551573/93‑1 (MQ=255)

                                 gccCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCa                                                        <  2:5171439/93‑1 (MQ=255)

                                  cccATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAt                                                       <  1:3602343/93‑1 (MQ=255)

                                  cccATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAt                                                       >  1:6999974/1‑93 (MQ=255)

                                     aTCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    <  1:1807854/93‑1 (MQ=255)

                                     aTCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    <  2:7657703/93‑1 (MQ=255)

                                     aTCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    >  1:954447/1‑93 (MQ=255)

                                      tCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCt                                                   >  1:1504160/1‑93 (MQ=255)

                                      tCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCt                                                   >  1:4306160/1‑93 (MQ=255)

                                      tCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCt                                                   >  1:2451625/1‑93 (MQ=255)

                                      tCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTCGCGCATTGCt                                                   >  1:2230170/1‑93 (MQ=255)

                                        cagACGCTTTAGGCCTCATTGA‑GCAGCGAAATGCAGGCATGGCGAGCGGCAGCATGGCCGATATGCCCTCTGGTATTATAGGCGCATTGCTTg                                                 <  2:2766298/91‑1 (MQ=255)

                                        cTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTg                                                 <  2:6105679/93‑1 (MQ=255)

                                         tGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTgg                                                >  1:3946893/1‑93 (MQ=255)

                                         tGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTgg                                                <  1:6889446/93‑1 (MQ=255)

                                          gACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTggg                                               >  1:3699851/1‑93 (MQ=255)

                                           aCGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGc                                              >  2:4704259/1‑93 (MQ=255)

                                            cGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCa                                             >  1:339212/1‑93 (MQ=255)

                                             gCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaa                                            >  1:3210201/1‑93 (MQ=255)

                                             gCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaa                                            >  2:5706380/1‑93 (MQ=255)

                                              cTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaa                                           >  1:2138119/1‑93 (MQ=255)

                                               tttAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaa                                          >  2:7121652/1‑93 (MQ=255)

                                               tctAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaa                                          >  2:6392125/3‑93 (MQ=255)

                                                ttAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaaa                                         >  1:720686/1‑93 (MQ=255)

                                                ttAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaaa                                         >  1:3570514/1‑93 (MQ=255)

                                                 tAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAAc                                        <  2:196892/93‑1 (MQ=255)

                                                 tAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAAc                                        <  1:6667320/93‑1 (MQ=255)

                                                  aGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACt                                       >  2:4468417/1‑93 (MQ=255)

                                                  aGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACt                                       <  1:277808/93‑1 (MQ=255)

                                                     cATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctctc                                                                        <  2:3723112/57‑1 (MQ=255)

                                                        tattGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGcc                                 >  2:1250013/1‑93 (MQ=255)

                                                        tattGATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGcc                                 >  1:1997865/1‑93 (MQ=255)

                                                            gATGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTat                             <  2:2487845/93‑1 (MQ=255)

                                                             aTGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatg                            <  2:230792/93‑1 (MQ=255)

                                                              tGCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatga                           <  2:6227934/93‑1 (MQ=255)

                                                               gCCG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatgat                          >  2:1869815/1‑93 (MQ=255)

                                                                 cG‑GAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGc                        >  1:1658666/1‑93 (MQ=255)

                                                                    gAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAg                      >  2:3050485/1‑93 (MQ=255)

                                                                        tGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtg                  <  1:5125027/93‑1 (MQ=255)

                                                                        tGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtg                  >  1:7293005/1‑93 (MQ=255)

                                                                           aGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc               >  1:3667290/1‑93 (MQ=255)

                                                                            ggCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa              >  1:4593140/1‑93 (MQ=255)

                                                                            ggCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa              >  1:3284342/1‑93 (MQ=255)

                                                                             gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             >  1:7656393/1‑93 (MQ=255)

                                                                             gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             >  1:3159086/1‑93 (MQ=255)

                                                                             gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             <  1:7785253/93‑1 (MQ=255)

                                                                              cATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAgg            >  2:4332267/1‑93 (MQ=255)

                                                                              cATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAgg            <  1:2536593/93‑1 (MQ=255)

                                                                               aTGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGc           <  1:443769/93‑1 (MQ=255)

                                                                               aTGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGc           >  2:3437495/1‑93 (MQ=255)

                                                                                tGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCt          >  1:1297098/1‑93 (MQ=255)

                                                                                  gCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt        <  1:1001241/93‑1 (MQ=255)

                                                                                  gCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt        <  1:5259194/93‑1 (MQ=255)

                                                                                  gCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt        <  1:4937488/93‑1 (MQ=255)

                                                                                   cGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgtg       <  2:1330519/93‑1 (MQ=255)

                                                                                   cGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgtg       <  1:545750/93‑1 (MQ=255)

                                                                                    gAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTgtgt      <  1:5797844/93‑1 (MQ=255)

                                                                                      gcggcaGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcg    <  2:4254244/93‑1 (MQ=255)

                                                                                      gcggcaGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcg    >  1:6008665/1‑93 (MQ=255)

                                                                                       cggcagCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc   <  1:5616583/93‑1 (MQ=255)

                                                                                        ggcagcATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgcg  >  2:5616240/1‑93 (MQ=255)

GCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCG‑GAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392868‑4393047

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