BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,545,868	A→G	K5K (AAA→AAG)	fimG →	FimG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,545,868	0	A	G	94.0%	744.9 / 27.0	234	K5K (AAA→AAG)	fimG	FimG
Reads supporting (aligned to +/- strand): ref base A (10/4); new base G (117/102); total (127/107)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 238 total aligned reads displayed.

ACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGC  >  NC_000913/4545778‑4545960

acaGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCg                                                                                            <  2:1426278/93‑1 (MQ=255)

  aGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTg                                                                                          >  1:4114496/1‑93 (MQ=255)

    gTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTggg                                                                                        <  1:554287/93‑1 (MQ=255)

      gCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgta                                                                                      <  1:874132/93‑1 (MQ=255)

      gCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgta                                                                                      <  1:5724788/93‑1 (MQ=255)

      gCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgta                                                                                      >  1:3224041/1‑93 (MQ=255)

       ccTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtat                                                                                     <  2:4687463/93‑1 (MQ=255)

        cTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatg                                                                                    <  2:4785183/93‑1 (MQ=255)

        cTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatg                                                                                    >  1:3142449/1‑93 (MQ=255)

         tGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgt                                                                                   >  2:5537239/1‑93 (MQ=255)

         tGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgt                                                                                   <  2:3432669/93‑1 (MQ=255)

         tGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgt                                                                                   >  1:5950343/1‑93 (MQ=255)

          gTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgta                                                                                  <  2:4903998/93‑1 (MQ=255)

          gTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgta                                                                                  >  1:879035/1‑93 (MQ=255)

          gTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGgtatgta                                                                                  >  1:760299/1‑93 (MQ=255)

          gTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGgtatgta                                                                                  >  1:1048147/1‑93 (MQ=255)

             aCTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTg                                                                               >  2:1538262/1‑93 (MQ=255)

             aCTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTg                                                                               >  2:1253447/1‑93 (MQ=255)

             aCTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTg                                                                               >  2:5221225/1‑93 (MQ=255)

             aCTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTg                                                                               >  2:4547400/1‑93 (MQ=255)

                 cGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTggcgg                                                                           <  1:1695060/93‑1 (MQ=255)

                    ggCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCaa                                                                        <  1:2765261/93‑1 (MQ=255)

                    ggCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCaa                                                                        >  2:1718315/1‑93 (MQ=255)

                    ggCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCaa                                                                        <  1:2648745/93‑1 (MQ=255)

                     gCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAat                                                                       <  1:5033796/93‑1 (MQ=255)

                      cATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTAGCGGCAata                                                                      >  1:3733684/1‑93 (MQ=255)

                       atatCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAatat                                                                     <  2:833743/93‑1 (MQ=255)

                       atatCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAatat                                                                     <  1:497570/93‑1 (MQ=255)

                        tatCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCATGCGTGGGTATGTATTGGCGGCAATAtt                                                                    >  1:3748951/1‑93 (MQ=255)

                        tatCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATAtt                                                                    <  2:5997474/93‑1 (MQ=255)

                         atCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTg                                                                   >  1:5471051/1‑93 (MQ=255)

                         atCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATGTGCGGCCATATTg                                                                   >  1:292230/1‑93 (MQ=255)

                          tCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTgg                                                                  <  1:5157580/93‑1 (MQ=255)

                           cAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGgc                                                                 <  1:3505784/93‑1 (MQ=255)

                              tGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTTgcacc                                                              >  1:2153401/1‑90 (MQ=255)

                              tGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCt                                                              <  2:5186643/93‑1 (MQ=255)

                               gCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTc                                                             <  1:1205007/93‑1 (MQ=255)

                               gCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGGGCAAGGGTGGGAAGGTATTGGGGGAAAAATTGGGGCtg                                                             >  2:5749264/1‑92 (MQ=255)

                                 cACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGc                                                           <  1:4815668/93‑1 (MQ=255)

                                  aCGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCa                                                          >  1:90119/1‑93 (MQ=255)

                                   cGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCaa                                                         <  1:1000695/93‑1 (MQ=255)

                                   cGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCaa                                                         >  2:5947710/1‑93 (MQ=255)

                                    ggCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAg                                                        <  1:4523878/93‑1 (MQ=255)

                                    ggCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAg                                                        >  1:1967724/1‑93 (MQ=255)

                                    ggCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAg                                                        >  2:249306/1‑93 (MQ=255)

                                     gCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAgg                                                       >  1:5870964/1‑92 (MQ=255)

                                     gCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGt                                                       >  1:3020077/1‑93 (MQ=255)

                                      cTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTg                                                      <  2:2695782/93‑1 (MQ=255)

                                      cTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTg                                                      >  2:3593264/1‑93 (MQ=255)

                                       tACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGc                                                     <  1:3968161/93‑1 (MQ=255)

                                       tACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGc                                                     >  2:625989/1‑93 (MQ=255)

                                        aCCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcg                                                    >  2:3104673/1‑93 (MQ=255)

                                         ccTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcga                                                   <  1:1947661/93‑1 (MQ=255)

                                         ccTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcga                                                   <  2:1848475/93‑1 (MQ=255)

                                            tCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcgacga                                                <  1:2538337/93‑1 (MQ=255)

                                              aCTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATa                                              >  2:4052851/1‑93 (MQ=255)

                                               ctctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATAc                                             >  2:4867126/1‑93 (MQ=255)

                                               ctctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATAc                                             <  1:1613829/93‑1 (MQ=255)

                                               ctctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATAc                                             >  1:2629797/1‑93 (MQ=255)

                                               ctctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATAc                                             >  1:3932336/1‑93 (MQ=255)

                                                tctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACa                                            <  2:4475876/93‑1 (MQ=255)

                                                tctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACa                                            <  2:413627/93‑1 (MQ=255)

                                                tctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACa                                            <  1:5793847/93‑1 (MQ=255)

                                                   tGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGc                                         <  1:4443933/93‑1 (MQ=255)

                                                   tGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGc                                         <  2:2554651/93‑1 (MQ=255)

                                                     aaTATCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAg                                       >  2:1660089/1‑93 (MQ=255)

                                                         tCAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGa                                   <  2:1635070/93‑1 (MQ=255)

                                                          cAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGAt                                  <  1:536436/93‑1 (MQ=255)

                                                          cAGTAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGAt                                  >  1:83510/1‑93 (MQ=255)

                                                             tAACTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACGGGCAGCCGATGTc                               >  2:2209202/1‑93 (MQ=255)

                                                               aCTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCAc                             <  1:1476806/93‑1 (MQ=255)

                                                               aCTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCAc                             >  2:595209/1‑93 (MQ=255)

                                                                cTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCAcc                            >  1:3098932/1‑93 (MQ=255)

                                                                cTGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCAcc                            <  1:1833437/93‑1 (MQ=255)

                                                                 tGGAGATGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCa                           <  2:4087185/93‑1 (MQ=255)

                                                                   gagaTGCTCATGAAATGGTGCAATCGTGGGTATGTATTGGCGGCAATAGTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATc                         <  2:3239613/93‑1 (MQ=255)

                                                                   gagaTGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATc                         >  2:240399/1‑93 (MQ=255)

                                                                   gagaTGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATc                         >  1:2558444/1‑93 (MQ=255)

                                                                    agaTGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCa                        <  2:3758144/93‑1 (MQ=255)

                                                                    agaTGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCa                        <  2:1042830/93‑1 (MQ=255)

                                                                      aTGCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACg                      >  1:462310/1‑93 (MQ=255)

                                                                      aTGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACg                      >  1:4350176/1‑93 (MQ=255)

                                                                        gCTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGt                    >  1:5702543/1‑93 (MQ=255)

                                                                         cTCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTg                   >  1:244459/1‑93 (MQ=255)

                                                                          tCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGa                  >  1:2865017/1‑93 (MQ=255)

                                                                          tCATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGa                  <  2:5584627/93‑1 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGaa                 <  2:4264768/93‑1 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGaa                 >  2:5396524/1‑93 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGaa                 >  2:696368/1‑93 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGaa                 >  2:906026/1‑93 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGGGaa                 >  1:141070/1‑93 (MQ=255)

                                                                           cATGAAATGGTGCAAGCGGGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGGGCGACGATACAGGCAGCCGATGTCACCATCACGGTGaa                 <  2:3924146/93‑1 (MQ=255)

                                                                            aTGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAAc                <  2:4158826/93‑1 (MQ=255)

                                                                            aTGAAATGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAAc                >  2:2161239/1‑93 (MQ=255)

                                                                                aaTGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTa            >  1:3316622/1‑93 (MQ=255)

                                                                                 aTGGTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTaa           >  2:2067714/1‑93 (MQ=255)

                                                                                 aTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTaa           <  1:4223163/93‑1 (MQ=255)

                                                                                   ggTGCAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAgg         <  2:2441934/93‑1 (MQ=255)

                                                                                       cAAGCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGgtcgt     <  1:4091730/93‑1 (MQ=255)

                                                                                          aCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGc  >  2:3009607/1‑93 (MQ=255)

ACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGC  >  NC_000913/4545778‑4545960

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