BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,391	+C	coding (474/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,391	1	.	C	98.2%	321.0 / ‑2.4	108	K158N (AAG→AAC)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (35/71); minor base T (2/0); total (37/71)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 139 total aligned reads displayed.

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAg                                                                               <  1:5632773/93‑1 (MQ=255)

cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAg                                                                               <  2:493051/93‑1 (MQ=255)

cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAg                                                                               <  2:2347803/93‑1 (MQ=255)

cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑GTTCAAAAAGCCATGCGCTGGAAg                                                                               >  2:2825328/1‑93 (MQ=255)

 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGGAgc                                                                              >  2:6624528/1‑92 (MQ=255)

 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  1:6231040/93‑2 (MQ=255)

 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgc                                                                              >  2:1620044/1‑92 (MQ=255)

 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑GTTCAAAAAGCCATGCGCTGGGAg‑g                                                                             >  1:1501405/1‑93 (MQ=255)

  cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                                             <  1:3847176/93‑3 (MQ=255)

  cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                                             <  2:7316603/93‑3 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTTGAAgcgg                                                                            >  2:2723678/1‑90 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  2:6620098/93‑4 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  1:2588691/93‑4 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  1:3139574/93‑4 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  1:1678421/93‑4 (MQ=255)

   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑GTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            >  1:825140/1‑90 (MQ=255)

    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  1:6940488/93‑5 (MQ=255)

    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  2:6272242/93‑5 (MQ=255)

    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  1:2023289/93‑5 (MQ=255)

    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  2:5206941/93‑5 (MQ=255)

    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcggg                                                                           >  1:6786195/1‑89 (MQ=255)

         gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTTA‑AAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                  >  1:6021564/1‑93 (MQ=255)

          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 <  2:3476810/93‑1 (MQ=255)

          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 <  2:1756023/93‑1 (MQ=255)

          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 <  2:5876197/93‑1 (MQ=255)

          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑GTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 >  2:7825112/1‑93 (MQ=255)

           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  2:6168524/1‑92 (MQ=255)

            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGGAGCGGTGGCG‑‑‑‑aaga                                                               >  1:3816411/1‑91 (MQ=255)

            aTCGCGTCAACGTCAGTGTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               <  2:5679147/93‑3 (MQ=255)

             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              <  2:6345123/93‑4 (MQ=255)

              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             >  2:6900205/1‑89 (MQ=255)

              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  1:6645153/93‑5 (MQ=255)

              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  2:7537857/93‑5 (MQ=255)

              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  1:6133067/93‑5 (MQ=255)

              cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  1:7277068/93‑5 (MQ=255)

                                              gcTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:2469917/92‑1 (MQ=255)

                                               cTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:1059237/93‑2 (MQ=255)

                                               cTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:6824227/93‑2 (MQ=255)

                                                tGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:5036371/93‑3 (MQ=255)

                                                tGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:6350927/93‑3 (MQ=255)

                                                 ggAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  2:137114/1‑90 (MQ=255)

                                                 ggAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  2:1261370/1‑90 (MQ=255)

                                                  gAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                             >  1:2317758/1‑89 (MQ=255)

                                                  gAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                             >  2:5683019/1‑89 (MQ=255)

                                                        aTGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  2:3179687/1‑93 (MQ=255)

                                                        aTGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:3772951/93‑1 (MQ=255)

                                                        aTGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  2:502144/1‑93 (MQ=255)

                                                        aTGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  1:5170210/93‑1 (MQ=255)

                                                         tGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:3177291/93‑1 (MQ=255)

                                                         tGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:2285960/93‑1 (MQ=255)

                                                         tGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:7051633/93‑1 (MQ=255)

                                                          ggcgggcgGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:4703577/93‑1 (MQ=255)

                                                          ggcgggcgGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:5824838/93‑1 (MQ=255)

                                                            cgggcggCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  1:2852182/93‑1 (MQ=255)

                                                             gggcggCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:4150200/93‑1 (MQ=255)

                                                             gggcggCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  >  1:7618833/1‑93 (MQ=255)

                                                             gggcggCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:6448749/93‑1 (MQ=255)

                                                             gggcggCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:6162743/93‑1 (MQ=255)

                                                              ggcggcG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:5450451/93‑1 (MQ=255)

                                                              ggcggcG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:2454381/93‑1 (MQ=255)

                                                              ggcggcG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:7078413/93‑1 (MQ=255)

                                                               gcggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  2:2387781/1‑93 (MQ=255)

                                                                cggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:1673796/93‑1 (MQ=255)

                                                                cggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:3005645/93‑1 (MQ=255)

                                                                cggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:184452/93‑1 (MQ=255)

                                                                cggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:4611868/93‑1 (MQ=255)

                                                                 ggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  1:4823526/1‑93 (MQ=255)

                                                                  gcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:3319815/1‑93 (MQ=255)

                                                                  gcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacc             >  1:5719364/1‑92 (MQ=255)

                                                                   cg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:3574992/93‑1 (MQ=255)

                                                                   cg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  2:3399617/1‑93 (MQ=255)

                                                                    g‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcg           >  1:1056564/1‑92 (MQ=255)

                                                                    g‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  1:2079819/93‑1 (MQ=255)

                                                                    g‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:6601284/93‑1 (MQ=255)

                                                                    g‑TTTCAAAAGGCCATGCGCTGGAAGCGGCGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcg           >  2:5160899/1‑92 (MQ=255)

                                                                       ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa         <  1:2136365/89‑1 (MQ=255)

                                                                        tcaaaaGGCCATGCGCTGGAAGTGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  1:4635586/4‑93 (MQ=255)

                                                                        tcaaaaGGCCATGCGCTGGAAGTGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  1:1052812/4‑93 (MQ=255)

                                                                        tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:1515050/4‑93 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:7352998/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:6859924/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  1:553209/3‑93 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:5401855/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:4920707/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:7277919/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:1777780/91‑1 (MQ=255)

                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:4433613/91‑1 (MQ=255)

                                                                          aaaaGGCCATGCGCTGGAAGCGGTTGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  1:975371/92‑1 (MQ=255)

                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4664414/2‑93 (MQ=255)

                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGCCGCCGAAATg     >  2:2449115/1‑93 (MQ=255)

                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:3147096/1‑93 (MQ=255)

                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:7112816/1‑93 (MQ=255)

                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     <  1:6346996/93‑1 (MQ=255)

                                                                            aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  2:7165808/1‑93 (MQ=255)

                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:48196/93‑1 (MQ=255)

                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:4224177/1‑93 (MQ=255)

                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:3129636/1‑93 (MQ=255)

                                                                              ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:6027781/93‑1 (MQ=255)

                                                                              ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:642075/93‑1 (MQ=255)

                                                                              ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:2850945/93‑1 (MQ=255)

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

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