BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,331	A→G	S404S (TCA→TCG)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,331	0	A	G	99.4%	533.0 / ‑5.6	175	S404S (TCA→TCG)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (60/114); total (60/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.66e-01

Only 100 of 197 total aligned reads displayed.

GCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTG  >  NC_000913/2208241‑2208420

 ccGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCaa                                                                                        <  2:786041/93‑1 (MQ=255)

 ccGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCaa                                                                                        <  2:4793090/93‑1 (MQ=255)

   gctaaGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagc                                                                                      <  2:3834926/90‑1 (MQ=255)

   gcgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCCAAATGTTTTGCAagg                                                                                      >  2:7428636/1‑92 (MQ=255)

    cgAAGTGAACAGGAAAGCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagca                                                                                     <  2:6327061/93‑1 (MQ=255)

    cgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagca                                                                                     <  2:135204/93‑1 (MQ=255)

     gAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTTGCAagcag                                                                                    >  1:7705176/1‑93 (MQ=255)

       aGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGGAAATTGTTTTGCAAgcagca                                                                                  >  1:4216474/1‑93 (MQ=255)

       aGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAgcagca                                                                                  <  2:3641363/93‑1 (MQ=255)

        gTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCaa                                                                                 <  2:3151288/93‑1 (MQ=255)

        gTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCaa                                                                                 <  1:2854846/93‑1 (MQ=255)

           aaCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAAc                                                                              <  1:1190859/93‑1 (MQ=255)

           aaCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAAc                                                                              <  1:5509386/93‑1 (MQ=255)

           aaCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAAc                                                                              <  1:3330335/93‑1 (MQ=255)

            aCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACg                                                                             <  2:3274526/93‑1 (MQ=255)

             cAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGc                                                                            <  1:1495184/93‑1 (MQ=255)

             cAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGc                                                                            <  1:505685/93‑1 (MQ=255)

              aGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCa                                                                           <  2:4780715/93‑1 (MQ=255)

              aGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCa                                                                           <  2:2952064/93‑1 (MQ=255)

               ggAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTTGCAAGCAGCAAAAAGCAg                                                                          >  2:7274232/1‑93 (MQ=255)

                gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGt                                                                         <  1:5714600/93‑1 (MQ=255)

                gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGt                                                                         <  2:2748766/93‑1 (MQ=255)

                   aaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCCGGCCAATTGTTTTGCAAGCAGCAAAAAGCAGTgtc                                                                      >  1:1699511/1‑91 (MQ=255)

                     cTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGa                                                                    <  1:5556723/93‑1 (MQ=255)

                     cTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGa                                                                    <  1:1733058/93‑1 (MQ=255)

                     cTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGa                                                                    <  1:252546/93‑1 (MQ=255)

                     cTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGa                                                                    <  1:6774928/93‑1 (MQ=255)

                       cACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATg                                                                  <  2:1783750/93‑1 (MQ=255)

                        cccGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATgg                                                                 <  1:3101974/92‑1 (MQ=255)

                         ccGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGc                                                                <  1:4929930/93‑1 (MQ=255)

                           ggCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGa                                                              <  2:7602537/93‑1 (MQ=255)

                              aGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCCAAACGCAGCGCCGCTGGCCAgct                                                           >  2:5003537/1‑93 (MQ=255)

                              aGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgct                                                           <  1:6515160/93‑1 (MQ=255)

                               gCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgctg                                                          >  2:6704204/1‑93 (MQ=255)

                                cTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgctgc                                                         <  1:7586685/93‑1 (MQ=255)

                                  aCGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCCATGGCGAGctgctg                                                       >  2:1644886/1‑93 (MQ=255)

                                     cgTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTa                                                    <  2:5593274/93‑1 (MQ=255)

                                     cgTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTa                                                    <  2:5701999/93‑1 (MQ=255)

                                      gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTTGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAg                                                   >  1:3693139/1‑93 (MQ=255)

                                      gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAg                                                   >  2:2654036/1‑93 (MQ=255)

                                      gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAg                                                   >  2:5885055/1‑93 (MQ=255)

                                       tCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGc                                                  <  1:7792840/93‑1 (MQ=255)

                                        cGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAAAGCAGTGCCGATGGCGAGCTGCTGTTAGCt                                                 >  2:3113270/1‑93 (MQ=255)

                                         ggCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTa                                                <  2:5751048/93‑1 (MQ=255)

                                          gCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAc                                               >  1:3510276/1‑93 (MQ=255)

                                          gCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAc                                               <  1:2117961/93‑1 (MQ=255)

                                          gCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAc                                               >  1:5469359/1‑93 (MQ=255)

                                          gCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAc                                               <  2:290283/93‑1 (MQ=255)

                                           cTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGGAGTGCCGATGGCGAGCCGCTGTTAGCTAcc                                              >  1:4647864/1‑93 (MQ=255)

                                           cTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAcc                                              >  1:3932540/1‑93 (MQ=255)

                                            tGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAccc                                             <  1:3279368/93‑1 (MQ=255)

                                            tGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAccc                                             >  1:3977436/1‑93 (MQ=255)

                                             gTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGATGTTAGCTACCCg                                            >  1:7866077/1‑93 (MQ=255)

                                              tttAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAAATGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGc                                           >  2:861004/1‑93 (MQ=255)

                                              tttAGTTTTCAGGCTGGCGCACTAGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGc                                           <  2:5651895/93‑1 (MQ=255)

                                                tAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCaa                                         <  2:1863041/93‑1 (MQ=255)

                                                tAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCaa                                         <  2:3510170/93‑1 (MQ=255)

                                                tAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCaa                                         >  2:317708/1‑93 (MQ=255)

                                                  gTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAAcc                                       >  2:4194424/1‑93 (MQ=255)

                                                    tttCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGc                                     <  1:6827697/93‑1 (MQ=255)

                                                       cAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTa                                  <  2:3383600/93‑1 (MQ=255)

                                                        aGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTaa                                 >  2:7120909/1‑93 (MQ=255)

                                                           cTGGCGCACTGGCGGGCGGGAAAATTGTTTCGCAAGCAGCAAATCGCAGTGCCGATGGCGCGCTGCTGTTAGCTACCCGCAACCGCTTAagca                              <  2:4604131/93‑1 (MQ=255)

                                                            gggCGCACTGGGGGGCGGGAAAATGGTTTCGCAAGCAGCAAACCGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAGCCGCTTAagcag                             <  2:7363548/92‑1 (MQ=255)

                                                             ggCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagcagc                            <  2:7869578/93‑1 (MQ=255)

                                                             ggCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagcagc                            <  2:4130289/93‑1 (MQ=255)

                                                              gcgcACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCg                           >  1:2701266/1‑93 (MQ=255)

                                                               ggcACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGt                          <  1:1569286/92‑1 (MQ=255)

                                                                gcACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGtt                         <  1:537983/93‑1 (MQ=255)

                                                                 cACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCGGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtg                        >  2:4670984/1‑93 (MQ=255)

                                                                 cACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtg                        >  2:6225548/1‑93 (MQ=255)

                                                                 cACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtg                        >  1:2398267/1‑93 (MQ=255)

                                                                  aCTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtgt                       >  1:5901463/1‑93 (MQ=255)

                                                                   cTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGAGAGCTGCTGTTTGCTACCCGCAACAGCTTAAGCAGCGCtgtg                      >  2:311187/1‑93 (MQ=255)

                                                                    tGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGc                     <  2:5362259/93‑1 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  >  2:3666609/1‑93 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  <  1:7241231/93‑1 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  <  2:3542687/93‑1 (MQ=255)

                                                                        gggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCt                 >  2:2577283/1‑93 (MQ=255)

                                                                         ggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTg                >  1:7467837/1‑93 (MQ=255)

                                                                         ggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTg                <  2:3038460/93‑1 (MQ=255)

                                                                              gcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGcc           >  1:2916570/1‑93 (MQ=255)

                                                                               cAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCt          <  2:5037587/93‑1 (MQ=255)

                                                                               cAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCt          <  1:4768560/93‑1 (MQ=255)

                                                                                aaaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTg         >  1:4953723/1‑93 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        <  1:6099385/93‑1 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        <  1:4295615/93‑1 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        <  2:1146294/93‑1 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        <  2:1535491/93‑1 (MQ=255)

                                                                                  aTTGTTTCGCAAGCAGCAAAATGCATTGCCGATGGCGAGCTGCTGGTAACTACCCGCAACCGCTTAAGCAGCGGTGTGCCGCTGTCGCCTGAt       >  2:7048071/1‑93 (MQ=255)

                                                                                   ttGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATg      >  1:3019049/1‑93 (MQ=255)

                                                                                   ttGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATg      >  1:6300759/1‑93 (MQ=255)

                                                                                    tGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGc     >  1:2367804/1‑93 (MQ=255)

                                                                                    tGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGc     >  1:6748057/1‑93 (MQ=255)

                                                                                     gTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGcc    <  2:5786909/93‑1 (MQ=255)

                                                                                     gTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGcc    <  2:4885754/93‑1 (MQ=255)

                                                                                     gTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGcc    <  1:6592259/93‑1 (MQ=255)

                                                                                      tttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCt   >  1:1002487/1‑93 (MQ=255)

                                                                                       ttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTg  <  2:6039319/93‑1 (MQ=255)

                                                                                       ttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTg  <  2:5370580/93‑1 (MQ=255)

GCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTG  >  NC_000913/2208241‑2208420

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