BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,661,865	T→A	L135Q (CTG→CAG)	suhB →	SuhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,661,865	0	T	A	99.6%	877.5 / ‑4.5	235	L135Q (CTG→CAG)	suhB	SuhB
Reads supporting (aligned to +/- strand): ref base T (1/0); new base A (132/102); total (133/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 247 total aligned reads displayed.

TTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTACG  >  NC_000913/2661775‑2661954

ttCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGa                                                                                         >  1:4808620/1‑93 (MQ=255)

 tCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGaa                                                                                        <  2:1056204/93‑1 (MQ=255)

  cGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCTTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAAc                                                                                       <  1:2279861/93‑1 (MQ=255)

  cGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAAc                                                                                       <  1:504549/93‑1 (MQ=255)

  cGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAAc                                                                                       <  1:756016/93‑1 (MQ=255)

  cGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAAc                                                                                       <  2:2469579/93‑1 (MQ=255)

  cGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAAc                                                                                       >  2:600831/1‑93 (MQ=255)

   gTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACg                                                                                      <  1:3489003/93‑1 (MQ=255)

   gTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACg                                                                                      >  2:1765459/1‑93 (MQ=255)

    tATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACgg                                                                                     >  2:562650/1‑93 (MQ=255)

     aTCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGc                                                                                    <  2:6713389/93‑1 (MQ=255)

      tCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCt                                                                                   <  1:1826180/93‑1 (MQ=255)

       cAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTa                                                                                  >  2:777510/1‑93 (MQ=255)

       cAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTa                                                                                  <  1:1703052/93‑1 (MQ=255)

        aaaGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTAc                                                                                 >  2:2160268/1‑93 (MQ=255)

         aaGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGACTatc                                                                                >  1:668802/1‑91 (MQ=255)

          aGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCg                                                                               >  2:132819/1‑93 (MQ=255)

           ggCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGa                                                                              >  2:571459/1‑93 (MQ=255)

           ggCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGa                                                                              >  2:6597208/1‑93 (MQ=255)

            gCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGAc                                                                             <  1:981035/93‑1 (MQ=255)

                 aCCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTgcgc                                                                        <  2:5263585/93‑1 (MQ=255)

                 aCCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTgcgc                                                                        >  1:3857721/1‑93 (MQ=255)

                 aCCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTgcgc                                                                        >  2:1692355/1‑93 (MQ=255)

                  ccGGAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGcgcg                                                                       <  2:5034658/93‑1 (MQ=255)

                    gAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGCCTGCGCGgc                                                                     >  2:73659/1‑93 (MQ=255)

                     aaGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGgca                                                                    >  1:2911621/1‑93 (MQ=255)

                      aGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGgcag                                                                   >  2:413287/1‑93 (MQ=255)

                      aGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGgcag                                                                   >  1:4721713/1‑93 (MQ=255)

                       gTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGgcagc                                                                  <  1:6482837/93‑1 (MQ=255)

                          gCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCAcc                                                               >  1:1293212/1‑93 (MQ=255)

                          gCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCAcc                                                               >  1:5143985/1‑93 (MQ=255)

                           cTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCg                                                              >  1:5917714/1‑93 (MQ=255)

                              tgGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTc                                                           >  2:448066/1‑93 (MQ=255)

                              tgGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTc                                                           <  2:4308457/93‑1 (MQ=255)

                               ggTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTcg                                                          <  2:3970968/93‑1 (MQ=255)

                               ggTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTcg                                                          <  1:1953447/93‑1 (MQ=255)

                               ggTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTcg                                                          <  2:6393410/93‑1 (MQ=255)

                                gTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTcgc                                                         <  2:5429695/93‑1 (MQ=255)

                                 tataCGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTcccg                                                        <  1:4936905/93‑1 (MQ=255)

                                   taCGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGAt                                                      <  1:3143271/93‑1 (MQ=255)

                                    aCGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGAtc                                                     <  1:4799789/93‑1 (MQ=255)

                                     cGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGAtct                                                    <  1:5427516/93‑1 (MQ=255)

                                     cGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGAtct                                                    >  2:5964197/1‑93 (MQ=255)

                                      caTCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGAtctc                                                   <  2:1957711/92‑1 (MQ=255)

                                       aTCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCg                                                  <  1:5248314/93‑1 (MQ=255)

                                       aTCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCg                                                  <  2:1192687/93‑1 (MQ=255)

                                       aTCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCg                                                  >  1:1781806/1‑93 (MQ=255)

                                        tCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGa                                                 >  1:4593009/1‑93 (MQ=255)

                                          cTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACg                                               <  1:730449/93‑1 (MQ=255)

                                           tATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACgg                                              >  1:2424807/1‑93 (MQ=255)

                                            aTGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGt                                             >  1:1441452/1‑93 (MQ=255)

                                            aTGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGt                                             <  1:4905506/93‑1 (MQ=255)

                                              gCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTat                                           >  2:3041807/1‑92 (MQ=255)

                                              gCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTAc                                           <  1:3547202/93‑1 (MQ=255)

                                               cGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACt                                          >  2:25563/1‑93 (MQ=255)

                                                gTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTa                                         >  1:5294184/1‑93 (MQ=255)

                                                 tAACGACCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTAt                                        <  1:6293817/93‑1 (MQ=255)

                                                 tAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTAt                                        >  1:1097841/1‑93 (MQ=255)

                                                 tAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTAt                                        >  2:380855/1‑93 (MQ=255)

                                                 tAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTAt                                        >  2:3480441/1‑93 (MQ=255)

                                                  aaCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTAtt                                       <  2:3536803/93‑1 (MQ=255)

                                                   aCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTc                                      >  2:2282248/1‑93 (MQ=255)

                                                   aCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTc                                      <  2:6477848/93‑1 (MQ=255)

                                                   aCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTc                                      >  2:49530/1‑93 (MQ=255)

                                                   aCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTc                                      <  2:4566963/93‑1 (MQ=255)

                                                   aCGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTc                                      <  1:5826868/93‑1 (MQ=255)

                                                    cGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCt                                     <  2:2828060/93‑1 (MQ=255)

                                                    cGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCt                                     >  2:3581532/1‑93 (MQ=255)

                                                      aaCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTgg                                   <  1:5086880/93‑1 (MQ=255)

                                                      aaCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTgg                                   >  1:2191379/1‑93 (MQ=255)

                                                      aaCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTgg                                   <  1:665394/93‑1 (MQ=255)

                                                      aaCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTgg                                   >  1:5754736/1‑93 (MQ=255)

                                                       aCTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGc                                  >  1:863375/1‑93 (MQ=255)

                                                        cTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCg                                 >  2:5510660/1‑93 (MQ=255)

                                                          gTTCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGAc                               <  2:1237109/93‑1 (MQ=255)

                                                           ttCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGAcc                              >  1:2644892/1‑93 (MQ=255)

                                                           ttCACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGAcc                              <  2:2048851/93‑1 (MQ=255)

                                                             cACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCgg                            <  1:4159942/93‑1 (MQ=255)

                                                             cACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCgg                            <  1:4244479/93‑1 (MQ=255)

                                                             cACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCgg                            >  1:3073017/1‑93 (MQ=255)

                                                             cACCGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCgg                            <  2:2951035/93‑1 (MQ=255)

                                                               ccGCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCt                          >  2:6197315/1‑93 (MQ=255)

                                                                 gCCACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTc                        >  1:1232822/1‑93 (MQ=255)

                                                                   cACTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTccc                      <  2:1969815/93‑1 (MQ=255)

                                                                    aCTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTTTGGCGACCGGCTTcgcg                     >  1:3628253/1‑93 (MQ=255)

                                                                    aCTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCg                     >  2:1120303/1‑93 (MQ=255)

                                                                     cTCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGt                    <  2:3291098/93‑1 (MQ=255)

                                                                      tCGCGGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGtt                   >  2:3102933/1‑93 (MQ=255)

                                                                       cgcgGTCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGGACTATTCTGGCGACCGGCTTCCCGTTc                  >  2:2387881/1‑93 (MQ=255)

                                                                            tCAGGGCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGc             <  1:1561707/93‑1 (MQ=255)

                                                                               gggCGCACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCaaa          >  2:5816046/1‑93 (MQ=255)

                                                                                 gcgcACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAAc        <  1:4876245/93‑1 (MQ=255)

                                                                                  cgcACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGCCCGGCTTCCCGTTCAAAGCAAAACa       <  2:2099522/93‑1 (MQ=255)

                                                                                   gcACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAg      <  2:6084001/93‑1 (MQ=255)

                                                                                   gcACAGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAg      <  1:2784818/93‑1 (MQ=255)

                                                                                    cacaGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGt     >  2:902854/1‑93 (MQ=255)

                                                                                     acaGCAGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTa    <  2:1303600/93‑1 (MQ=255)

                                                                                      cagcagAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTAc   <  1:2578111/93‑1 (MQ=255)

                                                                                      cagcagAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTAc   >  2:883378/1‑93 (MQ=255)

                                                                                      cagcagAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTAc   >  2:5221524/1‑93 (MQ=255)

TTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGCTCGCGATCTCGACGGTACTATTCTGGCGACCGGCTTCCCGTTCAAAGCAAAACAGTACG  >  NC_000913/2661775‑2661954

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