BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,453,403	T→C	G15G (GGT→GGC)	yjfG →	putative ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,453,403	0	T	C	99.7%	896.5 / ‑5.0	287	G15G (GGT→GGC)	yjfG	putative ligase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (132/154); total (132/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 294 total aligned reads displayed.

GACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCA  >  NC_000913/4453313‑4453492

gACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCt                                                                                         <  1:6664806/93‑1 (MQ=255)

gACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCt                                                                                         >  1:3597395/1‑93 (MQ=255)

 aCATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTg                                                                                        >  1:5195117/1‑93 (MQ=255)

  cATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTgg                                                                                       >  1:6198204/1‑93 (MQ=255)

  cATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTgg                                                                                       <  2:4434784/93‑1 (MQ=255)

     atGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGa                                                                                    <  2:5688506/93‑1 (MQ=255)

      tGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGAt                                                                                   >  2:6493285/1‑93 (MQ=255)

      tGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGAt                                                                                   >  2:119499/1‑93 (MQ=255)

       gCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATg                                                                                  <  2:1567338/93‑1 (MQ=255)

        cGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGc                                                                                 <  1:1440783/93‑1 (MQ=255)

          tAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTg                                                                               <  2:5557895/93‑1 (MQ=255)

           aaaaTGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTgg                                                                              <  1:1073270/93‑1 (MQ=255)

           aaaaTGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTgg                                                                              <  1:1796112/93‑1 (MQ=255)

           aaaaTGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTgg                                                                              <  2:2658136/93‑1 (MQ=255)

              aTGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGgcgc                                                                           <  1:4567993/93‑1 (MQ=255)

              aTGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGgcgc                                                                           <  1:3257214/93‑1 (MQ=255)

                gAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGgcgcgc                                                                         >  1:2837429/1‑93 (MQ=255)

                 aGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGcc                                                                        >  1:5530709/1‑93 (MQ=255)

                 aGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGcc                                                                        >  2:4986143/1‑93 (MQ=255)

                  gCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCa                                                                       >  2:4342933/1‑93 (MQ=255)

                    ggCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGt                                                                     <  2:3307998/93‑1 (MQ=255)

                      cAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTa                                                                   <  2:3830365/93‑1 (MQ=255)

                      cAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTa                                                                   <  2:1632240/93‑1 (MQ=255)

                       aGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAg                                                                  <  2:4995025/93‑1 (MQ=255)

                        gATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAgg                                                                 >  1:4846900/1‑93 (MQ=255)

                        gATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAgg                                                                 >  2:4690667/1‑93 (MQ=255)

                          ttAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGcc                                                               <  2:1488715/93‑1 (MQ=255)

                          ttAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGcc                                                               >  1:4314031/1‑93 (MQ=255)

                             aaaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATg                                                            <  2:426555/93‑1 (MQ=255)

                              aaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGGCATGa                                                           >  1:3521396/1‑93 (MQ=255)

                              aaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGa                                                           <  1:2730361/93‑1 (MQ=255)

                              aaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGa                                                           >  1:6030674/1‑93 (MQ=255)

                               aaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGaa                                                          <  1:231191/93‑1 (MQ=255)

                               aaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGaa                                                          <  1:3690667/93‑1 (MQ=255)

                               aaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGaa                                                          <  2:819613/93‑1 (MQ=255)

                                aaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAg                                                         <  2:1332670/93‑1 (MQ=255)

                                aaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAg                                                         >  1:3349526/1‑93 (MQ=255)

                                aaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAg                                                         >  2:2082677/1‑93 (MQ=255)

                                  ggATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTa                                                       >  2:1738922/1‑93 (MQ=255)

                                  ggATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTa                                                       >  1:3234017/1‑93 (MQ=255)

                                   gATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTaa                                                      >  1:697005/1‑93 (MQ=255)

                                   gATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGCTAGGCCATGAAGTaa                                                      >  1:6380670/1‑93 (MQ=255)

                                    aTAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAAc                                                     <  1:2691571/93‑1 (MQ=255)

                                     tAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACg                                                    >  1:155146/1‑93 (MQ=255)

                                      aGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACgg                                                   >  1:2962550/1‑93 (MQ=255)

                                        tGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGt                                                 >  2:2890756/1‑93 (MQ=255)

                                        tGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGt                                                 >  2:1929546/1‑93 (MQ=255)

                                          aCGTATGCGCAGTCATATTTTAGGAATTGGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTc                                               <  1:511095/93‑1 (MQ=255)

                                           cGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGCTCg                                              >  1:6268114/1‑93 (MQ=255)

                                             tATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGa                                            >  1:3403237/1‑93 (MQ=255)

                                               tGCGCATTCATATTTTAGGAATTTGTGGCTCATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACg                                          <  1:1024176/93‑1 (MQ=255)

                                               tGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACg                                          >  2:1871419/1‑93 (MQ=255)

                                                gcgcATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGc                                         >  2:3182615/1‑93 (MQ=255)

                                                  gcATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCa                                       <  2:5750582/93‑1 (MQ=255)

                                                   cATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCaa                                      <  2:5418359/93‑1 (MQ=255)

                                                    aTTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAt                                     >  1:5632010/1‑93 (MQ=255)

                                                    aTTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAt                                     >  1:6601281/1‑93 (MQ=255)

                                                    aTTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAt                                     >  2:1016434/1‑93 (MQ=255)

                                                     ttCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAtg                                    <  2:2978624/93‑1 (MQ=255)

                                                     ttCATATTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAtg                                    <  2:2824776/93‑1 (MQ=255)

                                                     ttCATATTTTAGGAATTTGTGGCACACTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAtg                                    >  2:2229732/1‑93 (MQ=255)

                                                       cATATTTTAGGAATTTGTGGCACAGTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAAtgtg                                  >  2:4861862/1‑93 (MQ=255)

                                                         tatTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTa                                <  1:2076002/93‑1 (MQ=255)

                                                         tatTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTa                                >  2:4103009/1‑93 (MQ=255)

                                                         tatTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTa                                <  1:857852/93‑1 (MQ=255)

                                                         tatTTTAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTa                                >  2:6271067/1‑93 (MQ=255)

                                                           ttttAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATc                              <  1:24754/93‑1 (MQ=255)

                                                             ttAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccg                            <  1:3087649/93‑1 (MQ=255)

                                                              tAGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccgc                           >  2:6321389/1‑93 (MQ=255)

                                                               aGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccgtc                          <  1:4206020/93‑3 (MQ=255)

                                                               aGGAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccgcc                          >  1:3845492/1‑93 (MQ=255)

                                                                ggAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccgccg                         >  1:1202230/1‑93 (MQ=255)

                                                                 gAATTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGa                        >  2:6008982/1‑93 (MQ=255)

                                                                   aTTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATg                      >  2:5279802/1‑93 (MQ=255)

                                                                   aTTTGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATg                      >  1:3938021/1‑93 (MQ=255)

                                                                    tttGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGa                     <  1:4968340/93‑1 (MQ=255)

                                                                    tttGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGa                     <  1:399877/93‑1 (MQ=255)

                                                                    tttGTGGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGa                     <  2:92470/93‑1 (MQ=255)

                                                                       gtgGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCa                  >  1:4118220/1‑93 (MQ=255)

                                                                         ggCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCAcc                <  2:51315/93‑1 (MQ=255)

                                                                         ggCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCAcc                <  1:4344717/93‑1 (MQ=255)

                                                                         ggCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCAcc                >  2:5353026/1‑93 (MQ=255)

                                                                          gCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCt               <  1:5456584/93‑1 (MQ=255)

                                                                           cacaTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCtt              <  2:4783856/93‑1 (MQ=255)

                                                                            acaTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTa             <  2:4405947/93‑1 (MQ=255)

                                                                            acaTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTa             >  2:3982628/1‑93 (MQ=255)

                                                                             caTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTAc            <  1:2181016/93‑1 (MQ=255)

                                                                             caTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTAc            <  2:2297752/93‑1 (MQ=255)

                                                                              aTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACt           <  1:1611149/93‑1 (MQ=255)

                                                                              aTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACt           <  1:770428/93‑1 (MQ=255)

                                                                              aTTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACt           <  2:2386452/93‑1 (MQ=255)

                                                                               tttATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACtt          >  2:3948544/1‑93 (MQ=255)

                                                                               tttATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACtt          <  2:674912/93‑1 (MQ=255)

                                                                                ttATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTg         <  1:3174704/93‑1 (MQ=255)

                                                                                ttATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTg         <  1:2319714/93‑1 (MQ=255)

                                                                                 tATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTga        <  1:1882144/93‑1 (MQ=255)

                                                                                    gggCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTACCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGaa     <  1:1928550/93‑1 (MQ=255)

                                                                                    gggCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGaa     <  1:167878/93‑1 (MQ=255)

                                                                                      gCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGc   <  1:5787129/93‑1 (MQ=255)

                                                                                       cGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCa  <  2:5446303/93‑1 (MQ=255)

GACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCA  >  NC_000913/4453313‑4453492

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