BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,411	+C	coding (494/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,410	1	.	C	99.5%	655.8 / ‑3.9	188	coding (493/516 nt)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (73/114); total (73/115)
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.87e-01

Only 100 of 194 total aligned reads displayed.

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AA‑‑‑‑CT‑GGCTACAG‑C‑TGAAGGATTG‑ATTGC‑GTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑GCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtt                                                                       <  1:5033459/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtt                                                                       >  1:337657/1‑93 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtt                                                                       <  2:2237438/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtt                                                                       <  2:2735328/93‑1 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gttt                                                                      <  2:6914722/93‑2 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gttt                                                                      <  1:5418012/93‑2 (MQ=255)

  tGGAAGTGATGGCGGGCTGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttg                                                                     <  2:7889599/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttg                                                                     <  2:6065576/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttg                                                                     <  2:5872788/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttg                                                                     >  1:4769163/1‑91 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttggt                                                                   <  1:414086/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttggt                                                                   <  1:5795688/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑Gtttggt                                                                   <  2:8637298/93‑5 (MQ=255)

           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGa                                                            <  2:6339319/93‑1 (MQ=255)

           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGa                                                            <  2:5721445/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGAt                                                           <  1:4509590/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGAt                                                           <  2:8755362/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGAt                                                           <  1:8788330/93‑1 (MQ=255)

            ggcgggcgCCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGAt                                                           <  2:421304/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATg                                                          <  2:1244670/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATg                                                          <  1:7587141/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATg                                                          <  1:2829196/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATg                                                          <  1:2455349/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAAGCTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATg                                                          <  1:8011368/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGa                                                         <  2:2568450/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGa                                                         <  2:8090470/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGa                                                         >  1:1004117/1‑93 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGaa                                                        <  1:2848457/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGaa                                                        <  2:497721/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGaa                                                        <  2:8370276/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGaa                                                        <  2:1912031/93‑1 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAac                                                       >  1:5543527/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAac                                                       >  2:3824349/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAacg                                                      >  2:3764502/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAacg                                                      <  1:6357237/93‑1 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAacga                                                     >  2:3745720/1‑93 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAacgac                                                    <  1:852079/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAacgac                                                    <  1:3665392/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGcc                                                 <  1:5184870/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGcc                                                 >  1:403474/1‑93 (MQ=255)

                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGa                                               >  2:7031658/5‑93 (MQ=255)

                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaa                                              >  2:144119/4‑93 (MQ=255)

                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaa                                              >  1:4387540/4‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaaa                                             <  2:412049/91‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaaa                                             >  2:3401072/3‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaaa                                             >  1:3273534/3‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGaaa                                             <  2:5036154/91‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAAt                                            >  2:1868207/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAAt                                            >  2:1302241/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAAt                                            <  2:3778566/92‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAAt                                            >  1:5353936/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAAt                                            >  2:5855352/2‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCGGAAATg                                           >  1:2685785/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATg                                           >  1:8732540/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATg                                           <  1:1458219/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGc                                          >  1:8320329/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGc                                          <  2:4556424/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGc                                          >  2:4465606/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGc                                          <  1:2268510/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCt                                         >  2:26644/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCt                                         <  1:2102219/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCt                                         <  2:5856790/93‑1 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                        <  2:7788741/93‑1 (MQ=255)

                                                                actgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  2:5879707/89‑1 (MQ=255)

                                                                actgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  1:6792944/89‑1 (MQ=255)

                                                                actgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  2:3994431/89‑1 (MQ=255)

                                                                actgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg         <  1:8139781/89‑1 (MQ=255)

                                                                actgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGAATTATggg         >  1:1000870/5‑93 (MQ=255)

                                                                 ttgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTTACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc        <  1:2174924/90‑1 (MQ=255)

                                                                 ctgggCTACAGTC‑TGAAGGAATGTATTGCGGTTG‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc        <  1:2414968/90‑1 (MQ=255)

                                                                  tgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa       <  2:3431550/91‑1 (MQ=255)

                                                                  tgggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa       <  2:1906386/91‑1 (MQ=255)

                                                                   gggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg      >  1:6255832/2‑93 (MQ=255)

                                                                   gggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg      <  1:1340236/92‑1 (MQ=255)

                                                                   gggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg      <  1:7719470/92‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  1:1217457/93‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  2:4277128/93‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  1:6302462/93‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  1:6443101/93‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  2:6623182/93‑1 (MQ=255)

                                                                    ggCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt     <  1:6995804/93‑1 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGTTGAACGACGCCGAAATGCTGTCGATGGC‑GGGGCAAG‑GCTGCATTATGGGCAGTg    >  1:7452218/1‑93 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    >  1:2234361/1‑93 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    >  1:31130/1‑93 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    >  2:6718306/1‑93 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    <  2:4887168/93‑1 (MQ=255)

                                                                     gCTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg    <  1:2656890/93‑1 (MQ=255)

                                                                      cTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGGAGTgc   >  2:1260452/1‑92 (MQ=255)

                                                                      cTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgg   >  2:1260453/1‑93 (MQ=255)

                                                                      cTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   <  2:3629974/93‑2 (MQ=255)

                                                                      cTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   >  2:6885042/1‑92 (MQ=255)

                                                                      cTACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   >  1:6911350/1‑92 (MQ=255)

                                                                      cTACAG‑CGTGAAGGATTGTATTGC‑GTTTGGTGAAGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc   <  2:8646167/93‑2 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAATGCTGCATTATGGGCAGTgcg  >  2:3150070/1‑91 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  >  1:4964817/1‑91 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  <  2:432277/93‑3 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  >  1:19048/1‑91 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  <  1:1945699/93‑3 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg  >  2:5406802/1‑91 (MQ=255)

                                                                       tACAGCC‑TGAAGGATTGTATTGC‑GTTTGGTAACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGCAAG‑GCTGCATTATGGGCAGTgcg  >  2:6011188/1‑91 (MQ=255)

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AA‑‑‑‑CT‑GGCTACAG‑C‑TGAAGGATTG‑ATTGC‑GTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑GCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492

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