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	98.6%	453.5 / ‑3.5	143	N651S (AAC→AGC)	yaeT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (67/74); total (67/76)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 146 total aligned reads displayed.

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

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

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

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      >  2:1939055/1‑93 (MQ=255)

   atGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCt                                                                                      >  2:105252/1‑93 (MQ=255)

     gACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTAt                                                                                    <  2:673909/93‑1 (MQ=255)

     gACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTAt                                                                                    <  1:1830869/93‑1 (MQ=255)

     gACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTAt                                                                                    <  1:4640730/93‑1 (MQ=255)

     gACGATCACAAATGGGTTGTGCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTAt                                                                                    <  2:4017980/93‑1 (MQ=255)

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

       cGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGc                                                                                  >  2:1166530/1‑93 (MQ=255)

       cGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGc                                                                                  >  1:4545890/1‑93 (MQ=255)

        gATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGcc                                                                                 <  2:920853/93‑1 (MQ=255)

        gATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGcc                                                                                 <  2:4189709/93‑1 (MQ=255)

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

          tCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGTTGGTTTAGGCGGCAAAGAGATGCCGTTCTAGGAGAGCTTCTATGCCgg                                                                               >  1:1209808/1‑93 (MQ=255)

           cacaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggt                                                                              >  2:3106241/1‑93 (MQ=255)

           cacaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggt                                                                              >  1:2606810/1‑93 (MQ=255)

           cacaAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGGGATGGCTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggg                                                                              >  2:395703/1‑92 (MQ=255)

              aaaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCggtggt                                                                           >  2:2988228/1‑93 (MQ=255)

               aaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGtt                                                                          <  1:4559895/93‑1 (MQ=255)

                 tGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTcc                                                                        <  2:2190426/93‑1 (MQ=255)

                 tGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGGGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTcc                                                                        >  1:3660861/1‑93 (MQ=255)

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

                     ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAgta                                                                    >  2:3324884/1‑91 (MQ=255)

                       gttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCAcc                                                                  <  1:2715432/93‑1 (MQ=255)

                       gttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCAcc                                                                  >  2:3408419/1‑93 (MQ=255)

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

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

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

                         tCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgt                                                                <  1:613966/93‑1 (MQ=255)

                          cTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCTGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtg                                                               <  1:1420543/93‑1 (MQ=255)

                           tGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgc                                                              >  2:2730963/1‑93 (MQ=255)

                            ggggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcg                                                             >  1:1108120/1‑93 (MQ=255)

                            ggggCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACcgtgcg                                                             >  2:3503898/1‑93 (MQ=255)

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

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

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

                                  tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCtt                                                       <  1:1906889/93‑1 (MQ=255)

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

                                   aCCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTc                                                      <  1:1021982/93‑1 (MQ=255)

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

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

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

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

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

                                        cTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTc                                                 >  2:1583949/1‑93 (MQ=255)

                                         tGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTcc                                                >  1:4173613/1‑93 (MQ=255)

                                         tGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTcc                                                >  2:2551228/1‑93 (MQ=255)

                                          ggggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCa                                               >  2:1513385/1‑93 (MQ=255)

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

                                            ggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAAt                                             <  1:3028903/93‑1 (MQ=255)

                                            ggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAAt                                             <  1:2861936/93‑1 (MQ=255)

                                            ggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAAt                                             >  2:3889525/1‑93 (MQ=255)

                                            ggTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAAt                                             >  1:2574882/1‑93 (MQ=255)

                                              ttATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATAc                                           >  1:4145962/1‑93 (MQ=255)

                                              ttATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATAc                                           >  1:572198/1‑93 (MQ=255)

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

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

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

                                                     gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGt                                    <  1:95409/93‑1 (MQ=255)

                                                      aTGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTc                                   >  1:4282949/1‑93 (MQ=255)

                                                      aTGGGTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTTCGTGGCTTCCAGTCCAATACCATTGGTc                                   <  2:1325509/93‑1 (MQ=255)

                                                       tGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTcc                                  <  1:3577682/93‑1 (MQ=255)

                                                       tGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTcc                                  <  1:3860372/93‑1 (MQ=255)

                                                        ggTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCg                                 <  1:3081218/93‑1 (MQ=255)

                                                        ggTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCg                                 <  2:2913222/93‑1 (MQ=255)

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

                                                         gTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGa                                <  1:4239905/93‑1 (MQ=255)

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

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

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

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

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

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

                                                              ggcggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCa                           <  2:830123/93‑1 (MQ=255)

                                                               gcgacAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAg                          <  2:2126388/93‑1 (MQ=255)

                                                                 ggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                        <  1:1190273/93‑1 (MQ=255)

                                                                 ggcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                        <  1:2926593/93‑1 (MQ=255)

                                                                  gcAAAGAGATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGttt                       >  2:2294094/1‑93 (MQ=255)

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

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

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

                                                                      agagATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACtt                   <  1:328613/93‑1 (MQ=255)

                                                                      agagATGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACtt                   >  2:2694123/1‑93 (MQ=255)

                                                                       gagaTGCCGTTCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTc                  <  1:4065216/93‑1 (MQ=255)

                                                                       gagaTGCCGTTCTACGAGAGCTGCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTc                  <  2:634008/93‑1 (MQ=255)

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

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

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

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

                                                                                 tCTACGAGAGCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAgg        >  2:2686743/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                                                       agagCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTa  >  2:4391095/1‑93 (MQ=255)

                                                                                       agagCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTa  >  2:1725213/1‑93 (MQ=255)

TCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTA  >  NC_000913/199789‑199968

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