BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	199,879	A→G	N651S (AAC→AGC)	yaeT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	199,879	0	A	G	99.4%	554.7 / ‑4.7	176	N651S (AAC→AGC)	yaeT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (84/91); total (85/91)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

Only 100 of 178 total aligned reads displayed.

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCt                                                                                         <  2:192602/93‑1 (MQ=255)

tcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCt                                                                                         >  2:1110164/1‑93 (MQ=255)

  gatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTc                                                                                       <  1:2924135/93‑1 (MQ=255)

  gatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTc                                                                                       <  1:941241/93‑1 (MQ=255)

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      <  1:382478/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGTTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  2:1642512/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  1:2621778/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  2:2355191/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  1:425376/93‑1 (MQ=255)

    tGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTa                                                                                     <  2:1560946/93‑1 (MQ=255)

      aCGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATg                                                                                   >  2:3254287/1‑93 (MQ=255)

      aCGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATg                                                                                   >  1:1526602/1‑93 (MQ=255)

        gATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGcc                                                                                 >  2:2678403/1‑93 (MQ=255)

         aTCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCg                                                                                >  2:4082012/1‑93 (MQ=255)

         aTCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCg                                                                                >  2:2336889/1‑93 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               >  1:109019/1‑93 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               <  2:3144495/93‑1 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               >  2:2847010/1‑93 (MQ=255)

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCgg                                                                               >  1:102811/1‑93 (MQ=255)

            acaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtg                                                                             >  1:2270371/1‑93 (MQ=255)

             caAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtgg                                                                            >  1:3215919/1‑93 (MQ=255)

              aaaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtggt                                                                           <  2:1233293/93‑1 (MQ=255)

               aaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAGGAGATGCCGTTCTACGAGAGCTTCTGTGCCGGGGGtt                                                                          >  2:3182624/1‑93 (MQ=255)

               aaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGtt                                                                          >  1:3447050/1‑93 (MQ=255)

                aTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTc                                                                         >  2:3596463/1‑93 (MQ=255)

                 tGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTcc                                                                        >  1:3897866/1‑93 (MQ=255)

                  gggTTGTTCTGGGGCGTTCCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCa                                                                       <  1:3668820/93‑1 (MQ=255)

                    gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGc                                                                     <  1:4129472/93‑1 (MQ=255)

                    gttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGc                                                                     <  2:1697106/93‑1 (MQ=255)

                        ttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcg                                                                 <  1:774637/93‑1 (MQ=255)

                         tCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgt                                                                >  1:2495911/1‑93 (MQ=255)

                          cTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtg                                                               <  1:1713622/93‑1 (MQ=255)

                           tGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgc                                                              <  2:4617920/93‑1 (MQ=255)

                            ggggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcg                                                             <  2:2488837/93‑1 (MQ=255)

                             gggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcgt                                                            <  2:3634902/93‑1 (MQ=255)

                             gggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcgt                                                            <  1:4198855/93‑1 (MQ=255)

                              ggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcgtg                                                           >  1:1649726/1‑93 (MQ=255)

                                cGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGc                                                         <  1:388325/93‑1 (MQ=255)

                                 gTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCt                                                        >  1:1505617/1‑93 (MQ=255)

                                  tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCtt                                                       >  1:3051771/1‑93 (MQ=255)

                                    cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTcc                                                     <  2:4560228/93‑1 (MQ=255)

                                     ccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCa                                                    <  2:1634184/93‑1 (MQ=255)

                                      cGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAg                                                   >  1:1615040/1‑93 (MQ=255)

                                       gCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGt                                                  >  1:2373805/1‑93 (MQ=255)

                                       gCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGt                                                  <  2:167214/93‑1 (MQ=255)

                                         tGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTcc                                                <  1:2599356/93‑1 (MQ=255)

                                           gggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCaa                                              >  2:4434412/1‑93 (MQ=255)

                                           gggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCaa                                              <  1:2797411/93‑1 (MQ=255)

                                           gggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCaa                                              >  1:3973681/1‑93 (MQ=255)

                                             gTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATa                                            >  1:299023/1‑93 (MQ=255)

                                             gTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATa                                            <  2:329928/93‑1 (MQ=255)

                                                aTGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCa                                         <  1:4667908/93‑1 (MQ=255)

                                                 tGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCAt                                        >  1:175451/1‑93 (MQ=255)

                                                 tGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCAt                                        <  1:34983/93‑1 (MQ=255)

                                                 tGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCAt                                        <  2:3893396/93‑1 (MQ=255)

                                                 tGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCAt                                        <  1:3392714/93‑1 (MQ=255)

                                                  ggTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCAtt                                       <  2:190695/93‑1 (MQ=255)

                                                   gTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTg                                      <  1:1193869/93‑1 (MQ=255)

                                                   gTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTg                                      >  1:1185434/1‑93 (MQ=255)

                                                    tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTgg                                     >  1:3427747/1‑93 (MQ=255)

                                                    tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTgg                                     <  1:2333320/93‑1 (MQ=255)

                                                    tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTgg                                     >  2:2116558/1‑93 (MQ=255)

                                                    tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTgg                                     >  1:1330699/1‑93 (MQ=255)

                                                     gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGt                                    >  1:117319/1‑93 (MQ=255)

                                                     gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGt                                    >  1:1246093/1‑93 (MQ=255)

                                                         gTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGa                                >  1:702786/1‑93 (MQ=255)

                                                         gTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGa                                >  2:4333371/1‑93 (MQ=255)

                                                          tttAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGaa                               >  1:115723/1‑93 (MQ=255)

                                                          tttAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGaa                               >  1:1436219/1‑93 (MQ=255)

                                                           ttAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGaaa                              <  2:4389570/93‑1 (MQ=255)

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAACCa                           <  2:662996/93‑1 (MQ=255)

                                                               gcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAg                          <  1:655101/93‑1 (MQ=255)

                                                                cggcAAATAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGt                         <  1:2265873/93‑1 (MQ=255)

                                                                cggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGt                         <  2:496937/93‑1 (MQ=255)

                                                                cggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGt                         <  1:3219777/93‑1 (MQ=255)

                                                                   cAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTa                      >  2:884014/1‑93 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATCCCATTGGTCCGAAAGCAGTTTAc                     <  2:1930547/93‑1 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     >  1:4426079/1‑93 (MQ=255)

                                                                    aaaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTAc                     >  1:59862/1‑93 (MQ=255)

                                                                     aaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACt                    <  2:691830/93‑1 (MQ=255)

                                                                     aaGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACt                    >  2:1299132/1‑93 (MQ=255)

                                                                          aTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCg               >  1:3525900/1‑93 (MQ=255)

                                                                          aTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATCGGTCCGAAAGCAGTTTACTTCCCg               <  2:2304156/93‑1 (MQ=255)

                                                                           tGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGc              <  2:2752646/93‑1 (MQ=255)

                                                                            gCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCa             <  2:3193908/93‑1 (MQ=255)

                                                                             ccGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCAt            <  2:538300/93‑1 (MQ=255)

                                                                              cGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATc           <  1:4601927/93‑1 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCGTCa          >  2:3418355/1‑93 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          >  1:4353533/1‑93 (MQ=255)

                                                                               gTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCa          >  1:2450268/1‑93 (MQ=255)

                                                                                ttCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg         <  1:4469739/93‑1 (MQ=255)

                                                                                ttCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg         >  1:2288674/1‑93 (MQ=255)

                                                                                ttCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg         >  2:1423962/1‑93 (MQ=255)

                                                                                  cTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc       <  1:1993129/93‑1 (MQ=255)

                                                                                   tACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGcc      >  2:2995291/1‑93 (MQ=255)

                                                                                    aCGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCa     <  2:2536959/93‑1 (MQ=255)

                                                                                    aCGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCa     >  1:203277/1‑93 (MQ=255)

                                                                                     cGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAg    >  1:4648566/1‑93 (MQ=255)

                                                                                     cGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAg    <  2:1298285/93‑1 (MQ=255)

                                                                                      gagaGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGt   >  1:2466518/1‑93 (MQ=255)

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

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