BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,865,091	T→C	T162A (ACC→GCC)	rpoS ←	sigma S factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,865,091	0	T	C	88.1%	657.7 / 71.9	261	T162A (ACC→GCC)	rpoS	sigma S factor
Reads supporting (aligned to +/- strand): ref base T (9/22); new base C (96/134); total (105/156)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 268 total aligned reads displayed.

TTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTA  >  NC_000913/2864999‑2865183

ttCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGt                                                                                              >  2:3476479/1‑93 (MQ=255)

  cATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGctt                                                                                            <  2:5919214/93‑1 (MQ=255)

  cATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGctt                                                                                            <  2:4227926/93‑1 (MQ=255)

  cATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGctt                                                                                            >  1:7299727/1‑93 (MQ=255)

   aTGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTg                                                                                           >  2:5368372/1‑93 (MQ=255)

   aTGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTg                                                                                           <  2:3703904/93‑1 (MQ=255)

   aTGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTg                                                                                           <  1:6428356/93‑1 (MQ=255)

    tGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTgg                                                                                          <  2:586374/93‑1 (MQ=255)

       tCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTc                                                                                       >  1:5054124/1‑93 (MQ=255)

       tCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTc                                                                                       <  1:3691331/93‑1 (MQ=255)

       tCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTc                                                                                       <  2:156329/93‑1 (MQ=255)

          aGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATa                                                                                    <  2:3305991/93‑1 (MQ=255)

           gCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATaa                                                                                   <  2:2059986/93‑1 (MQ=255)

           gCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATaa                                                                                   >  1:40163/1‑93 (MQ=255)

            cTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAAt                                                                                  >  1:1026046/1‑93 (MQ=255)

            cTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAAt                                                                                  >  1:763729/1‑93 (MQ=255)

             ttATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATc                                                                                 <  1:6888354/93‑1 (MQ=255)

             ttATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATc                                                                                 >  2:6369911/1‑93 (MQ=255)

               aTGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGc                                                                               <  1:5352906/93‑1 (MQ=255)

               aTGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGc                                                                               >  1:5663080/1‑93 (MQ=255)

                tGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGcc                                                                              <  1:2927378/93‑1 (MQ=255)

                 gggACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGccc                                                                             <  2:1223192/93‑1 (MQ=255)

                  ggACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCg                                                                            <  1:6425585/93‑1 (MQ=255)

                     cAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTc                                                                         >  2:1919271/1‑93 (MQ=255)

                     cAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTc                                                                         <  1:532065/93‑1 (MQ=255)

                       aCTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCaa                                                                       >  1:4149352/1‑93 (MQ=255)

                        cTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAAt                                                                      <  2:1413326/93‑1 (MQ=255)

                        cTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAAt                                                                      >  1:7155287/1‑93 (MQ=255)

                        cTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAAt                                                                      <  2:2286069/93‑1 (MQ=255)

                        cTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAAt                                                                      <  1:1083897/93‑1 (MQ=255)

                         tCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATc                                                                     >  2:6071783/1‑93 (MQ=255)

                          cACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCg                                                                    >  2:7299415/1‑93 (MQ=255)

                          cACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCg                                                                    <  2:3328916/93‑1 (MQ=255)

                          cACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCg                                                                    >  1:3071366/1‑93 (MQ=255)

                            cGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTc                                                                  <  2:7340004/93‑1 (MQ=255)

                             gTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCt                                                                 <  2:3612907/93‑1 (MQ=255)

                              tGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTg                                                                >  1:459812/1‑93 (MQ=255)

                              tGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTg                                                                <  2:6848837/93‑1 (MQ=255)

                                cGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGc                                                              >  2:5039121/1‑93 (MQ=255)

                                 ggTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCg                                                             <  2:4876735/93‑1 (MQ=255)

                                  gTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGcta                                                            >  2:7189414/1‑91 (MQ=255)

                                      gCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATcc                                                        >  2:4345467/1‑93 (MQ=255)

                                       cAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACTGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATcca                                                       <  2:3139019/93‑1 (MQ=255)

                                       cAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATcca                                                       <  2:1667786/93‑1 (MQ=255)

                                       cAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATcca                                                       <  1:3189089/93‑1 (MQ=255)

                                        aGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATccac                                                      <  1:4351068/93‑1 (MQ=255)

                                         ggTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATccacc                                                     <  1:1387526/93‑1 (MQ=255)

                                         ggTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATccacc                                                     >  2:6820438/1‑93 (MQ=255)

                                         ggTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATccacc                                                     <  1:2233000/93‑1 (MQ=255)

                                          gTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATccacca                                                    >  1:3769373/1‑93 (MQ=255)

                                             aaCGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTGTGGCGAATCCACCAGGt                                                 >  2:7468416/1‑93 (MQ=255)

                                              aCGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGtt                                                <  2:1935162/93‑1 (MQ=255)

                                              aCGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGtt                                                <  2:172884/93‑1 (MQ=255)

                                              aCGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGtt                                                >  1:5483328/1‑93 (MQ=255)

                                               cGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTg                                               >  2:5593631/1‑93 (MQ=255)

                                                gTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGc                                              >  2:6690814/1‑93 (MQ=255)

                                                gTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGc                                              >  2:5206176/1‑93 (MQ=255)

                                                 ttCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCg                                             <  1:1471499/93‑1 (MQ=255)

                                                 ttCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCg                                             >  1:4813304/1‑93 (MQ=255)

                                                 ttCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCg                                             >  1:6307410/1‑93 (MQ=255)

                                                  tCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGt                                            >  1:254489/1‑93 (MQ=255)

                                                   cAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTa                                           <  2:2421982/93‑1 (MQ=255)

                                                   cAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTa                                           >  1:2853257/1‑93 (MQ=255)

                                                   cAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTa                                           <  2:6214503/93‑1 (MQ=255)

                                                    aGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTAt                                          >  1:38660/1‑93 (MQ=255)

                                                    aGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTAt                                          >  2:4355789/1‑93 (MQ=255)

                                                      cTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGt                                        <  2:1445140/93‑1 (MQ=255)

                                                       tCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGtt                                       >  2:4546562/1‑93 (MQ=255)

                                                       tCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGCATGtt                                       <  1:1368767/93‑1 (MQ=255)

                                                        ccTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTg                                      <  2:4710535/93‑1 (MQ=255)

                                                        ccTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTg                                      >  2:2546771/1‑93 (MQ=255)

                                                         cTTTACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTga                                     >  2:689702/1‑93 (MQ=255)

                                                           ttACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTgaga                                   <  2:1122135/93‑1 (MQ=255)

                                                            tACGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGaa                                  <  2:4458591/93‑1 (MQ=255)

                                                              cGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGc                                <  2:553432/93‑1 (MQ=255)

                                                              cGATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGc                                <  2:523879/93‑1 (MQ=255)

                                                               gATGTGAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCg                               >  1:1826444/1‑93 (MQ=255)

                                                                aTGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCgg                              <  2:6989097/93‑1 (MQ=255)

                                                                 tgtgAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGa                             <  2:2317184/93‑1 (MQ=255)

                                                                  gtgAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGaa                            >  1:1672643/1‑93 (MQ=255)

                                                                   tgAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGaaa                           >  1:3387370/1‑93 (MQ=255)

                                                                    gAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAc                          <  2:6559958/93‑1 (MQ=255)

                                                                    gAATCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAc                          <  1:3233232/93‑1 (MQ=255)

                                                                      aTCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCa                        >  1:6765735/1‑93 (MQ=255)

                                                                      aTCGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCa                        <  1:5756266/93‑1 (MQ=255)

                                                                        cGGCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACg                      <  1:6651941/93‑1 (MQ=255)

                                                                         ggCAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGt                     >  1:128999/1‑93 (MQ=255)

                                                                           cAAACGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTc                   >  1:2751709/1‑93 (MQ=255)

                                                                            aaaCGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTcc                  <  1:2329607/93‑1 (MQ=255)

                                                                             aaCGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCg                 <  2:394909/93‑1 (MQ=255)

                                                                             aaCGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCg                 <  2:2886837/93‑1 (MQ=255)

                                                                              aCGAATAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCgg                <  2:1337508/93‑1 (MQ=255)

                                                                                 aaTAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTc             >  2:5744452/1‑93 (MQ=255)

                                                                                 aaTAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTc             >  2:6266932/1‑93 (MQ=255)

                                                                                 aaTAGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTc             <  2:2685821/93‑1 (MQ=255)

                                                                                  aTAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCa            >  2:247673/1‑93 (MQ=255)

                                                                                    aGTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCaaa          >  1:2607321/1‑93 (MQ=255)

                                                                                     gTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAAc         >  1:2123765/1‑93 (MQ=255)

                                                                                     gTACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAAc         <  1:355831/93‑1 (MQ=255)

                                                                                      tACGGGCTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACt        <  1:1967429/93‑1 (MQ=255)

TTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTA  >  NC_000913/2864999‑2865183

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