BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	547,645	A→G	H22R (CAC→CGC)	b4507 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	547,645	0	A	G	83.5%	590.2 / 108.6	266	H22R (CAC→CGC)	b4507	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (21/23); new base G (109/113); total (130/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 272 total aligned reads displayed.

TAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCA  >  NC_000913/547555‑547735

  gAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctgg                                                                                        <  1:7104620/93‑1 (MQ=255)

    aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctggct                                                                                      >  1:5711942/1‑93 (MQ=255)

    aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctggct                                                                                      >  1:3222340/1‑93 (MQ=255)

    aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggct                                                                                      >  1:1219942/1‑93 (MQ=255)

     cGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACgctggctg                                                                                     <  2:337783/93‑1 (MQ=255)

      gTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctggctgg                                                                                    >  2:4721801/1‑93 (MQ=255)

      gTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctggctgg                                                                                    >  2:3065547/1‑93 (MQ=255)

      gTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGctggctgg                                                                                    >  2:1322519/1‑93 (MQ=255)

        ttGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGAt                                                                                  <  1:6471933/93‑1 (MQ=255)

         tGCAATAAGGGGTTCCCTTGTTTCAATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGAtg                                                                                 <  1:5620188/93‑1 (MQ=255)

          gCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGAtgt                                                                                <  2:1235374/93‑1 (MQ=255)

          gCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGAtgt                                                                                >  1:2411802/1‑93 (MQ=255)

            aaTAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGc                                                                              >  2:6644859/1‑93 (MQ=255)

            aaTAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGc                                                                              >  1:7496580/1‑93 (MQ=255)

            aaTAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGc                                                                              >  2:85500/1‑93 (MQ=255)

             aTAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCa                                                                             <  1:100903/93‑1 (MQ=255)

              tAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCaa                                                                            >  2:6262723/1‑93 (MQ=255)

               aaGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAAc                                                                           <  1:5555540/93‑1 (MQ=255)

                aGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAAcc                                                                          >  2:3255192/1‑93 (MQ=255)

                   ggTTCCCTTGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGc                                                                       <  1:6181964/93‑1 (MQ=255)

                   ggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGc                                                                       <  1:2846820/93‑1 (MQ=255)

                    gTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGct                                                                      >  2:3610806/1‑93 (MQ=255)

                     ttCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGctt                                                                     >  2:2308195/1‑93 (MQ=255)

                      tCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttc                                                                    <  2:1488820/93‑1 (MQ=255)

                      tCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttc                                                                    <  1:6302791/93‑1 (MQ=255)

                      tCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttc                                                                    >  1:5845788/1‑93 (MQ=255)

                       cccATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttct                                                                   <  2:3961921/93‑1 (MQ=255)

                        ccATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttctt                                                                  >  2:1146922/1‑93 (MQ=255)

                        ccATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttctt                                                                  <  2:1884747/93‑1 (MQ=255)

                        ccATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGcttctt                                                                  >  1:4429048/1‑93 (MQ=255)

                          aTGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCa                                                                >  2:3096025/1‑93 (MQ=255)

                          aTGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCa                                                                <  1:4085076/93‑1 (MQ=255)

                          aTGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCa                                                                <  1:168408/93‑1 (MQ=255)

                          aTGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCa                                                                >  2:4570998/1‑93 (MQ=255)

                           tGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCAc                                                               <  1:2155171/93‑1 (MQ=255)

                            gTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACt                                                              <  1:4542605/93‑1 (MQ=255)

                             tttACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACtt                                                             <  2:259380/93‑1 (MQ=255)

                              ttACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTa                                                            <  2:2423794/93‑1 (MQ=255)

                              ttACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTa                                                            >  1:6026483/1‑93 (MQ=255)

                                aCATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATc                                                          <  1:1967124/93‑1 (MQ=255)

                                 cATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCa                                                         <  1:1433991/93‑1 (MQ=255)

                                   tCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGc                                                       >  2:1743754/1‑93 (MQ=255)

                                          cgcAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaa                                                <  1:3126525/93‑1 (MQ=255)

                                          cgcAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaa                                                <  1:7269783/93‑1 (MQ=255)

                                          cgcAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaa                                                <  2:7284238/93‑1 (MQ=255)

                                           gcAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaaa                                               <  1:507290/93‑1 (MQ=255)

                                           gcAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaaa                                               <  1:6111117/93‑1 (MQ=255)

                                             aaGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACg                                             >  2:347877/1‑93 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            >  2:1589763/1‑93 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            >  2:597038/1‑93 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            <  2:1300585/93‑1 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            <  1:6959779/93‑1 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            <  2:4264067/93‑1 (MQ=255)

                                              aGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGa                                            >  1:2459318/1‑93 (MQ=255)

                                               gCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGag                                           <  2:3547546/93‑2 (MQ=255)

                                               gCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGag                                           <  2:5768179/93‑2 (MQ=255)

                                               gCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGag                                           <  2:2920269/93‑2 (MQ=255)

                                                ccAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGagg                                          <  1:2965420/93‑3 (MQ=255)

                                                ccAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGagg                                          >  2:4537509/1‑91 (MQ=255)

                                                 cAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAggg                                         <  2:4905083/93‑1 (MQ=255)

                                                  aaTGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGc                                        >  2:2064717/1‑93 (MQ=255)

                                                  aaTGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGc                                        >  2:2784058/1‑93 (MQ=255)

                                                      ctgcGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                    >  1:5177427/1‑93 (MQ=255)

                                                      ctgcGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                    <  1:2768213/93‑1 (MQ=255)

                                                       tgcGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                   <  1:5139603/93‑1 (MQ=255)

                                                       tgcGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                   >  1:4302/1‑93 (MQ=255)

                                                        gcGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                  >  1:712766/1‑93 (MQ=255)

                                                         cGGTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                 <  1:2259224/93‑1 (MQ=255)

                                                           gTTATCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg                               >  1:1369840/1‑93 (MQ=255)

                                                            ttATGGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                              >  2:7371435/1‑93 (MQ=255)

                                                              aTCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                            <  2:3837924/93‑1 (MQ=255)

                                                              aTCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                            >  2:7202092/1‑93 (MQ=255)

                                                              aTCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                            <  1:4143446/93‑1 (MQ=255)

                                                               tCGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                           <  1:3496715/93‑1 (MQ=255)

                                                                cGAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                          >  2:5872394/1‑93 (MQ=255)

                                                                 gAACAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                         <  1:4238596/93‑1 (MQ=255)

                                                                 gAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                         <  1:5384783/93‑1 (MQ=255)

                                                                  aaCAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACg                        >  1:2097340/1‑93 (MQ=255)

                                                                  aaCAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACg                        >  2:1600263/1‑93 (MQ=255)

                                                                   aCAAATTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                       <  1:6777511/93‑1 (MQ=255)

                                                                     aaaTTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                     >  2:447912/1‑93 (MQ=255)

                                                                     aaaTTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                     >  1:6005187/1‑93 (MQ=255)

                                                                      aaTTCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                    >  2:5470086/1‑93 (MQ=255)

                                                                        ttCGTCGCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgc                  <  1:1735860/93‑1 (MQ=255)

                                                                         tcgtcgCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                 <  1:4606991/93‑1 (MQ=255)

                                                                           gtcgCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc               >  1:2514165/1‑93 (MQ=255)

                                                                            tcgCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc              >  2:7008701/1‑93 (MQ=255)

                                                                            tcgCGCTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc              <  1:3719071/93‑1 (MQ=255)

                                                                                gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg          <  1:4729678/93‑1 (MQ=255)

                                                                                 cTCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc         <  1:3618902/93‑1 (MQ=255)

                                                                                  tCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg        >  1:6866641/1‑93 (MQ=255)

                                                                                  tCGTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg        >  2:1005516/1‑93 (MQ=255)

                                                                                  tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGAGCCGCCAATgcg        >  1:6573614/1‑93 (MQ=255)

                                                                                    gTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt      <  2:6436545/93‑1 (MQ=255)

                                                                                    gTCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt      >  1:3873012/1‑93 (MQ=255)

                                                                                     tCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg     >  1:1162876/1‑93 (MQ=255)

                                                                                     tCCACGCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg     >  1:7002178/1‑93 (MQ=255)

                                                                                     tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg     >  2:5301995/1‑93 (MQ=255)

                                                                                       cacgCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGc   >  2:4916974/1‑93 (MQ=255)

                                                                                        acacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCa  >  2:2036131/1‑93 (MQ=255)

TAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCA  >  NC_000913/547555‑547735

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