BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,566	A→G	G515G (GGT→GGC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,566	0	A	G	97.4%	346.1 / ‑0.5	115	G515G (GGT→GGC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (47/65); total (49/66)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 119 total aligned reads displayed.

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                         <  2:564418/93‑1 (MQ=255)

 ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCa                                                                                        >  1:2236506/1‑93 (MQ=255)

  tCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAt                                                                                       <  2:168477/93‑1 (MQ=255)

   ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                      <  2:2330007/93‑1 (MQ=255)

   ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                      <  2:118722/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                     >  1:2098333/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:3726569/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  1:1366243/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:3726567/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:1285554/1‑93 (MQ=255)

         aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                                                                >  1:2600323/1‑93 (MQ=255)

         aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCAATTACTGAc                                                                                >  1:1594512/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATCGCCCGTTGCATTAAGGGTGGCGCAATTCCTGACAt                                                                              >  1:2947616/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                              >  2:2155465/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                              >  1:2947615/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                              <  2:104890/93‑1 (MQ=255)

            gTCCGTTGTAATATTAACTATCCCCATACGGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                             <  2:2501728/93‑1 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                             <  1:1985362/93‑1 (MQ=255)

            gTCCGTTGAAAGATTAACTACCCCCTGACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                             <  1:1985363/93‑1 (MQ=255)

             tCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                                                                            <  2:4393989/93‑1 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                          <  1:443065/93‑1 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        >  2:2584698/1‑93 (MQ=255)

                  tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                                                                       <  2:4200546/93‑1 (MQ=255)

                  tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                                                                       <  2:4657032/93‑1 (MQ=255)

                    aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                     <  2:3483046/93‑1 (MQ=255)

                      atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                                                                   >  2:285233/1‑93 (MQ=255)

                      atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                                                                   >  1:2432424/1‑93 (MQ=255)

                        atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGCGGCGCCATTACTGACATTTACGATACCGCt                                                                 <  1:2605306/93‑1 (MQ=255)

                         ttAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTa                                                                >  2:1150809/1‑93 (MQ=255)

                           aaCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAcc                                                              >  1:3481357/1‑93 (MQ=255)

                            aCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa                                                             >  2:768685/1‑93 (MQ=255)

                            aCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa                                                             >  2:3559486/1‑93 (MQ=255)

                             cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa                                                            >  2:2914461/1‑93 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGACCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                           >  2:405965/1‑93 (MQ=255)

                                tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt                                                         <  1:3147934/93‑1 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       <  1:3758896/93‑1 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       <  1:4418314/93‑1 (MQ=255)

                                  cccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       <  2:1588466/93‑1 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                      <  1:4094287/93‑1 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                      <  2:1592319/93‑1 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                      <  1:4250926/93‑1 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     <  1:85332/93‑1 (MQ=255)

                                     ttACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                    >  2:534990/1‑93 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   <  2:2023309/93‑1 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   <  2:1687580/93‑1 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   <  2:693126/93‑1 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGAGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   >  1:3715612/1‑93 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                 >  2:3070308/1‑93 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                 <  2:2101743/93‑1 (MQ=255)

                                        ccGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                 >  2:4044077/1‑93 (MQ=255)

                                         cGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt                                                <  1:3264305/93‑1 (MQ=255)

                                          gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                               <  2:4544457/90‑1 (MQ=255)

                                          gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                               >  2:3484533/4‑93 (MQ=255)

                                          gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                               >  1:3718977/4‑93 (MQ=255)

                                          gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                               <  2:3452540/90‑1 (MQ=255)

                                            gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                             >  1:4161114/2‑93 (MQ=255)

                                             ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                            <  2:2571838/93‑1 (MQ=255)

                                             ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                            >  2:870483/1‑93 (MQ=255)

                                              gCATTTCGTCCGATAAATCCATAGCCCGTTGAATTAAGGTTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATa                                           <  2:3231765/93‑1 (MQ=255)

                                               cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                          <  2:1424465/93‑1 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     <  1:1721975/93‑1 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     <  1:3010497/93‑1 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     <  1:3796596/93‑1 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                    >  1:3399202/1‑93 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 >  2:2354141/1‑93 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 <  1:1173318/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 <  1:70761/93‑1 (MQ=255)

                                                         gATAAATCCATAGCCCTTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                <  2:1532196/93‑1 (MQ=255)

                                                         caTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                <  2:2368200/92‑1 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                               <  2:2078747/93‑1 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                               >  1:4506423/1‑93 (MQ=255)

                                                           tAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                              >  2:853556/1‑93 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                             <  2:2108009/93‑1 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                             <  2:3779716/93‑1 (MQ=255)

                                                             aaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTGCGGTACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                            <  1:2515977/93‑1 (MQ=255)

                                                             aaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                            <  1:2515976/93‑1 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           >  1:4310611/1‑93 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           >  2:675285/1‑93 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           >  1:2248013/1‑93 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                          <  2:3988107/93‑1 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                          >  1:2830731/1‑93 (MQ=255)

                                                                ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         <  1:4157059/93‑1 (MQ=255)

                                                                ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         <  2:2379556/93‑1 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                       <  2:3398700/93‑1 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                       <  2:4187408/93‑1 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                       <  2:4187409/93‑1 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATGACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                       >  2:2905558/1‑93 (MQ=255)

                                                                    aGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGatgatg                     >  2:237033/1‑93 (MQ=255)

                                                                        cGTTGAATTAAGGGTGGCGCCATTACTGACATTTATGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                 >  1:526042/1‑93 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACCTTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGCTGTCGca                >  2:277791/1‑93 (MQ=255)

                                                                          ttGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcat               >  1:3494922/1‑93 (MQ=255)

                                                                             aaTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca            <  2:2720433/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa         <  1:1562254/93‑1 (MQ=255)

                                                                                 aaGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg        >  1:2851708/1‑93 (MQ=255)

                                                                                   gggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt      >  1:3691378/1‑93 (MQ=255)

                                                                                     gTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc    <  2:4348819/93‑1 (MQ=255)

                                                                                      tGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa   <  2:4741023/93‑1 (MQ=255)

                                                                                      tGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa   <  2:2564922/93‑1 (MQ=255)

                                                                                       ggCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg  >  1:2652326/1‑93 (MQ=255)

                                                                                       ggCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg  <  1:1206549/93‑1 (MQ=255)

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

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