BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,684	+A	coding (400/501 nt)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,681	1	.	A	99.2%	1401.2 / 1.7	392	coding (397/501 nt)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base A (179/210); total (180/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑C‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

gAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c                                                                                           <  1:2492341/93‑1 (MQ=255)

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c‑a                                                                                         >  1:6766349/1‑93 (MQ=255)

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c‑a                                                                                         >  1:4642332/1‑93 (MQ=255)

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c‑a                                                                                         >  1:3207653/1‑93 (MQ=255)

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c‑a                                                                                         >  2:8803877/1‑93 (MQ=255)

     aGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaa‑c‑a‑aaac                                                                                    <  1:2796584/93‑3 (MQ=255)

          aTTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACg                                                                                <  2:1140285/93‑1 (MQ=255)

           tttGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGa                                                                               <  2:4063744/93‑1 (MQ=255)

            ttGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGAt                                                                              <  2:3669303/93‑1 (MQ=255)

             tGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGAtt                                                                             >  2:6907713/1‑93 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTga                                                                           >  1:1729788/1‑92 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTa                                                                          >  2:5165303/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTa                                                                          >  1:3385599/1‑93 (MQ=255)

                 tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTAt                                                                         >  1:5088246/1‑93 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCt                                                                     <  1:8071135/93‑1 (MQ=255)

                      agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTg                                                                    <  1:6382892/93‑1 (MQ=255)

                       gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTgg                                                                   <  1:7777813/93‑1 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGa                                                                 >  1:7916393/1‑93 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGa                                                                 >  1:3908839/1‑93 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGa                                                                 <  1:7399783/93‑1 (MQ=255)

                          aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGat                                                                >  2:2824316/1‑93 (MQ=255)

                            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGatat                                                              <  1:137931/93‑1 (MQ=255)

                             cccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATAtt                                                             >  2:6124370/1‑93 (MQ=255)

                              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTc                                                            >  2:737923/1‑93 (MQ=255)

                              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑C‑AGAACAGACGATTGTTTGCCTGGattattt                                                            >  1:4437758/1‑92 (MQ=255)

                               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCt                                                           >  2:2717602/1‑93 (MQ=255)

                               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCt                                                           >  2:1274899/1‑93 (MQ=255)

                                gtAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTa                                                          >  2:3475929/2‑93 (MQ=255)

                                aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTa                                                          >  2:904933/1‑93 (MQ=255)

                                aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTa                                                          <  2:6456480/93‑1 (MQ=255)

                                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACa                                                        <  2:4863509/93‑1 (MQ=255)

                                    cAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGa                                                      <  1:6546895/93‑1 (MQ=255)

                                     aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATACTACAGaa                                                     >  1:3738886/1‑93 (MQ=255)

                                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaa                                                   <  2:7783517/93‑1 (MQ=255)

                                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaa                                                   <  2:4807706/93‑1 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                 <  2:7912458/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                >  2:3133847/1‑93 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                >  1:8777420/1‑93 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  1:2360797/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               >  2:5467404/1‑93 (MQ=255)

                                           agAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  1:5873084/92‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:3568793/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:1023889/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:2425854/93‑1 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑TAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:4365159/93‑1 (MQ=255)

                                               aataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAATA‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:1279587/92‑1 (MQ=255)

                                                 taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                         >  2:5988488/1‑93 (MQ=255)

                                                 taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                         >  2:1882791/1‑93 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       <  2:6307042/93‑1 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACGGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      <  2:1633809/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    <  2:1425384/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    <  1:7321381/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    <  2:4274433/93‑1 (MQ=255)

                                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                   >  1:5175784/1‑93 (MQ=255)

                                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                   >  1:1603956/1‑93 (MQ=255)

                                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                   >  1:303229/1‑93 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  >  1:7676220/1‑93 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  >  2:1017668/1‑93 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                 >  1:5532204/1‑93 (MQ=255)

                                                          tatGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                <  2:5835694/92‑1 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  2:7581475/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaag                              >  1:6581022/1‑92 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  1:6190479/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGACTGACATACCGaaaa                              >  1:5134921/1‑93 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             <  1:5416098/93‑1 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             <  2:6795706/93‑1 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             >  1:6059660/1‑93 (MQ=255)

                                                               aTTTCAAAAACCGTTAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                           >  1:8565156/1‑93 (MQ=255)

                                                               aTTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAggt                           >  1:5338279/1‑92 (MQ=255)

                                                               aTTTCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                           <  2:3261379/93‑1 (MQ=255)

                                                               aTTGCAAAAACCGTGAAATTCCGTAAGAT‑CAA‑AACAGACGATTGTATGCCTGGGATATTCAACAGAAAAAACAGTCTGACATTCCGAAAAGGa                           <  1:7844585/93‑1 (MQ=255)

                                                                tttCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt                          <  1:8355378/93‑1 (MQ=255)

                                                                  tCAAAAACCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa                        <  1:6118513/93‑1 (MQ=255)

                                                                    aaaatCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGAGATATTCTACAGAAAAACCAGTCTGGCATACCGAAAAGGATTaaa                      <  1:1213817/93‑1 (MQ=255)

                                                                    aaaaaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                      <  2:8350712/93‑1 (MQ=255)

                                                                    aaaaaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                      >  2:2571990/1‑93 (MQ=255)

                                                                    aaaaaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                      <  1:4776098/93‑1 (MQ=255)

                                                                     aaaaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                     >  2:5531131/1‑93 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  2:6628342/93‑1 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  1:5679409/93‑1 (MQ=255)

                                                                        aCCGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGtt                  <  1:2605288/93‑1 (MQ=255)

                                                                         ccGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTta                 <  2:2789461/93‑1 (MQ=255)

                                                                         ccGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTta                 >  1:4132224/1‑93 (MQ=255)

                                                                          cGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                <  1:5251398/93‑1 (MQ=255)

                                                                          cGTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                <  1:4986351/93‑1 (MQ=255)

                                                                           gTGAATTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata               <  1:8700821/93‑1 (MQ=255)

                                                                              aaTTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGTAAAGGATTAAAGGTTATATTc            >  1:3606082/1‑93 (MQ=255)

                                                                              aaTTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTc            >  1:667519/1‑93 (MQ=255)

                                                                              aaTTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTc            >  1:1856973/1‑93 (MQ=255)

                                                                              aaTTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTc            <  2:7975010/93‑1 (MQ=255)

                                                                               aTTTCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCa           <  2:2252431/93‑1 (MQ=255)

                                                                                tttCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          >  2:4478690/1‑93 (MQ=255)

                                                                                  tCCGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  2:6864963/93‑1 (MQ=255)

                                                                                    cGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAgg      <  1:8419614/93‑1 (MQ=255)

                                                                                    cGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAgg      >  2:8256769/1‑93 (MQ=255)

                                                                                    cGTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAgg      <  2:8582869/93‑1 (MQ=255)

                                                                                     gTTAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc     >  2:2976723/1‑93 (MQ=255)

                                                                                      ttAGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg    >  2:7684429/1‑93 (MQ=255)

                                                                                        aGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  >  2:2643108/1‑93 (MQ=255)

                                                                                        aGAA‑CAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  <  1:7036346/93‑1 (MQ=255)

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAA‑C‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

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