BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,610,025	A→C	D24A (GAC→GCC)	dcrB →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,610,025	0	A	C	98.7%	520.3 / ‑2.3	152	D24A (GAC→GCC)	dcrB	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (83/65); total (86/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 155 total aligned reads displayed.

TGGCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGACACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTGCCA  >  NC_000913/3609928‑3610119

tGGCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGcca                                                                                              >  1:182911/1‑100 (MQ=255)

 ggCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGccac                                                                                             <  1:1738615/100‑1 (MQ=255)

  gCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACt                                                                                            <  1:14684/100‑1 (MQ=255)

  gCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACt                                                                                            >  2:1328302/1‑100 (MQ=255)

     gTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAAc                                                                                         <  1:2146657/100‑1 (MQ=255)

     gTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAAc                                                                                         <  1:3006866/100‑1 (MQ=255)

     gTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAAc                                                                                         >  2:3973148/1‑100 (MQ=255)

      tttATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACg                                                                                        <  2:1561809/100‑1 (MQ=255)

      tttATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACg                                                                                        <  2:3838591/100‑1 (MQ=255)

      tttATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACg                                                                                        <  1:2334234/100‑1 (MQ=255)

      tttATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACg                                                                                        <  1:2099688/100‑1 (MQ=255)

       ttATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGc                                                                                       <  2:222223/100‑1 (MQ=255)

        tATCGTAAGGAAGAAGTCGATGTGCTATTTGGTTAAAGATGTCGGAATAGGCGGGCTGGTTATGGGGCTTGCGGCCTGTGATGAGAAAGGCATTAACGCt                                                                                      >  2:683402/1‑100 (MQ=255)

         aTCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTa                                                                                     >  1:627913/1‑100 (MQ=255)

         aTCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTa                                                                                     >  1:1226221/1‑100 (MQ=255)

          tCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTAc                                                                                    >  2:1122571/1‑100 (MQ=255)

             aaaGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGgc                                                                                 <  2:2384582/100‑1 (MQ=255)

               aGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGgcgc                                                                               >  1:2045119/1‑100 (MQ=255)

                ggAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCa                                                                              >  2:3532964/1‑100 (MQ=255)

                 gaagaaGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAg                                                                             <  2:3575644/100‑1 (MQ=255)

                 gaagaaGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAg                                                                             >  2:2382525/1‑100 (MQ=255)

                 gaagaaGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAg                                                                             >  1:2913968/1‑100 (MQ=255)

                   agaagTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAggg                                                                           <  1:1791711/100‑1 (MQ=255)

                      agTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGGGGCCTGTGATGATAAAGGCACTAACGCTACGGCGCAGGGTTc                                                                        >  2:3871065/1‑100 (MQ=255)

                      agTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTc                                                                        <  1:439033/100‑1 (MQ=255)

                      agTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTc                                                                        >  2:2096403/1‑100 (MQ=255)

                      agTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTc                                                                        >  2:1068511/1‑100 (MQ=255)

                       gTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCg                                                                       >  1:1398613/1‑100 (MQ=255)

                        tCAATGCTCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGTCTTGCGGCCTGTGAGGATAAAGCCACTGACGCTACGGCGCAGGGCTCgg                                                                      >  2:1399755/1‑100 (MQ=255)

                          aaTGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTc                                                                    >  1:739911/1‑100 (MQ=255)

                            tGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTcgc                                                                  >  2:3065258/1‑100 (MQ=255)

                            tGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTcgc                                                                  >  1:219743/1‑100 (MQ=255)

                             gcgcAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTcgcg                                                                 >  1:2684592/1‑100 (MQ=255)

                             gcgcAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTcgcg                                                                 <  1:2496643/100‑1 (MQ=255)

                             gcgcAATCTGGTTAAATATGTCGGAATTGGCCGGCTGGTTATGGGGTTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTcgcg                                                                 <  2:1922590/100‑1 (MQ=255)

                              cgcAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCgg                                                                <  2:3925771/100‑1 (MQ=255)

                              cgcAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCgg                                                                <  1:991159/100‑1 (MQ=255)

                              cgcAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCgg                                                                >  2:719207/1‑100 (MQ=255)

                                cAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGaa                                                              >  1:796430/1‑100 (MQ=255)

                                cAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCCGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGaa                                                              <  2:761137/100‑1 (MQ=255)

                                  aTCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAg                                                            <  1:687862/100‑1 (MQ=255)

                                  aTCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAg                                                            <  1:3312121/100‑1 (MQ=255)

                                     tGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTaa                                                         >  1:3730468/1‑100 (MQ=255)

                                     tGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGGGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTaa                                                         >  1:699706/1‑100 (MQ=255)

                                        ttAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGc                                                      <  2:724714/100‑1 (MQ=255)

                                         tAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCt                                                     >  1:2225750/1‑100 (MQ=255)

                                          aaaTATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTa                                                    >  1:1522277/1‑100 (MQ=255)

                                            atatGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTAcc                                                  <  1:990210/100‑1 (MQ=255)

                                                gTCGGAATTGGTCTGCTGGTTATGGGGCCTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTCCGCTCGCGGAAAGTAACGCTACCgggc                                              >  2:3295968/1‑99 (MQ=255)

                                                gTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACGGGGa                                              >  2:3796903/1‑100 (MQ=255)

                                                gTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGa                                              >  1:2132248/1‑100 (MQ=255)

                                                 tCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGaa                                             >  1:801689/1‑100 (MQ=255)

                                                  cGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAAt                                            >  2:2653431/1‑100 (MQ=255)

                                                   ggAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATc                                           <  1:1683243/100‑1 (MQ=255)

                                                   ggAATTGGCCGGCTGGTTATGGGGCTTGCGGACTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTGCGGTCGCGGAAAGTAACGCTACCGGGAATc                                           >  2:446302/1‑100 (MQ=255)

                                                    gAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATcc                                          <  2:158960/100‑1 (MQ=255)

                                                     aaTTGGCCTGCTGGTTATGGGGCTTGCGGCCTGGGATGATAAAGCCACTAACGCGACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATccc                                         <  1:4056739/100‑1 (MQ=255)

                                                        tGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTc                                      <  2:216188/100‑1 (MQ=255)

                                                         ggCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCa                                     >  1:317141/1‑100 (MQ=255)

                                                         ggCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCa                                     >  1:3680569/1‑100 (MQ=255)

                                                         ggCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCa                                     <  1:4039740/100‑1 (MQ=255)

                                                         ggCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCa                                     <  2:1710746/100‑1 (MQ=255)

                                                         ggCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCa                                     <  2:1856214/100‑1 (MQ=255)

                                                          gCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCaa                                    >  1:2350750/1‑100 (MQ=255)

                                                          gCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCaa                                    >  2:2863683/1‑100 (MQ=255)

                                                          gCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCaa                                    >  2:706428/1‑100 (MQ=255)

                                                            ctgctgGTTATGGGGCTTGCGGCCTGTGATGATGAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAAcc                                  >  2:2538105/1‑100 (MQ=255)

                                                            ctgctgGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAAcc                                  <  1:2707609/100‑1 (MQ=255)

                                                             tgctgGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCt                                 <  2:2821066/100‑1 (MQ=255)

                                                             tgctgGTTATGGGGCTTGCGGACTGTGATGATAAAGCCGCTAACGCTACGGCGCGCGGATCGGTCTCGGAAAGTAACGCTACCGCGAATCCCGTCAACCt                                 >  2:3810370/1‑100 (MQ=255)

                                                              gctgGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTg                                >  1:3384313/1‑100 (MQ=255)

                                                                tgGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCt                              <  1:954590/100‑1 (MQ=255)

                                                                tgGTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCt                              <  1:2616873/100‑1 (MQ=255)

                                                                 ggTTATGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCtt                             >  1:924909/1‑100 (MQ=255)

                                                                 cgatctGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCtt                             <  2:2889489/95‑1 (MQ=255)

                                                                  gttggggggCTTGCGGCCTGTGATGGTAGAGCCACGAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTg                            <  2:3285477/95‑1 (MQ=255)

                                                                      tGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATgg                        >  1:1261081/1‑100 (MQ=255)

                                                                      tGGGGCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATgg                        <  1:687555/100‑1 (MQ=255)

                                                                       ggggCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGc                       >  1:2512860/1‑100 (MQ=255)

                                                                       ggggCTTGCGGCCTGTGATGAGAAAGCAACTATCGGTACGGCGCAGGGTTGGGTAGGGGAAAGTAACGTTACCGGGAATCCCGTCGCCCTGCTTTATGGc                       <  2:27883/100‑1 (MQ=255)

                                                                        gggCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCa                      <  2:587291/100‑1 (MQ=255)

                                                                        gggCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCa                      <  2:3649381/100‑1 (MQ=255)

                                                                         ggCTTGCGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCaa                     >  1:2640388/1‑100 (MQ=255)

                                                                               cGGCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCaagttaag               >  1:4050222/1‑100 (MQ=255)

                                                                                ggCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCaagttaaag              >  2:2017497/1‑98 (MQ=255)

                                                                                 gCCTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCaagttaagtt             >  1:536990/1‑100 (MQ=255)

                                                                                  ccTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGttt            >  1:3820483/1‑100 (MQ=255)

                                                                                   cTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtc           <  2:3028925/100‑1 (MQ=255)

                                                                                   cTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtc           <  1:602760/100‑1 (MQ=255)

                                                                                   cTGTGATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtc           >  2:1940798/1‑100 (MQ=255)

                                                                                     gtgATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtctc         >  2:2635163/1‑100 (MQ=255)

                                                                                     gtgATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtctc         >  2:1854894/1‑100 (MQ=255)

                                                                                     gtgATGATAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTtctc         >  2:3014865/1‑100 (MQ=255)

                                                                                      tgatgaTAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCg        >  1:3620581/1‑100 (MQ=255)

                                                                                       gatgatAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGc       >  1:1355720/1‑100 (MQ=255)

                                                                                         tgatAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTg     >  1:2020662/1‑100 (MQ=255)

                                                                                         tgatAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTg     <  1:1763489/100‑1 (MQ=255)

                                                                                         tgatAAAGCCACTAACGCTACGGCGCAGGGTTCGGTAGCGGAAAGTAACGCTACCGGGAATCCCGTCGACCTGCTTGAGGGCAAGTTAAGTTTCTCGCTg     <  2:3498399/100‑1 (MQ=255)

                                                                                           atAAAGCCACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTGcc   <  1:892070/100‑1 (MQ=255)

                                                                                            tAAAGCCACTAACGCTACGGTGCAGGGGTCGGTCGCGGAAAGTAACGCGACCGGTAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTGcga  >  2:3983267/1‑98 (MQ=255)

TGGCCGTTTATCGAAAGGAAGAAGTCAATGCGCAATCTGGTTAAATATGTCGGAATTGGCCTGCTGGTTATGGGGCTTGCGGCCTGTGATGATAAAGACACTAACGCTACGGCGCAGGGTTCGGTCGCGGAAAGTAACGCTACCGGGAATCCCGTCAACCTGCTTGATGGCAAGTTAAGTTTCTCGCTGCCA  >  NC_000913/3609928‑3610119

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