BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,383	T→C	S475S (AGT→AGC)	ytfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,383	0	T	C	99.7%	1030.2 / NA	330	S475S (AGT→AGC)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (121/207); total (123/207)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

TGTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGTATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTGGA  >  NC_000913/4441293‑4441472

tGTTACACGCGGCACGCTGGGCTGTATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCAt                                                                                         <  1:228795/93‑1 (MQ=255)

tGTTACACGAGACACGCTGAGCTGGATTGAAACAGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCAt                                                                                         <  2:7747148/93‑1 (MQ=255)

 gTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCAtt                                                                                        <  1:606040/93‑1 (MQ=255)

   tACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCg                                                                                      <  1:7178718/93‑1 (MQ=255)

   tACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCg                                                                                      >  1:1625950/1‑93 (MQ=255)

   tACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCg                                                                                      <  2:6511671/93‑1 (MQ=255)

    acacGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGt                                                                                     <  1:4608468/93‑1 (MQ=255)

     cacGCGGCACGCTGGGCTGTATTGAAACCGGTGATTGCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                    <  2:5196287/93‑1 (MQ=255)

     cacGCGGCACGCTGGGCTGGATTGAAGCCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                    <  1:6212159/93‑1 (MQ=255)

     cacGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                    >  1:1416434/1‑93 (MQ=255)

     cacGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGAACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                    >  1:6772988/1‑93 (MQ=255)

      acGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATATGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTgg                                                                                   >  1:948547/1‑93 (MQ=255)

        gcgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCt                                                                                 <  1:7598011/93‑1 (MQ=255)

         cgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCta                                                                                <  1:6969795/93‑1 (MQ=255)

          ggCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtac                                                                               <  2:6350649/93‑1 (MQ=255)

          ggCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtac                                                                               <  1:1816678/93‑1 (MQ=255)

           gCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtaca                                                                              <  1:7105247/93‑1 (MQ=255)

              cGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaat                                                                           <  1:1212719/93‑1 (MQ=255)

              cGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaat                                                                           <  2:6480274/93‑1 (MQ=255)

               gctgggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaata                                                                          <  1:6318945/93‑1 (MQ=255)

               gctgggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaata                                                                          <  1:4771307/93‑1 (MQ=255)

                   ggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaatacaaa                                                                      <  1:2777964/93‑1 (MQ=255)

                   ggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaatacaaa                                                                      <  2:6687930/93‑1 (MQ=255)

                    gctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTacaaatacaaat                                                                     <  1:3530895/93‑1 (MQ=255)

                     ctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATc                                                                    <  2:1821818/93‑1 (MQ=255)

                      tggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCt                                                                   <  1:1396739/93‑1 (MQ=255)

                       ggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTa                                                                  <  2:3446327/93‑1 (MQ=255)

                         aTTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATc                                                                <  2:3548707/93‑1 (MQ=255)

                         aTTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATc                                                                <  2:6148104/93‑1 (MQ=255)

                         aTTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATc                                                                <  1:7638451/93‑1 (MQ=255)

                          ttGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCg                                                               >  1:6875799/1‑93 (MQ=255)

                           tGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGc                                                              <  2:1437809/93‑1 (MQ=255)

                              aaCCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTcc                                                           >  2:4476907/1‑93 (MQ=255)

                              aaCCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTcc                                                           <  2:6201006/93‑1 (MQ=255)

                                ccGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGa                                                         >  1:432380/1‑93 (MQ=255)

                                  ggTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGaaa                                                       <  1:7770677/93‑1 (MQ=255)

                                  ggTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGaaa                                                       <  2:1242940/93‑1 (MQ=255)

                                  ggTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGaaa                                                       <  1:4898559/93‑1 (MQ=255)

                                   gTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAAt                                                      >  2:5138980/1‑93 (MQ=255)

                                    tGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATa                                                     >  1:5876565/1‑93 (MQ=255)

                                    tGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATa                                                     >  1:1291582/1‑93 (MQ=255)

                                    tGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATa                                                     <  1:822299/93‑1 (MQ=255)

                                    tGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATa                                                     <  1:3458850/93‑1 (MQ=255)

                                      aTTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACg                                                   <  1:6609497/93‑1 (MQ=255)

                                       tttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCTGGGGCGACCGCAGTATTCGTGTCTACAATTACAAATCTATCGCGCCGAAATACGc                                                  >  2:7012224/1‑93 (MQ=255)

                                       tttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGc                                                  <  2:2576453/93‑1 (MQ=255)

                                        ttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGcc                                                 <  2:1679053/93‑1 (MQ=255)

                                          cGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCaa                                               <  2:875428/93‑1 (MQ=255)

                                           gACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAAc                                              <  1:3047346/93‑1 (MQ=255)

                                           gACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAAc                                              >  1:5558884/1‑93 (MQ=255)

                                            aCAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACg                                             >  2:4287603/1‑93 (MQ=255)

                                            aCAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACg                                             <  2:2067505/93‑1 (MQ=255)

                                             cAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACgg                                            >  2:3736292/1‑93 (MQ=255)

                                              aaaGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGt                                           >  1:1011377/1‑93 (MQ=255)

                                              aaaGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGt                                           <  2:7191119/93‑1 (MQ=255)

                                                 gTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGAc                                        >  2:335259/1‑93 (MQ=255)

                                                   aCCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCt                                      >  1:7916016/1‑93 (MQ=255)

                                                   aCCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCt                                      <  1:4109763/93‑1 (MQ=255)

                                                    ccgccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTg                                     >  1:4923062/1‑93 (MQ=255)

                                                     cgccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGa                                    <  2:5484107/93‑1 (MQ=255)

                                                     cgccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGa                                    <  2:1277558/93‑1 (MQ=255)

                                                       ccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGaaa                                  >  2:794224/1‑93 (MQ=255)

                                                        cgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAg                                 <  1:5519556/93‑1 (MQ=255)

                                                        cgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAg                                 <  2:1138133/93‑1 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                                <  1:3839032/93‑1 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                                <  2:4589715/93‑1 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                                <  2:4063954/93‑1 (MQ=255)

                                                          gATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAggg                               >  1:5934775/1‑93 (MQ=255)

                                                           aTCTGCGTTTCTTCGCCGGGGGCAACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgggg                              <  2:6616125/93‑1 (MQ=255)

                                                            tCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGc                             <  1:7399845/93‑1 (MQ=255)

                                                              tGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCt                           >  1:727804/1‑93 (MQ=255)

                                                                cGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCg                         <  1:3137762/93‑1 (MQ=255)

                                                                 gTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGa                        >  2:3784646/1‑93 (MQ=255)

                                                                 gTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGa                        <  1:5138248/93‑1 (MQ=255)

                                                                  tttCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGaa                       <  2:5930116/93‑1 (MQ=255)

                                                                  tttCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGaa                       >  1:5368475/1‑93 (MQ=255)

                                                                  tttCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGaa                       <  1:3895941/93‑1 (MQ=255)

                                                                    tcttcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGt                     <  1:5704363/93‑1 (MQ=255)

                                                                      ttcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTg                   >  2:5688591/1‑93 (MQ=255)

                                                                       tcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGa                  >  2:7463673/1‑93 (MQ=255)

                                                                       tcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGa                  <  2:2367880/93‑1 (MQ=255)

                                                                          ccGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATaa               >  2:7517963/1‑93 (MQ=255)

                                                                            gggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAAcc             >  2:28366/1‑93 (MQ=255)

                                                                             ggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCg            >  2:3937128/1‑93 (MQ=255)

                                                                             ggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCg            <  1:1540331/93‑1 (MQ=255)

                                                                               ggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGa          <  2:1938742/93‑1 (MQ=255)

                                                                               ggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGa          >  1:2627909/1‑93 (MQ=255)

                                                                                gCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGAt         <  2:304508/93‑1 (MQ=255)

                                                                                gCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGAt         <  1:243793/93‑1 (MQ=255)

                                                                                   gccGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc      <  1:2164273/92‑1 (MQ=255)

                                                                                   aCCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc      >  1:7466848/1‑93 (MQ=255)

                                                                                   aCCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc      >  2:2984896/1‑93 (MQ=255)

                                                                                   aCCGCAGCATTCGTGGCTACAAATACAAAACTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc      >  2:4839561/1‑93 (MQ=255)

                                                                                    ccGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCt     >  2:5788919/1‑93 (MQ=255)

                                                                                    ccGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCt     >  1:6439825/1‑93 (MQ=255)

                                                                                    ccGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCt     >  2:2647376/1‑93 (MQ=255)

                                                                                     cGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTg    <  2:4994593/93‑1 (MQ=255)

                                                                                      gcagcaTTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTgg   >  2:7254440/1‑93 (MQ=255)

                                                                                      gcagcaTTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTgg   >  1:5025886/1‑93 (MQ=255)

                                                                                       cAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTGGa  >  1:2002519/1‑93 (MQ=255)

TGTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGTATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTGGA  >  NC_000913/4441293‑4441472

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